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This review presents an overview of the thermal properties of mesoscopic structures. The discus- 
sion is based on the concept of electron energy distribution, and, in particular, on controlling and 
probing it. The temperature of an electron gas is determined by this distribution: refrigeration 
is equivalent to narrowing it, and thermometry is probing its convolution with a function char- 
acterizing the measuring device. Temperature exists, strictly speaking, only in quasiequilibrium 
in which the distribution follows the Fermi-Dirac form. Interesting nonequilibrium deviations can 
occur due to slow relaxation rates of the electrons, e.g., among themselves or with lattice phonons. 
Observation and applications of nonequilibrium phenomena are also discussed. The focus in this 
paper is at low temperatures, primarily below 4 K, where physical phenomena on mesoscopic 
scales and hybrid combinations of various types of materials, e.g., superconductors, normal met- 
als, insulators, and doped semiconductors, open up a rich variety of device concepts. This review 
starts with an introduction to theoretical concepts and experimental results on thermal properties 
of mesoscopic structures. Then thermometry and refrigeration are examined with an emphasis on 
experiments. An immediate application of solid-state refrigeration and thermometry is in ultra- 
sensitive radiation detection, which is discussed in depth. This review concludes with a summary 
of pertinent fabrication methods of presented devices. 
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I. INTRODUCTION 

Solid state mesoscopic electronic systems provide a 
micro-laboratory to realize experiments on low tem- 
perature physics, to study quantum phenomena such 
as fundamental relaxation mechanisms in solids, and a 
way to create advanced cryogenic devices. In a broad 
sense, mesoscopic here refers to micro- and nanostruc- 
tures, whose size falls in between atomic and macro- 
scopic scales. The central concept of this Review is 
the energy distribution of mesoscopic electron systems, 
which in thermal equilibrium (Fermi-Dirac distribution) 
determines the temperature of the electron gas. The 
non-Fermi distributions are discussed in depth, since 
they are often encountered and utilized in mesoscopic 
structures and devices. This Review aims to discuss 
the progress mainly during the past decade on how 
electron distributions can be controlled, measured and 
made use of in various device concepts. When appro- 
priate, earlier developments are reviewed as well. The 
central devices and concepts to be discussed are elec- 
tronic refrigerators, thermometers, radiation detectors, 
and distribution-controlled transistors. Typically the 
working principles or resolution of these detectors rely 
on phenomena that show up only at cryogenic tempera- 
tures, i.e., at temperatures of the order of a few kelvin 
and below. A practical threshold in terms of temperature 
is set by liquefaction of helium at 4.2 K. This also sets 
the emphasis in this Review: the devices and principles 
working mostly at temperatures above 4.2 K are at times 
mentioned only for reference. 

Section ^] of this Review introduces formally the cen- 
tral ingredients; the relevant theoretical results are ei- 
ther derived or given there. We also review some of the 
new developments concerning the thermoelectric effects 
in mesoscopic systems. Although the theoretical analysis 
of the effects in the later sections is based on Sec.m the 
main messages can be understood without reading it in 
detail. Section IIIII explains how the electronic temper- 
ature is typically measured and what is required of an 
electronic thermometer. Accurate and fast thermome- 
ters can be utilized for thermal radiation detection as 
explained in Sec. IIVI which reviews such detectors. The 
resolution of these devices is ultimately limited by the 
thermal noise, which can be lowered by refrigeration. In 
Sec. we show how the electron temperature can be 
lowered via electronic means, and discuss the direct ap- 
plications of this refrigeration. Section IVII explains the 
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FIG. 1 Schematic picture of the system considered in this 
review. We describe an electron system in a diffusive wire, 
connected to two normal-metal or superconducting reservoirs 
via contacts of resistance Rm- The reservoirs are further 
connected to the macroscopic measurement apparatus (see 
Subs. III.HII . The heat flows and thermal conductances be- 
tween the studied electron system and the external driving, 
the electromagnetic environment, and the phonons in the lat- 
tice fSubs. in.Cl and Hl.DII are indicated with the arrows. The 
description of phonons in the lattice can further be divided in 
the film phonons, substrate phonons and finally the heat bath 
on which the substrate resides (Subs. HI. Gil . If the system is 
used as a radiation detector, it also couples to the radiation 
field, typically via some matching circuit (Sec. llVjl . 



main methods used in the fabrication of mesoscopic elec- 
tronic devices, and in Sec. lVIII we briefly discuss some of 
the main open questions in the field and the prospects 
of practical instruments based on electronic refrigeration 
and using the peculiar out-of-equilibrium energy distri- 
butions. 



II. THERMAL PROPERTIES OF MESOSCOPIC SCALE 
HYBRID STRUCTURES AT SUB-KELVIN 
TEMPERATURES 

The schematic picture of a setup studied in typical ex- 
periments described in this Review is shown in Fig. ^ 
The main object is a diffusive metal or heavily-doped 
semiconductor wire connected to large electrodes acting 
as reservoirs where electrons thermalize quickly to the 
surroundings. The electrons in the wire interact between 
themselves, and are coupled to the phonons in the film 
and to the radiation and the fluctuations in the electro- 
magnetic environment. The temperature Tph of the film 
phonons can, in a non-equilibrium situation, differ from 
that of the substrate phonons, Tgub and this can even dif- 
fer from the phonon temperature Tq in the sample holder 
that is externally cooled. Under the applied voltage, the 
energy distribution function f{E) of electrons depends 
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on each of these couphngs, and on the state (e.g., super- 
conducting or normal) of the various parts of the system. 
In certain cases detailed below, f{E) is a Fermi function 



exp[(E - fl)/ (fcsTejJ + 1 

characterized by an electron temperature Te and poten- 
tial ^. One of the main goals of this review is to explain 
how Te, and in some cases also Tp^, can be driven even 
much below the lattice temperature Tq, and how this low 
Te can be exploited to improve the sensitivity of appli- 
cations relying on the electronic degrees of freedom. We 
also detail some of the out-of-equilibrium effects, where 
f{E) is not of the form of Eq. (Q. In some setups, the 
specific form of f{E) can be utilized for novel physical 
phenomena. 

Throughout the Review, we concentrate on wires 
whose dimensions are large enough to fall in the quasi- 
classical diffusive limit. This means that the Fermi wave- 
length Xf, elastic mean free path £ei and the length of 
the wire L have to satisfy Xp ^ iei ^ L. In this regime, 
the electron energy distribution function is well defined, 
and its space dependence can be described by a diffusion 
equation (Eq. l|2Jl). In most parts of the Review, we as- 
sume the capacitances C of the contacts large enough, 
such that the charging energy Ec — e^/2C is less than 
any of the relevant energy scales and can thus be ignored. 

Our approach is to describe the electron energy distri- 
bution function /(r, E) at a given position r of the sam- 
ple and then relate this function to the charge and heat 
currents and their noise. In typical metal structures in 
the absence of superconductivity, phase-coherent effects 
are weak and often it is enough to rely on a semiclassical 
description. In this case, /(r, E) is described by a diffu- 
sion equation, as discussed in Subs. III.AI The electron 
reservoirs impose boundary conditions for the distribu- 
tion functions, specified in Subs. Ill.ijI The presence of 
inelastic scattering due to electron-electron interaction, 
phonons or the electromagnetic environment can be de- 
scribed by source and sink terms in the diffusion equa- 
tion, specified by the collision integrals and discussed in 
Subs. III. CI In the limit when these scattering effects are 
strong, the distribution function tends to a Fermi func- 
tion /eq(i?; Te(r), /i(r)) throughout the wire, with a po- 
sition dependent potential fj,{r) and temperature Te(r). 
In this quasiequilibrium case, detailed in Subs. III.DI it 
thus suffices to find these two quantities. Finally, with 
the knowledge of f{r;E), one can obtain the observable 
currents and their noise as described in Subs. III. El 

In many cases, it is not enough to only describe the 
electrons inside the mesoscopic wire, assuming that the 
surroundings are totally unaffected by the changes in this 
electron system. If the phonons in the film are not well 
coupled to a large phonon bath, their temperature is in- 
fluenced by the coupling to the electrons. In this case, it 
is important to describe the phonon heating or cooling in 
detail (see Subs. ^^3). Often also the electron reservoirs 



may get heated due to an applied bias voltage, which 
has to be taken into account in the boundary conditions. 
This heating is discussed in Subs. III. HI 

At the temperature range considered in this Review, 
many me tals undergo a t ransition to the superconduct- 
ing state llTinkhamLll99l . This gives rise to several new 
phenomena that can be exploited, for example, for ther- 
mometry (see Sec. IIIip . for radiation detection fSec. IIV|I 
and for electron cooling (Sec. 0). The presence of su- 
perconductivity modifies both the diffusion equation (in- 
side normal- metal wires through the proximity effect, 
see Subs. III.A|) and especially the boundary conditions 
(Subs. III.BII . Also the relations between the observ- 
able currents and the distribution functions are modified 
rSubs. HlEll . 

Once the basic equations for finding /(r, E) are out- 
lined, we detail its behavior in different types of normal- 
metal - superconductor hctcrostructurcs in Subs. III.FI 

A. Bottzmann equation in a diffusive wire 

The semiclassic al Bo ltzmann equation 

llAshcroft and MerminL Il976l : Smith and Jensenl 
119891) describes the average number of particles, 
/(r, p)(ir(ip/(27r)^, in the element {dr,dp} around the 
point {r, p} in the six-dimensional position-momentum 
space. The kinetic equation for /(r,p) is a continuity 
equation for particle flow, 

+ V • + eE • /(r, p;t) = Uf] + /;„[/]. (2) 

Here E is the electric field driving the charged particles 
and the elastic and inelastic collision integrals Ie\ and 
/in, functionals of /, act as source and sink terms. They 
illustrate the fact that scattering breaks translation sym- 
metries in space and time — the total particle numbers 
expressed through the momentum integrals of / still re- 
main conserved. 

In the metallic dif fusive limit, Eq. ^ may be 
simplified as follows llNaeaevl Il992t iRammeil Il998t 
Sukhorukov and Lossl Il999(l . The electric field term can 
be absorbed in the space derivative by the substitu- 
tion E — ep + /i(r) in the argument of the distribu- 
tion function, such that E describes both the kinetic 
Ep and the potential energy of the electron. There- 
fore, we are only left with the full r-dependent derivative 

V • V/ = V • drf + eE • dpf on the left-hand side of 
Eq. In the diffusive regime, one may concentrate on 
length scales much larger than the mean free path iei. 
There, the particles quickly lose the memory of the di- 
rection of their initial momentum, and the distribution 
functions become nearly isotropic in the direction of v. 
Therefore, we may expand the distribution function / in 
the two lowest spherical harmonics in the dependence on 

V = v/u, /(v) = /o + V • Sf, and make the relaxation- 
time approximation to the elastic collision integral with 
the elastic scattering time r, i.e., Ic\ = — v • Sf/r. In 
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the hmit where the time dependence of the distribution 
function takes place in a much slower scale than t, this 
yields the diffusion equation with a source term, 



{dt-DVl)fQ{r;E,t)^Ii,,[fo] 



(3) 



Here we assume that the particles move with the Fermi 
velocity, i.e., v — vp'v. As a result, their diffusive motion 
is characterized by the diffusion constant D = VpT/3. In 
what follows, we will mainly concentrate on the static 
limit, i.e., assume dtfo{r;E,t) = and lift the subscript 
from the angle-independent part /o of the distribution 
function. 

Equation |(2Jl can also be derived rigorously from the 
microscopic theory through the us e of the quasiclassi- 
cal Keldysh formalism ( Rammer an d Smithl . [r9 86). With 
such an approach, one can also take into acco unt super- 
condu cting effects, such as Andreev reflectio n llAndreeyl 
Il964b() and the proximity effect (Bclzig et al l Il999|) . 
In the diffusiv e limit, one obtains the Usadel equation 
llUsadelLll97nl) . which in the static case is 

• / = [-tEf, + A(r) + Sin(r), G(r; E)] . (4) 

Here G{r; E) is the isotropic part of the Keldysh Green's 
function in the Keldysh (g) Nambu space, A and a are 
the cross section and the normal-state conductivity of the 
wire, 73 = l(g)f3 is the third Pauli matrix in Nambu space, 
A = 1 (g) A is the pair potential matrix, and Ejn describes 
the inelastic scattering that is not contained in A. Us- 
adcl equat ion d escribes the matrix current / = crAGVG 
l|Nazarovill999|) . whose components integrated over the 
energy yield the physical currents. In the Keldysh space, 
G is of the form 



G = 



/G« G^ 
^ G-^ 



where each component is a 2 x 2 matrix in Nambu 
particle-hole space. Equation Q has to be augmented 
with a normalization condition G^ = 1. This im- 
plies (G^/^Y — 1 Eind allows a parametrization G^ — 
G^h — hG^ , where h is a distribution function ma- 
trix with two free parameters. The equations for the 
retarded/ advanced functions G^/^ ((Ijl) Si,nd (2,2) - 
Keldysh components of Eq. Q) describe the behavior 
of the pairing amplitude. The solutions to these equa- 
tions yield the coefficients for the kinetic equations, i.e., 
the (1,2) or the Keldysh part of Eq. This describes 
the symmetric and antisymmetric parts of the energy 
distribution function with respect to the chemical po- 
tential of the superconductors. The latter is assumed 
everywhere equal to allow a time-independent descrip- 
tion. A common choice i s a diagonal h = + f'^Ts 
(ISchmid and SchonL ^97^, where f^{E) = f{-E) ~ 
f{E) is the antisymmetric and f'^{E) = l-f{E)-f{-E) 
the symmetric part of the energy distribution func- 
tion f{E). With this choice, inside the normal metals 



where A = , we get two kinetic equations of th e form 
l|Belzig et a/.l . ll999tlvrrtanen and Heikkilall2004aj) 

DS/ ■ f = Sf„, = aA(PiV/^ - r-"v/^ + jsH , 

(5a) 

(5b) 



Here = -jTr \{t\ describes the spectral en- 

ergy current, and = jTt [(ti (g) t^)!] the spectral 
charge current. The inelastic effects are described by 
the collision integrals = ^Tr [(ti (g) i)[Si„ , G]] and 
= iTr [(ti (g)f3)[Si„, G]]. The kinetic coefficients are 



Vl = ^Tr [1-G^G^] 

Vt = ^Tt [l-G^faG^fa] 

T^" = ^Tr [G^G^fg] 

js = \tt [(G^VG« - G^VG^)f3]. 

Here, Vl and Vt are the spectral energy and charge 
diffusion cocfhcients, and js is t he spectral den sity of 
the supercurrent-carrying states ijHeikkila et all . |2002). 
The cross-term T^" is usually small but not completely 
negligible. In a normal-metal wire in the absence of a 



proximity effect, G^ = fg and G^ 



-T3. Then we 



obtain Vl = Vt — 1, T^" — is ~ ^ and the kinetic 
equations H5b|l reduce to Eq. ||2J) in the static limit. 



B. Boundary conditions 

The quasiclassical kinetic equations cannot directly de- 
scribe constrictions whose size is of the order of the Fermi 
wavelength, such as tunnel junctions or quantum point 
contacts. However, such contacts can be describe d by 
the boundary conditions derived bv iNazarovl l)l999() . 

iL=iR = Z{]^[GL.GR] + )[GR.GL]-, (6) 



where 



2p2 , 7- 



h ^ 2 + rJx- 1)' 
Here, [GL,Gfl]± = GlGi^, ± GrGl. and 4 



and 



Gl{r) are the matrix current and the Green's func- 
tion at the left (right) of the constriction, evaluated 
at the interface and flowing towards the right. The 
constriction is described by a set {7^} of transmission 
eigenvalues through the function Z[x). For large con- 
strictions, it is typically enough to transform the sum 
over the eigenvalues to an integral over the transmission 
probabilities T, weighted by their probability distribu- 
tion p{T). In the case of a tunnel barrier, 7^ ^ 1, 
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and thus Z{x) 



'Yl,n%i/h = For a bal- where 



hstic contact 7^ = 1 and Z{x) — Gi^/{x + 1). For 
other types of contacts, it is typically useful to find 
the observable for arbitrary T and weight it with p{T), 
e.g., p(T) = h- K GNl\ {2e^)TVT^] for a diffusi ve con- 
tact llNa,7a,rovl ^ij, p{ T) = HGn /\e^T ^ /^y/l^ ] 
for a dirty interface l)ScheD and Baueil Il997(l or 
p(T) = 2hGN/\e'^\/T(l -T)] for a chaotic cavity 
l|Bar anger and MelioL Il994l) . This way, the observables 



can be related to the normal-state conductance Gn of 
the junction. 

Equation © yields a boundary condition both for the 
"spectral" functions G^/^ and for the distribution func- 
tions. In the absence of superconductivity, we simply 
have G^/"^ = ifs, and the boundary condition for the 
distribution functions becomes independent of the type 
of the constriction. 



(7) 



In this case, the two currents can be included in the same 
function by defining ji±E) = {j^{E) ± f{E))/2. This 
yields the spectral current through the constriction 



j{E)^GM{fL{E)- fniE)), 



(8) 



where Jl/r is the energy distribution function on the 
left/right side of the constriction. 

Another interesting yet tractable case is the one where 
a superconducting reservoir (on the "left" of the junc- 
tion) is connected to a normal metal (on the "right") 
and the proximity effect into the latter can be ignored. 
The latter is true, for example, if we are interested in the 
distribution function at energies far exceeding the Thou- 
less energy of the normal- metal wire, or in the presence 
of strong depairing. In this case, the spectral energy and 
charge currents are 



(9a) 



(9b) 



Here E = \E\ — A and 9{E) is the Heaviside step func- 
tion, and the energy-dependent coefficients are 

2 ((2 - T)\E\n + Tn^) 
MdT) = J ' ^, „^ (10a) 



{{2^T)n + T\E\Y 
2TA^ 

4(T-l)£:2 + (r-2)2A2 
2\E\ 



{2-T)n + T\E\ 



(10b) 
(10c) 



Here we defined il = \/E'^ — A'^. In the tunneling limit 
T < 1, we get 



f/^ = GME)if'n''^-f'J^), 



(11) 



Ns{E) 



Re 



E + iT 



A2 



d{E)\E\/n 



(12) 

is the reduced superconducting density of states (DOS). 
The first form of Eq. (|12|l assumes a finite pair-breaking 
rate F, which turns out to be relevant in some cases dis- 
cussed in Sec. I^C.l Unless specified otherwise, we as- 
sume that the superconductors are of the conventional 
weak-coupling type and the superconducting energy gap 
A at T = is r elated to the cri tical temperature Tc by 
A « 1.764:kBTc llTinkhamL[l99l . 



C. Collision integrals for inelastic scattering 

The collision integral /jn in Eq. Q is due to electron- 
electron, electron-phonon interaction and the interaction 
with the photons in the electromagnetic environment. 



1. Electron-electron scattering 

For the electron-electron interaction, the collision in- 
tegral is of the form 



(13) 

where a and n^f}^ depend on the type of scattering and 
the "in" and "out" collisions are 



^coU 



[1 - f{E)][l - f{E')]f{E ~ nw)f{E' + hij) 

(14a) 

Cn = fiE)f{E')[l ~ f{E - hu^)][l - f{E' + fio,)]. 

(14b) 

Electron-electron scattering ca n be either due to 
the d irect Coulomb interaction l)Altshuler and Aronovl 
mil, or mediated through magnetic impurities 
which can flip their spin in a scattering process 
l)Kaminski and GlazmanLl200l|) or other types of impuri- 
ties with internal dynamics. In practice, both o f these ef- 
fects c ontribute to t he ene rgy relaxation l)Anthore et aO. 
l2003t iPierre et all l2000(l . Assuming the electron- 
electron interaction is local on the scale of the variations 
in the distribution function, the direct interaction yields 
l|Altshuler and Aronovl fl985l) Eq. lO with a = ~2+d/2 
for a d-dimensional wire. In a diffusive wire, the effec- 
tive dimensionality of the wire is determined by com- 
paring the dimensions to the energy-dependent length 
Le = yJhD/E. The prefactor K^f^^ for a d-dimensional 
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sample is 



fective when ivTe-e ~ 1- This yields a length scale 



.(1) 



1 



i Huard et a/.l.l2004 (15a) 



= — faammer and Smithl . Il98fil) (15b) 
SEpT 

1 



87r2\/2n2j,^£,3/2' 



I Rammer 



199 



i) (15c) 



where lyp = v{Ef) is the density of states at the Fermi 
energy Ep and A is the wire cross-section. 

In the cas e of relaxation due to magnet ic impurities, 
one expects l)Kaminski and Glazmanl . l200l|) a = — 2 and 



2 hvp ' 



S{S + 1) [in (^) 



Here 



is the 



concentration, S is the spin, and Tr- is the Kondo tem- 
perature of the magnetic impurities responsible for the 
scattering. This form is valid for T > Tk- For a more 
detailed account of the magnetic-impurity effects, see 



fCoppert et al, 2002; Goppcr t and Grabcrt, 200l ] l 2003t 
iKroha and Zawadowskil l2002t iUisaghv et all l2004j) and 

the references therein. 

For d — 3, and for small deviations Sf from 
equilibri um, the colli sion integral can be approxi- 
mated ( Rammerl Il998>") by —Sf/Te-e, where Te-e — 
SVS^iVS- l)C{?,/2){kBTf'^ /{IQkpleiVh^Ep) is the re- 
laxation time (C(3/2) « 2.612), r = leijvp is the elastic 
scattering time and hp Vo the Fermi momentum. In the 
case when a < —1/2, the usual relaxation-time approach 
does not work for the electron-electron interaction as the 
expressio n for the relaxation time ha s an infrared di - 
vergence llMtshuler and Arnnnvl [T985t lRa,mme l Il99j^ . 
Therefore, one has to solve the full Boltzmann equation 
with the proper collision integrals. To obtain an estimate 
for the length scale at which the electron-electron inter- 
action is effective, we can proceed differently. Introduc- 
ing dimensionless position x = x/L and energy variables 
E' = E' IE* and cZ) = hbjjE*, we get 

a|/ = -Ke-eicoll. 

Here the dimensionless integral /coii characterizes the de- 
viation in the shape of the distribution function from the 
Fermi function and Ke~e depends on the specific system. 
For a quasi-ld wire with bare Coulomb interaction. 




(16) 



where Rk = h/{2e'^), Rp, = L/{aA) is the resistance of 
the wire and Ep = fiD jl? is the Thouless energy. In the 
case when the wire terminates in a point contact with 
resistance Rp, the resistance i?_D in Eq. l[THl) sh ould be 
replac ed with the total resistance Ro + Rp (fekola et al\ . 
I2004a|) . Typically the energy scale characterizing the 
deviation from (quasi)equilibrium is E* — eV. At 
eV ^ fcsT, electron-electron collisions start to be ef- 



£:_,^\ RKAa 



2hD 
~eV' 



(17) 



where A is the cross section of the wire and a its conduc- 
tivity. Using a wire with resistance Rp = 10 and 
Ep » 10 ueV for L = 1 fim (close to the values in 
l)Huard et al\ . \2004} ). and a voltage V = 100 fiV, we get 
£*_g « 24 ^m. Increasing the temperature, /coii becomes 
smaller, and this eff ective length also de creases. The ex- 
perimental results of lHuard et al\ l)2n04|) indicate at least 
an order of magnitude larger Kg-e and thus smaller ^*_e 
than predicted by this theory. At present, the reasons 
for this discrepancy have not been found. 



2. Electron-phonon scattering 

Another source of inelastic scattering is due to 
phonons, for which the collision integr al is of the form 
i)Rammeil [1993 IWellstood et aZ.LIl994|) 



y-C-ph 

^coll 



27T / d^a'Fiu:) iZ, [E, uj) - [E, lo) 



(18) 



Here 



(19) 



+ f{E)]f{E^nw)n^^{Lo)) 
aUE,u;) =fiE)[l - fiE - hLj)][l + nph(c.)] 
+ f{E)[l-f{E + nij)]np^{uj). 

The kernel a'^F{ijj) (the Eliashberg function) of the in- 
teraction depends on the type of considered phonons 
(longitudinal or transverse), on the relation between the 
phonon wavelength Xph and the electron mean free path 
iei, o n the dimensionality of the electron and phonon sys- 
tem l)Sergeev et al. . l2005^■ and on the c haracteristics of 
the Fermi surface i Prunnila et a?J . l2005j) . At sub-kelvin 
temperatures and low voltages, the optical phonons can 
be neglected, and one can only concentrate on the acous- 
tic phonons. In what follows, we also neglect phonon 
quantization effects which may be important in restricted 
geometries. Moreover, the phonon distribution func- 
tion nph(w) is considered to be in (quasi)equilibrium, 
i.e., described by a Bose distribution function n-p\i(ijj) = 
ricq(w) = [eyi\){fvjj I (kpT)) — (for phonon relaxation 
processes, see Subs. Urm) . 

When the phonon temperature Tph is much lower than 
the Debye temperature Tp, the phonon dispersion rela- 
tion is linear and one can estimate the phonon wavelength 
using Xph — hvs/kBTph- For typical metals, the speed 
of sound is U5 ~ 3 . . . 5 km/s, which yields a wavelength 
Xph - 100 . . . 200 nm at Tp,, = 1 K and Ap^ - 1 . . . 2 ^m 
at Tph — 100 mK. In the clean limit Xph <C iei, approxi- 
mating the electron-phonon coupling with a scalar defor- 
mation potential, only the longitudinal phonons are cou- 
pled to the electrons. In this case for lo <^ kpTp/h, kpvs 
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Ag 



Al 



Au 



AuCu 



Cu 



Mo 



Si 



Ti 



Te (mK) 



50.. .400 



35.. .130 
200.. .300 



80.. .1200 



50.. .120 



25.. .800 
100.. .500 
50.. .200 



980 



120. ..400 
173.. .450 
320. ..410 



300.. .800 



Tph 
(mK) 



50... 400 



35 
200 



80.. .1000 



20.. .120 



25. ..320 
280.. .400 
50.. .150 



80... 980 



175.. .400 
173 

320.. .410 



500.. .800 



(W 

3K-5) 



0.5 10" 



a2i(r 

0.3 10^ 



2.4 10 



2.4 10 



2.0 10 
0.9. ..4 lO'^ 
2.0 lO'' 



0.9 10" 



0.1 10 
0.04 lO'* 
0.1 lO'' 



1.3 10 



Measured in 




tz et all 

'Mcschkc et aL, 200^ 



( Echt ernach et al.. 

19921) 

"( Wel lstood et all 
1989) 



iRoukcs et al., 1985) 
I Lcivo^^/jl_^9^^^ 
:Meschke ef adl2004l) 



Sav in and Pekolal 
2005) 



(Savin et al., 2001) 
( Pm nnila et a/., .2002ih 
( Buonomo et aZi 120031) 
^Ma nninen et all 
1999) 



TABLE I Measured electron-phonon coupling constant E for 
different materials. The second and third columns indicate 
the temperature ranges (electron and phonon temperatures, 
respectively) of the measurements. 



amme 



\M\' 



(20) 



where |Mp is the square of the matrix element for the 
deformation potential. Generally this is inversely pro- 
portional to the mass density of the ions, but its pre- 
cise microscopic form depends on the details of the lat- 
tice structure. Therefore it is useful to present \M\ in 
terms of a separately measurable quantity, e.g., the pref- 
actor E of the power P — SVT^ dissipated to the lat- 
tice of volume V in the quasiequilibrium limit (see Table 
Uand Subs. ITLD^ : |M|2 = TTh^v%i:/ {l2C{b)k%), where 
C(5) « 1.0369. 

In the dirty limit \ph S> ^el^ the power of uj in the 
Eliashberg function can be either 1 or 3, for the cases 
of static or vibrating disorder, respectively. For further 
details about the dirty limi t, we refer to feel itz, laait 
iRammer and Schm id. 1986: ISergeev and Mitinl . (200(1) ""^ 

The relaxation rate for electron-phonon scattering is 

given by l/ve-ph = -{5h-ph[I{E)]/5f{E)}\f!^E)=u^(E): 
where fcq{E) now is a Fermi function at the la ttice tem- 
perature. With this definition at E = Ep, ijRammeil 
[l998) 



47r 



-ph 



duj- 



sinh 



\kBT J 



(21) 



Thus, in the clean case for fcsT <^ 2hkpvs we ob- 
tain 1/Te-ph = aT3, a = 7C(3)E/(24C(5)fc|;/i.) » 
0.3'iT,/{kliyp). With typical values for Cu, S = 2 • 10^ 



WK 



-5, 



and i^i. = l.e-lO-*^ J 



we get Te-ph = 45 



ns at T = IK. Assuming Xp <^ t 
electron-phonon scattering length is £. 



el < 



-ph: 



the 



ph — \/ DTe—ph ■ 

For the above values and a typical diffusion constant 
D = 0.01 m^/s, te-ph w 21 ^m at T = 1 K and 
te-ph ~ 670 yum at T = 100 mK. 

In the disordered limit Xph ^ lei, the temperature de- 
pendence of the electron-phonon scattering rate is ex- 
pected to follow either the or laws, dependin g on 
the nature of the disorder ijSergeev and Mitinl l200(]|) . 

It seems that although most of the experiments are 
done in the limit where the phonon wavelength at least 
slightly exceeds the electron mean free path, in major- 
ity of the cases the results have fitted to the clean-limit 
expressions, i.e., the scattering rate oc and the heat 
current flowing into the phonon system oc Tp, see E^. (I24|) 
below (for an exception, see ijKarvonen et al\ . 12005^) ). 
Finding the correct exponent is not straightforward, as 
the film phonons are also typically affected by the mea- 
surement, and because of the reduced dimensionality of 
the phonon system (see Subs. ITlUIi . In this Review, we 
concentrate on the clean-limit expressions. 



D. Quasiequilibrium limit 

The shape of the distribution function at a given 
position of the wire strongly depends on how the 
inelastic scattering length compares to the length L 
of the wire. For L <C Im [nonequilihrium limit), we may 
neglect the inelastic scattering altogether. In this case, 
the distribution function is a solution to either Eqs. (O 
or Eq. (PJ, where the collision integrals/self energies for 
inelastic scattering can be neglected. As a result, the 
shape of the electron distribution functions inside the 
wire at a finite bias voltage eV > fepT may strongly de- 
viate from a F er mi di str ibution l[ Giazotto_e|_a|.L 12004 
Heikkila et all . 1200^ iHeshnga and Kla.pwiikL Il99c ; 



Pekola et all T2004at fpierre all . 120011: iPothier et al. 
^Q^F The nonequilihrium shape shows up in most 
of the observable properties of the system, including 
the I — V characteristics, the current noise or the 
supercurrent. In general, it can only be neglected in the 
I — V characteristics if the charge transport process is 
energy independent as in the case of purely normal- metal 
samples. Even in this case the form of f{E) can be 
observed in the current noise. 

The kinetic equations can be greatly simplified in the 
limit where kn for one type of scattering is much smaller 
than L. In the quasiequilibrium limit, the energy relax- 
ation length due to electron-elec tron scattering is much 
shorter than the wire, £e-e ^ L ijNagaevl Il995j) . In this 
case, the local distribution function is a Fermi function 
characterized by the temperature Te{r,t) and potential 
/i(r,t). Mathematically, this can be seen by considering 
the Boltzmann equation Q with the electron-electron 
collision integral, Eq. (|13|l in the limit where the pref- 
actor of the latter becomes very large. As the left-hand 
side of Eq. Q is not strongly dependent on the form of 
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f{E, r) as a function of energy, the equation can only 
be satisfied if the collision integral without the prefactor 
becomes small. It can be easily shown that the latter 
vanishes for f{E) — fcq{E). Thus, the deviations from 
the Fermi-function shape will be at most of the order of 
kn/L, and can be neglected in the quasiequilibrium limit. 

In this limit, we are still left with two un- 
knowns, Te{r,t) and fJ.{r,t). Substituting f{E,r) — 
/eq(^;re(r,t),/i(r,t)) in Eq. ® yields 

{dt - DW^jf = {dt - DV^){T,dTj + 

^[(VM)'a^/) + {VT,fdlj + 2VT,V^i^TM = Icon[f], 

where /coii [/] contains the other types of inelastic scatter- 
ings, e.g., those with the phonons. In the right hand side 
of the upper line, the differential operators dt and V act 
only on Tg and fi. Integrating this over the energy and 
multiplying by i^pE and then integrating over E yields 

{dt - DV^)fi{r) - 0, (22) 
Ce{r,t)dtTe - V(K(r,t)VTe) - cT{Vfi/ef = 7coU. (23) 

We assumed that the energy integral over /coii[/] van- 
ishes.^ Here Ce(r, i) = Tr^iypkgTe{r,t)/3 is the electron 
heat capacity, a — Dvpe^ is the Drude conductivity, 
K{v,t) = aLoTe{r,t) is the electron heat conductivity, 
Lo = 7r2/fc|/(3e2) w 2.45 • lO'^ WnK-^ is the Lorenz 
number and /coii(7e,M) contains the power per unit vol- 
ume emitted or absorbed by other excitations, such as 
phonons or the electromagnetic radiation field. The last 
term on the left hand side of Eq. H23|l describes the Joule 
heating due to the applied voltage. In what follows, 
we write the volume explicitly in the collision integral 
by averaging over a small volume V around the point 
r where T{t) is approximately constant, thus defining 

-Pcoii(r) = V/coii- 

Fo r the electron-phonon scattering l)Wellstood et all 
PcoU reads in the clean case (see also Table Pi 

P:„7'' = SV(T,(r)5-Tp,(r)5). (24) 

For the dirty limit specified below Eq. (|20|l . the Eliash- 
berg functions scaling with translate into temper- 
ature dependences scalin g as T"+'^, i.e., and 
l)Sergeev and Mitinl . l2000(l . 

The electrons can also be heated due to the thermal 
noise in their electromagnetic environment unless proper 
filtering is realized to prevent this heating. If one aims 
to detect the electromagnetic environment as discussed in 
Sec. II VI this discussion can of course be turned around to 
find the optimal coupling to the radiation to be observed. 
A model for such coupling in the quasiequilibrium limit 



In the diffusive limit wliere tlie inelastic scattering rates are lower 
than l/r, this is related to the particle number conservation and 
is thus generally valid. 



was considered bv lSchmidt et g/l l|2004a|) . who obtained 
an expression for the emitted/absorbed power due to the 
external noise in the form 

P:or = r^iT!~T^). (25) 

Here r = AR^R^ / {Re + -^7)^ is the coupling constant, 
T-y is the (noise) temperature of the environment, and 
Re and i?-y are the resistances characterizing the ther- 
mal noise in the electron system and the environment, 
respectively. This expression assumes a frequency in- 
dependent environment in the relevant frequency range. 
Fo r some examples on the frequency dependence, we refer 
to lSchmidt et al\ l)2004a(l . 



E. Observables 

1. Currents 

In the noncquilibrium diffusive limit, the charge cur- 
rent in a normal-metal wire in the absence of a proximity 
effect or any such interference effects as weak localization 
is obtained from the local distribution function by 

/oo 
dED{E)v{E)\/ f{x]E) (26) 
-00 

and the heat current from a reservoir with potential 11 is 



Q = -A dE{E ^ ^i)D{E)v{E)Vf{x\E). (27) 

Here we included the possible energy dependence of the 
diffusion constant D{E) and of the density of states i^{E), 
due to the energy dependence of the elastic scattering 
time, or due to the nonlinearities in th e quasiparticle 
dispersion relation. If the Kondo effect ijVavilov et all 
2003) can be neglected, such effects are very small in 
good metals at temperatures of the order of 1 K or less. 

In the quasiequilibrium limit, assuming D{E) =const. 
and v{E) = vp, Eqs. (|26I27|I can be simplified to 

/ = -ctAV^/c (28a) 
Q = -kAVT. (28b) 

When a diffusive wire of resistance Rp) is connected 
to a reservoir through a point contact characterized by 
the transmission eigenvalues {^}, the final distribution 
function is obtained after solving the Boltzmann equation 
^ or Eqs. with the boundary conditions given by 
Eq. ®, ©I or by ®. However, when Rp, is much less 
than the normal-state resistance 1/Gn of the contact, 
we can ignore the wire and obtain the full current by a 
direct integration over Eq. (jS)) or Eq. For example, 
for NIS or SIS tunnel junctions, the expressions for the 
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charge and heat currents from the left side of the junction 
become 



1 



Q 



1 

c'^Rt 



dENL{E)NR{E)[fL{E) ^ fniE)] (29) 
dEENLiE)NRiE)[fLiE)- fniE)]. (30) 



Here E = E - eV, Rt = 1/Gn, Nl/r{E) = v{E)/vf 
is the reduced density of states for the left/right wire, 
N{E) = 1 for a normal metal and N{E) = Ns{E) for a 
superconductor. Furthermore, if the two wires constitute 
reservoirs, /l/r are Fermi functions with potentials fJ.L = 
—eV, fiR = 0. The resulting NIS or SIS charge current is 
a sensitive probe of temperature and can hence be used 
for thermometry or radiation detection, as explained in 
Sees. IIII.AI and IIVI respectively. Moreover, analysis of 
the heat current Eq. H^Ofl shows that the electrons can 
in certain situations be cooled in NIS/SIS structures, as 
discussed in Sec. IV.C.ll 

In the presence of a proximity effect, the equations 
for the charge and thermal currents in the quasiclassical 
limit, i.e., ignoring the energy dependence of the diffusion 
constant and the normal-metal density of states, are 



J dEf 
f dE{Ef - 



(31a) 
(31b) 



Here fj, is the potential of the reservoir from which the 
heat current is calculated. As in Eq. Q, these currents 
can be separated into quasiparticle, anomalous, and su- 
percurrent parts. 



2. Noise 

Often one can express the zero-frequency current 
noise in terms of the l ocal distribution function 
l|Blanter and Biittikeilf200(]() . The noise is characterized 
by the correlator 



5 = 2 



dt'{Si{t + t')Si{t)), 



where 61 = I — (/) and / is the current operator. In a 
stationary system S is independent of t. In a normal- 
metal wire of length L in the noneq uilibrium limit , the 
current noise S can be expressed as (|Nagaevl 1199^ 



S = 



L2 



dx / dEfiE,x)[l~ fiE,x)]. (32) 



In the quasiequili brium regime, this equation simplifies 
to ijNagaevl [T995|) 



5 = 



l2 



dxTe{x). 



(33) 



When the resistance of a point contact dominates that 
of the wire, the noise power can be expressed through 
l|Blanter and Biittikeil . |2000) 

Snn I [/^(l - + /^(l - /fl)] 

n 

+ %,{l-%,)(h-fRf} 



for a normal-metal contact and l)de Jong and Beenakkeil 
EM) 



2e2 ^ f 



r2 



2fL{E)[i~fLm 
ifir 



(2 - r„)2 

, i6r„2(i-rj^^^^2 



(2 - T„) 

for an incoherent NS contact at £' < A. Here /^/^ 
are the distribution functions in the left/right (nor- 
mal/superconducting in the latter case) side of the con- 
tact and fliE) = 1 — //.(i?) — fL^—E), the symmetric 
part w.r.t. the S potential. If the scattering probabili- 
ties are independent of energy and if the two sides are 
reservoirs with the same temperature, these expressions 
simplify to 

Snn^'^Y. [^^bTT^ + eV coth {v) T„(l - r„)] 



(2 - r„)2 



2fcsT 



i6T„2(i-r„) 

(2 - TaY 



(2eV"coth {v)-AkBT) 



(34a) 
(34b) 



where v = eV/{2kBT). For eV = (thermal Johnson- 
Nyquist noise), S^n/s = '^kBTG^/S} where GnRk = 
Y^^Tn and GsRk = E„'4V(2 - T^f are the conduc- 
tances of the point contact in the NN and NS cases, re- 
spectively. In the opposite limit eV ^ fc^T (shot noise), 
one obtains S = ieF^/slN/s^ where In/s = G^^/sV 
is the average current through the junction, and i^jv/s 
are the Fano factors: F^GnRk = Yn'^ni^ — %i) and 
i^sGsi?if = E„ 16^'(1 - ^n)/(2 - T„)f ^ 

In the presence of the superconducting proximity 
effect, the e xpression for the n o ise b ecomes more 
complicated l)Houzet and Pistolesil l2004|) . In gen- 
eral, it can be found by employing the counting-field 
technique developed by N azarov and coworkers, see 
l)Nazarov and Bagretsl l2002(l and the references therein. 
This technique can also be applied to study the full 
counting statistics of the transmitted currents through a 
given sample within a given measurement time ijNazarovl 
120031) . 

Apart from the charge current, also the heat current 
in electric circuits fluctuates. For example, the zero- 
frequency heat current noise from the "left" of a tunnel 
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FIG. 2 Nonequilibrium quasiparticle energy distribution 
function in a diffusive normal-metal wire in the absence of 
inelastic scattering: (a) wire placed between two normal- 
metal reservoirs and (b) wire placed between a normal-metal 
{x = L) and a superconducting (a; = 0) reservoir. In the 
latter picture, we assume hD/L^ <^ eV <^ A, such that 
the proximity effect and the states above the gap can be ne- 
glected. In both pictures, the lattice temperature was fixed 
to ksT = eV/W. 



5; 0.6 
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FIG. 3 (Color in online edition): Quasiparticle energy distri- 
bution function in the center of a normal-metal wire placed 
between two normal-metal reservoirs and biased with a volt- 
age eV = lOfeflT. The three lines correspond to three ex- 
treme limits: L <^ £e~e, le-ph (nonequilibrium limit, black 
solid line), ie-e <^ L <^ f-e-ph (quasiequilibrium limit, red 
dashed line), and £e-e,£e-ph ^ L (equilibrium limit, blue 
dotted line). 



contact biased with voltage V is given by 

At low voltages V ^ ksT/e, the heat current noise obeys 
the fluctuation-dissipation result Sq = AkBT^Gth, where 
Gth is the thermal conductance. This quantity is related 
to the noise equivalent power (NEP) discussed in the lit- 
erature of thermal detectors by Sq =NEP^. The total 
NEP contains contributions not only from the electrical 
heat current noise, but also from other sources, such as 
the direct charge current noise and electron-phonon cou- 
pling. A detailed discussion of various NEP sources is 
presented in Sec. IIVI 

Another important quantity is the cross-correlator be- 
tween the current and heat current fluctuations. At zero 
frequency, this is given by 

Si, — ^JdEENUE~eV)N,iE) ^^^^ 

[/fl(l-/L) + /L(l -/«)]• 

These types of fluctuations have to be taken into ac- 
count for example when analy zing the NEP of bolometers 
(|Golubev and Kiizminl . l200l t) . Recently, also the general 
statistics of the heat c urrent fluctuations have be e n theo - 
retically addressed by iKindermann and Pilgranil l)2004j) . 



F. Examples on different systems 

Below, we detail the solutions to the kinetic equations, 
© or (jSJ in some example systems. The aim is first to 
provide an understanding of the general behavior of the 
distribution function in these systems, but also to show 
the generic features, such as the electron cooling in NIS 
junctions. 



1. Normal-metal wire between normal-metal reservoirs 

The simplest example is a quasi-one-dimensional 
normal- metal diffusive wire of resistance Rd = L/{Aa), 
connected to two normal-metal reservoirs by clean con- 
tacts. In the full nonequilibrium limit, we find (see Fig.|5] 
(a) ; the coordinate x follows Fig. ^| 

f{E,x) = h{E) + IfniE) - hiE)]^, (37) 

where and fji are the (Fermi) distribution functions 
in the left and right reservoirs with temperatures and 
Tn and potentials hl and = hl + eV , respectively. 
The resulting tw o-step form illustrate d in Fig. Et- was 
first measured bv Pothi er et ol\ l|l997b() . 

In the quasiequilibrium limit, we get f(E,x) = 
fcq{E;^i{x),Te{x)), where 

X 

fi{x) — fii^ + eV— (38a) 

T,(.)^ = r|-H(r|-r|)|-f£|(i-|). (38b) 

In both cases, the current is simply given by / = V/Rd 
and the heat current in the limit V = Qhy Q = Lq{T-^ — 
Tf^)/Ru. For y ^ 0, the resistor dissipates power and 
Q is not conserved. The electron distribution function in 
the center of the wire, x ~ L/2 is plotted in Fig.Ofor the 
nonequilibrium, quasiequilibrium and local equilibrium 
(strong electron-phonon scattering) Hmits. 

To obtain estimates for the thermoelectric effects due 
to the particle-hole symmetry breaking, let us lift the as- 
sumption of energy independent diffusion constant and 
density of states, expanding them as D{E) « Dq -I- 
cd ^e^'' and i^iE) « vp + cjM ^^^^ ■ For hnear response, 
the resulti ng express ions for the charge and heat currents 
are ((Cutler and Mott. 1969) 

/ = -GLVfi/e + GLaVT (39a) 

Q = -nCLVfi/e + GthLVT. (39b) 
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FIG. 4 Electron cooling at NIS junctions, (a) Sketch of the 
energy band diagram of a voltage biased NIS junction. Upon 
biasing the structure, the most energetic electrons (e) can 
most easily tunnel into the superconductor. As a result the 
electron gas in the N electrode is cooled, (b) Maximum cooling 
power surface density Qa vs interface transmissivity T at 
different temperatures calculated for a NS contact. 

Here G = c^i'fDqA/ L is the Drude conductance, 
Gth = LqGT = kA/L {Wiedemann- Franz law) is the 
heat conductance, a — eLoG'T/G {Mott law) is the 
Seebeck coefficient describing the thermoelectric power, 
n = aT {Onsager relation) is the Peltier coefBcient, and 
G" = e'^{cDi^F + DqCm)A/{LEf) describes effects due 
to the particle-hole symmetry breaking. We see that the 
thermoelectric effects are in general of the order kBT/Ep; 
in good metals at temperatures of the order of 1 K they 
can hence be typically ignored. Therefore, the Peltier re- 
frigerators discussed in Subs. rvTBl relv on materials with 
a low Ep. 



2. Superconducting tunnel structures 

Consider a NIS tunnel structure, coupling a large su- 
perconducting reservoir with temperature T^^s to a large 
normal-metal reservoir with temperature T^^n via a tun- 
nel junction with resistance Rt- Let us then assume 
a voltage V applied over the system. In this case, the 
heat current (cooling power) from the normal metal is 
given by Eq. ^ with Nl{E) = 1, Nr{E) = NsiE), 
fL{E) = U^[E - eV,T,^N)_ and fR{E) = f,^{E,T,^s)- 
For small pair breaking inside the superconductor, i.e., 
r <C A, Qnis is positive for eV < A, i.e., it cools 
the normal metal. It is straightforward to show that 
QmsiV) — Qnis(— ^^)- This is in contrast with Peltier 
cooling, where the sign of the current determines the di- 
rection of the heat current. For eV > A, the current 
through the junction increases strongly, resulting in Joule 
heating and making Qnis negative. The cooling power is 
maximal near eV ~ A. 

In order to understand the basic mechanism for cool- 
ing in such systems, let us consider the simplified energy 
band diagram of a NIS tunnel junction biased at voltage 
V, as depicted in Fig. Efa). The physical mechanism 
underlying quasiparticle cooling is rather simple: owing 
to the presence of the superconductor, in the tunneling 




FIG. 5 NS point contact characteristics, (a) Heat current 
from the N side as a function of voltage for different trans- 
parencies T (from top to bottom, tunneling limit T —> 0, 
T = 0.005, T = 0.01, T = 0.05, T = 0.5 and T = 1) at 
ksT = 0.3A. For T < 0.05, the heat current is positive, cor- 
responding to cooling, (b) Current-voltage characteristics for 
the same values of T and T as in (a) (now T increases from 
bottom to top). The first four curves lie essentially on top 
of each other. The corresponding NIN curves are shown with 
the dashed lines. 



process quasiparticles with energy E < A cannot tunnel 
inside the forbidden energy gap, but the more energetic 
electrons (with E > A) are removed from the N elec- 
trode. As a consequence of this "selective" tunneling of 
hot particles, the electron distribution function in the N 
region becomes sharper: the NIS junction thus behaves 
as an electron cooler. 

The role of barrier transmissivity in governing 
heat flux across the NI S structure was analyzed by 
iBardas and Averinl l)l995j) . They pointed out the inter- 
play between single-particle tunneling and Andreev re- 
flection ifAndreev... J|j964a) on the heat current. In the 
following it is useful to summarize their main results. 

The cooling regime requires a tunnel contact. The ef- 
fect of transmissivity is illustrated in Fig. 01b), which 
shows the maximum of the heat current density (i.e., the 
specific cooling power) Qa versus interface transmissiv- 
ity T at different temperatures. This can be calculated 
for a generic NS junction using Eqs. ©, HlUa|) and Ij31b|l . 
The quantity Qa is a non-monotonic function of interface 
transmissivity, vanishing both at low and high values of 
T. In the low transparency regime, Qa turns out to 
be linear in T, showing that electron transport is domi- 
nated by single particle tunneling. Upon increasing bar- 
rier transmissivity, Andreev reflection begins to dominate 
quasiparticle transport, thus suppressing heat current ex- 
traction from the N portion of the structure. The heat 
current Qa is maximized between these two regimes at 
an optimal barrier transmissivity, which is temperature 
dependent. Furthermore, by decreasing the latter leads 
to a reduction of both the optimal T and of the trans- 
missivity window where the cooling takes place. In real 
NIS contacts used for cooling applicati ons, the ave r age T 
is ty pically in the range 10~^...10~* ijLeivo et al\ . 1 19961 
Nahu m et all ^94) corresponding to junction specific 
resistances Rc (i.e., the product of the junction normal 
state resistance and the contact area) from tens to sev- 
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FIG. 6 Nonequilibrium distribution inside the SINIS island, 
calculated from Eq. (42): (a) No inelastic scattering, Xcoii = 
0, T = O.lTc, and using a depairing strength F — 10~'*A 
inside the superconductors. The voltage runs from zero to 
V = 3A/e in steps of A/(2e) and spans three regimes: (i) 
anomalous heating regime {eV = ... A, solid lines , widen - 
ing with an increasing V) discussed in (Pc kola et a/.L l2004a^ . 
(ii) cooling regime {eV — 1.5 A and eV — 2 A, dashed lines, 
narrowing with an increasing V) and the strong nonequi- 
librium heating regime {eV — 2.5A and eV — 3A, dotted 
lines, widening with an increasing V). (b) Distribution at 
eV — 2.5A for different strengths of the electron-electron in- 
teraction, parametrized by the parameter K^-e from Eq. 11611 
with Rd replaced by Rt and E* by A, as in this case these 
describe the collision integral. 



eral thousands O/Ltrn^. This limits the achievable Qa 
to some pW//zm^. From the above discussion it appears 
that exploiting low- i?c tunnel contacts is an important re- 
quirement in order to achieve large cooling power through 
NIS junctions. However, in real low-i?c barriers, pinholes 
with a large T appear. They contribute with a large Joule 
heating (see Fig.|31 which shows the heat and charge cur- 
rents through the NIS interface as functions of voltage for 
different T), and therefore tend to degrade the cooling 
performance at the lowest temperatures, and to overheat 
the superconductor at the junction region due to strong 
power injection. Barrier optimization in terms of both 
materials and technology seems to be still nowadays a 
challenging task (see also Sec. IVI.F.ip . 

In a SlNlS system in the quasiequilibrium limit, the 
temperatures T^^n and T^^s of the N and S parts may 
differ. If the superconductors are good reservoirs, Te.s 
equals the lattice temperature. In the experimentally 
relevant case when the resistance Ro of the normal-metal 
island is much lower than the resistances Rl and Rji of 
the tunnel junctions, the normal-metal temperature T(,^n 
is determined from the heat balance equation 



QSINIS(^; Te,SiTe 



N 



coll, 



(40) 



where Qsinis = "^Qnis and Pcoii describes inelastic scat- 
tering due to phonons (Eq. (|24|) ') and/or due to the elec- 
tromagnetic environment (Eq. (|25|) '). 

For further details about NIS/SINIS coohng in the 

Jiuasiequilibrium limit^ we refer to Sec. IV.C.ll and 
Anghcl and Pckola, 2Q0l|). 

Let us now consider the limit of full nonequilibrium, 
neglecting the proximity effect from the superconductors 



on the normal-metal island.^ Then the distribution func- 
tion inside the normal metal may be obtained by solving 
the kinetic equation Eq. J^J) in the static case along with 
the boundary conditions given by Eq. JSJ . For simplicity, 
let us assume Rl ^ Rr = Rt ^ Rd- In this limit, the 
distribution function /n^E) in the normal-metal island 
is almost independent of the spatial coordinate. Then, in 
Eq. we can use d^f w {dxf{L) - d.j.J{0))/L, where L 
is the length of the N wire. From Eq. Hll() we then get 

^ (NuifN - In) - Nl ih - /jv)) = tdXcoII, (41) 

where td = Lj^^/D is the diffu sion time through 
the island. As a result (,Giazotto et all l2004b: 
iHeslinea and Klapwiikl Il993|) . the distribution function 
in the central wire can be expressed as 



N 



fnNn + /lNl + RiVe'^iyFlcoU 

Nr + Nl 



(42) 



Here V = LA is the volume of the island. In the 
presence of inelastic scattering, this is an implicit equa- 
tion, as Icoii is a functional of /jv- Examples of the 
effect of inelastic scattering have been considered in 
l|Heslinga and Klap wiik|.|l993|) in the r elaxation time ap- 
proach. and in Iciazotto al\ . l2004bl) including the full 
electron-electron scattering collision integral. The dis- 
tribution function Jn in Eq. (|42|) is plotted in Fig. El for 
a few values of the voltage for Xcoii = and for a few 
strengths of electron-electron scattering at eV = 2.5A. 



3. Superconductor-normal-metal structures with transparent 
contacts 

If a superconductor is placed in a good electric con- 
tact with a normal metal, a finite pairing amplitude 
is induced in the normal metal within a thermal co- 
herence length = \/ hD/ {271 ksT) near the inter- 
face. This superconducting proxi mity effect modifies 
the properties of the normal metal l)Belzig et all Il999t 
iLam bert and RaimondiL Il998(l and makes it possible to 
transport supercurrent through it. It also changes the 
local distribution functions by modifying the kinetic co- 
efficients in Eq. jSJ- Superconducting proximity effect 
is generated by Andreev reflection (Andrcev, 1964bJ at 
the normal-superconducting interface, where an electron 
scatters from the interface as a hole and vice versa. An- 
dreev reflection forbids the sub-gap energy transport into 
the superconductor, and thus modifies the boundary con- 
ditions for the distribution functions. In certain cases, 
the proximity effect is not relevant for the physical ob- 
servables whereas the Andreev reflection can still con- 
tribute. This incoherent regime is realized if one is in- 
terested in length scales much longer than the extent of 



This is justified in the limit Rj^, Rn S> Rd- 
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the proximity effect. Below, we first study such "pure" 
Andreev-reflection effects, and then go on to explain how 
they are modified in the presence of the proximity effect. 

In the incoherent regime, for i? < A, Andreev reflec- 
tion can be describ ed through the boundary conditions 
l)Pierre et all l200l|) 

f{^Jis + E)^l- f{ps - E) (43a) 
n- V[/(/.s-S)-/(ms + ^)] -0, (43b) 

evaluated at the normal-superconducting interface. Here 
n is the unit vector normal to the interface and /is is the 
chemical potential of the superconductor. The former 
equation guarantees the absence of charge imbalance in 
the superconductor, while the latter describes the van- 
ishing energy current into it. For E > IS., the usual 
normal-metal boundary conditions are used. Note that 
these boundary conditions are valid as long as the resis- 
tance of the wire by far exceeds that of the interface; in 
the general case, one should apply Eqs. ®. 

Assume now a system where a normal-metal wire is 
placed between normal- metal (at x — L, potential /ijv 
and temperature T — T^^n) and superconducting reser- 
voirs (at a; = 0, /15 = 0, T = T^^s)- Solving the Boltz- 
mann equation then yields IjNagaev and Buttiker. 
120011) 

f(E) = / Hi + f ) fNiE) + i (1 - f ) fNiE), E<A 
^ ' \(l-f)/eq(£;;A* = 0,T,.s) + f/Ar(i?), E>A 

(44) 

Here fN{E) = fcq{E; fJ,N,Te,N) and fN{E) = 
feq{E;—^N,Te^N)- This function is plotted in Fig. [21 
(b). In the quasiequilibrium regime, the problem can 
be analytically solved in the case eV^ksT <C A, i.e., 
when Eqs. H43I) apply for all relevant energies. Then the 
boundary conditions are /i = /X5 = and n • VTg — 
at the NS interface. Thus, the quasiequilibrium distri- 
bution function is given by f{E,x) — f{E;^{x),Te{x)) 
with fJ.{x) = /iiv(l — x/E) and 



Te{x) = 




Note that this result is independent of the temperature 
Tg^S of the superconducting terminal. 

In this incoherent regime, the electrical conductance is 
unmodified compared to its value when the superconduc- 
tor is replaced by a normal-metal electrode: Andreev re- 
flection effectively doubles both the length of the normal 
conducto r and the conduct ance for a single transmission 
channel ijBeenakkeii and thus the total conduc- 

tance is unmodified. Wiedemann-Franz law is violated 
by the Andreev reflection: there is no heat current into 
the superconductor at subgap energies. However, the 
sub-gap current induces Joule heating into the normal 
metal and this by far overcompensates any cooling effect 
from the states above A (see Fig.|5}. 



The SNS syst em in the incoherent regi me has also 
been a nalyzed bv lBezudvi et all l|200(]|) and lPierre et al\ 
l|200ll) . Using the boundary conditions in Eqs. at 
both NS interfaces with different potentials of the two su- 
perconductors leads to a set of recursion equations that 
determine the distribution functions for each energy. The 
recursion is terminated for energies above the gap, where 
the distribution functions are connected simply to those 
of the superconductors. This process is called the mul- 
tiple Andreev reflection: in a single coherent process, a 
quasiparticle with energy E < —A entering the normal- 
metal region from the left superconductor undergoes mul- 
tiple Andreev reflections, and its energy is increased by 
the applied voltage during its traversal between the su- 
perconductors. Finally, when it has Andreev reflected 
~ (2 A/el/) times, its energy is increased enough to over- 
come the energy gap in the second superconductor. The 
resulting energy distribution function is a staircase pat- 
tern, and it is described in detail in (Pierre et ai, 2001). 
The width of this distribution is approximately 2A, and 
it thus corresponds to extremely strong heating even at 
low applied voltages. 

The superconducting proximity effect gives rise to two 
types of important contributions to the electrical and 
thermal properties of the metals in contact to the super- 
conductors: it modifies the charge and energy diffusion 
constants in Eqs. (|3J), and allows for finite supercurrent 
to flow in the normal- metal wires. 

The simplest modification due to the proximity effect is 
a correction to the conductance in NN'S systems, where 
N' is a phase-coherent wire of length L, connected to a 
normal and a superconducting reservoir via transparent 
contacts. In this case, js = — and the kinetic 
equation for the charge current reduces to the conser- 
vation of j"^ = VtOxJ^ ■ This can be straightforwardly 
integrated, yielding the current 

/ = GN/e J dEf^^iE; V)Dt{E), (46) 

where f^q{E;V) = {tfinh[{E + eV)/{2kBT)]-tiiiih[{E ~ 
eV) / {2kBT)]} /2, T is the temperature in the normal- 
metal reservoir, and we assumed the normal metal in po- 
tential —eV. The proximity effect can be seen in the 
term Dt = iJo ^5^) / L)-' ■ For T = 0, the dif- 
ferential conductance is dl/dV = DT{eV)GN ■ A de- 
tailed investigation of Dt{E) requires typically a numer- 
ical solution o f the retarded/advaii ced part of the Us- 
adel equation ijGolubov et alV [l997j) . In general, it de- 
pends on two energy scales, the Thouless energy Et = 
hD/L^ of the N' wire and the superconducting energy 
gap A. The behavior of the differential conductance 
as a function of the voltage and of the linear conduc- 
tance as a function of the temperature a re qualitatively 
similar, exhibiting the reentrance effect jCharla t et al\ . 
119961: iGolubov 6^0,/! Il997t Iden Hartog et oiJ . ll997D : for 
eV,kBT < Et and for max(ey, fc^T) » Et, they tend 
to the normal-state value whereas for intermediate 
voltages/temperatures, the conductance is larger than 
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Gn, showing a maximum for max(eV, ksT) of the order 
of Et- 

The proximity-effect modification to the conductance 
can be tuned in an Andreev interferometer, where a 
normal-metal wire is connected to tw o normal-metal 
reser v oirs and two sup e rconductors ij ^^l^^v e^^J. 
119971: iNazarov and Stooi Il996t iPothier et all Il994() . 
This system is schematized in the inset of Fig. [7| Due 
to the proximity effect, the conductance of the normal- 
metal wire is approximatively of the form G{(f>) — Gpf + 
6G{1 + cos(0))/2, where ip is the phase difference be- 
tween the two superconducting contacts and 6G is a posi- 
tive temperature and voltage-dependent correction to the 
conductance. Its magnitude for typical geometries is at 
maximum some O.lGjv- The proximity-induced conduc- 
tance correc tion is widely studied in the literature and 
we re fer to l)Belzig et all Il999t iLambert and Raimondl 
Il998|) for a more detailed list of references on this topic. 

The thermal conductance of Andreev interferometers 
has been studied very recently. The formulation of the 
problem is very similar as for the conductance correction. 
For i? <C A, there is no energy current into the supercon- 
ductors. Therefore, it is enough to solve for the energy 
current = V^Oxf^ in the wires 1, 2 and 5. This yields 

Q = Gn/c^ J dEEDLiE)if^\iE; T^) ~ f^^{E- T,)), 

where ft^{E;T) = ia.-ah[E / {2kBT)] and the thermal 
conductance correction is obtained from Dl{E) = 
(/p x>l'{x-e) /^^~^- Here the spatial integral runs 
along the wire between the two normal-metal reser- 
voirs. In the general case, has to be calcu- 
lated numerically. The thermal conductance corre ction 
has been analyzed byfe czuglvi and Vinokur l|2003() and 
Ijiang and Chandrasekhad (|2005bj). They found that it 
can be strongly modulated with the phase 0: in the 
short-junction limit where Et S> A, for = 0, Dl al- 
most vanishes, whereas for (j) = n, Dl approaches unity 
and thus the thermal conductance approaches its normal- 
state value. For a long junction with Et < A, the effect 
becomes smaller, but still clearly observable. The first 
measurements (Jiang and Chan drase khar, 2005a) of the 
proximity-induced correction to the thermal conductance 
show the predicted tendency of the phase-dependent de- 
crease of K compared to the normal-state (Wiedemann- 
Franz) value. 

Prior to the experiments on heat conductance 
in proximity structures, the thermoelectric power 
a was ex perimentally studied in Andreev in terfer- 
ometers lioikin g/i |2002albl: lEom a/1 Il998 : 
iJiang and Chandrasekha 

3 120041: I Parsons et all 
l2003allJ) . The observed thermopower was surpris- 
ingly large, of the order of 100 neV/K — at least one to 
two orders of magnitude larger than the thermopower 
in normal-metal samples. Also contrary to the Mott 
relation (c.f., below Eqs. this value depends 

nonmonotonically on the temperature ( Eom et all 11993) 
at the temperatures of the order of a few hundred mK, 
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FIG. 7 Supercurrent-induced N-S thermopower ajvs for the 
system depicted in the inset. Solid line: qjvs at Ti « T2. Dot- 
ted line: approximation 1471 . The correction 14811 accounts 
for most of the difference. Dashed line: ajvs at fcsTi = 3.6Et 
with varying T2 . Dash-dotted line: the corresponding approx- 
imation. Inset: Andreev interferometer where thermoelec- 
tric effects are studied. Two superconducting reservoirs with 
phase difference (j) are connected to two normal-metal reser- 
voirs through diffusive norma l-metal wires. Adapted from 
(IVhtanen and Heikkill 12004b') . 

and ev en a sign change could be found ijParsons et all 
l2003bj) . Moreover, the thermopower was found to oscil- 
late as a function of the phase (j). The symmetry of the 
oscillations has in most cases been found to be antisym- 
metric with respect to (f> = 0, i.e., a vanishes for (j) — 0. 
However , in some measure m ents, the Chandrasekhar's 
group ijEom et all Il998t iJiang and Chandrasekhail 
I2OO4I) have found symmetric oscillations, i.e., a had the 
same phase as the conductance. 

The observed behavior of the thermopower is not com- 
pletely understood, but the major features can be ex- 
plained. The first theor e tical p redictions were given by 
IClaughton and Lamberti l)l996() . who constructed a scat- 
tering theory to describe the effect of Andreev reflection 
on the thermoelectric properties of proximity systems. 

Based on their work, it was shown (Heikkila et al 200d|) 

that the presence of Andreev reflection can lead to a vio- 
lation of the Mott relation. This means that a finite ther- 
mopower ca n ari se even in the prese nce o f electron-hole 
symmetry. ISeviour and Volkov (2000bV iKogan et al\ 
l|2002|) . and lVirtanen and Heikkila. (j2004a<bi) showed that 
in Andreev interferometers carrying a supercurrent, a 
large voltage can be induced by the temperature gra- 
dient both between the normal-metal reservoirs and be- 
tween thfi_jioraml_jnetaI^^ con- 
tacts. I Virtanen and Heikkilal(l2004bD showed that in long 
junctions, the induced voltages between the two normal- 
metal reservoirs and the superconductors can be related 
to the temperature dependent equilibrium supercurrent 
/s(T) flowing between the two superconductors via 

_ 1 ^5(2i?4/3 + R^)RVA [IsjTl) ~ /s(T2)] ^^^^ 

Here i?NNN = -Ri + i?2 + ^5, Rsns — R3 + R4 + R5 
and Rk are the resistances of the five wires defined 
in the inset of Fig. [71 This is in most situations 
the dominant term and it can be also phenomenolog- 
ically argued based on the temperature dependence of 
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the supercurrent and the conservation of total current 
(supercurrent plus quasipartic le current) in the circuit 
/irtanen anc A similar result can also 

be obtained in the quasiequilibrium limit for the linear- 
response thermopower l|Virtanen and Hcikkila . u2Q04a) ■ 
In addition to this term, the main correction in the long- 
junction limit comes from the anomalous coefficient T*^", 



eVl 



/2 



Here {T^' 



i?NNN 



an \ 
/2) 



RnnnRsns 



(48) 



L^Jo -JO dErr{E)[ft^{E,T,) - 
f^^{E,T2)] and Lj. is the length of wire k. One finds 
that for a "cross" system without the central wire (i.e., 
i?5 = 0), this term dominates V^^^. Further corrections 
to t he result lITTIl are discussed b y iKogan et al\ l|2002l) 
and lVirtanen and Heikkilal l)2004a|) . 

The above theoretical results explain the observed 
magnitude and temperature dependence of the ther- 
mopower and also predict an induced voltage oscillating 
wi th the phase (j). However, th e thermopower calcu lated 
bvlSeviour and Volkovl |2000lJl. iKoean et al\ ll2002() and 
IVirtaneiranT^eikkila ( 2004alb ) is always an antisym- 
metric function of 0, and it vanishes for a vanishing su- 
percurrent in the junction (including all the correction 
terms). Therefore, the symmetric oscillations of a can- 
not be explained with this theory. 

The presence of the supercurrent breaks the time- 
reversal symmetry and hence the Onsager relation H = 
Ta (see Eqs. and below) need not to be valid for 
(/) ^ 0. Heikkila, et al. predic t ed a nonequilibrium 
Peltier-type effect l|Heikkila et all 1200,^) where the su- 
percurrent controls the local effective temperature in an 
out-of-equilibrium normal-metal wire. However, it seems 
that the induced changes are always smaller than those 
due to Joule heating and thus no real cooling can be re- 
alized with this setup. 

Recently, the thermoelectric effects in coherent 
SNS Josephson point contacts have been analyzed in 
ijZhao et ai!.Ll2nO,'il2nol . In such point contacts, the heat 
transport is strongly influenced by the Andreev bound 
states forming between the two superconductors. 



G. Heat transport by phonons 

When the electrons are thermalized by the phonons, 
they may also heat or cool the phonon system in the film. 
Therefore, it is important to know how these phonons 
further thermalize with the substrate, and ultimately 
with the heat bath on the sample holder that is typically 
cooled via external means (typically by either a dilution 
or a magnetic refrigerator). Albeit slow electron-phonon 
relaxation often poses the dominating thermal resistance 
in mesoscopic structures at low temperatures, the poor 
phonon thermal conduction itself can also prevent full 
thermal equilibration throughout the whole lattice. This 



is particularly the case when insulating geometric con- 
strictions and thin films separate the electronic structure 
from the bulky phonon reservoir (see Fig. QJ. In the 
present section we concentrate on the thermal transport 
in the part of the chain of Fig.^bcyond the sub-systems 
determined by electronic properties of the structure. 

The bulky three-dimensional bodies cease to conduct 
heat at low temperatures according to the well appreci- 
ated K (X law in crystalline solids arising from Debye 
heat capacity via 



K = Cus4i,ph/3, 



(49) 



where k is the thermal conductivity and C is the heat 
capacity per unit volume, vs is the speed of sound and 
fpi.nh is the mean free path of phonons in the solid 
l|Ashcroft and MerminL ll976V Thermal conductivity of 
glasses follows the universal ( x law, as was discov- 
ered bv lZeller and Pohl l)l97lj) . which dependence is ap- 
proximately followed by non-crystalline materials in gen- 
eral (FobcU, 1996). These laws are to be contrasted 
to cx T thermal conductivities of pure normal metals 
(Wicdemann-Franz law, c.f., below Eq. H^yfl). Despite 
the rapid weakening of thermal conductivity toward low 
temperatures, the dielectric materials in crystalline bulk 
are relatively good thermal conductors. One important 
observation here is that the absolute value of the bulk 
thermal conductivity in clean crystalline insulators at 
low temperatures does not provide the full basis of ther- 
mal analysis without a proper knowledge of the geom- 
etry of the structure, because the mean free paths of- 
ten exceed the dimensions of the structures. For ex- 
ample, in pur e silicon crystals, measu red at sub-kelvin 
temperatures (jKlitsner and Poiilll987tl . thermal conduc- 
tivity K ~ 10 Wm.-^K-'^T^, heat capacity C ~ 0.6 
JK^"'m~^T^ and velocity vs — 5700 m/s imply by Eq. 
(|49|l a mean free path of ~ 10 mm, which is more than 
an order of magnitude longer than the thickness of a typ- 
ical silicon wafer. Therefore phonons tend to propagate 
ballistically in silicon substrates. 

What makes things even more interesting, but at the 
same time more complex, e.g., in terms of practical 
thermal design, is that at sub-kelvin temperatures the 
dominant thermal wavelength of the phonons, Xph — 
hvs/k^T, is of the order of 0.1 /im, and it can exceed 
1 /im at the low temperature end of a typical experiment 
(see discussion in Subs. llLC.2p . A direct consequence of 
this fact is that the phonon systems in mesoscopic sam- 
ples cannot typically be treated as three-dimensional, but 
the sub-wavelength dimensions determine the actual di- 
mensionality of the phonon gas. Metallic thin films and 
narrow thin film wires, but also thin dielectric films and 
wires are to be treated with constraints due to the con- 
finement of phonons in reduced dimensions. 

The issue of how thermal conductivity and heat ca- 
pacity of thin membranes and wires get modified due to 
geometrical constraints on the scale of the thermal wave- 
length of phono ns has been ad d ressed by several authors 
experimentally ijHolmes et al\ . Il998t iLeivo and Pekolal 
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FIG. 8 The suspended silicon nitride b ridge structure on the 
left, which was used bv lSchwab et ali |2000) to measure the 
quantum of thermal conductance shown in the graph on the 
right: thermal conductance levels off at the value 16 go at 
temp eratures well below 1 K. Adapted from iSchwab et all 



systems in the two subsystems cannot be considered as 
independent. We are not aware of direct experimental 
investigations on this problem. 

On the device level reduced dimensions can be ben- 
eficial, e.g., in isolating thermally those parts of the 
devices to be refrigerated from those of the surround- 
ing heat bath. This has been the method by which 
NIS based phonon coolers, i.e., refrigerators of the lat- 
tice have been realized experimentally utilizing thin sil- 
icon nitride films and narrow silicon nitride bridges 
(IClark et all 12005: Fisher et ai, 1999; Luukancn et all 
120001: I Manninen et all Il997|l . These devices are dis- 
cussed in detail in Section V.C. 

Detectors utilizing phonon engineering are discussed in 
Sec.Hvl 



ll998HWoodcraft et all\200(t and theor etically (see, e.g., 
l)Anghel et a/.l . ll998HKuhn et alU2004} ). The main con- 
clusion is that structures with one or two dimensions 
d < Xph restrict the propagation of (ballistic) phonons 
into the remaining "large" dimensions, and thereby re- 
duce the magnitude of the corresponding quantities k and 
C at typical sub-kelvin temperatures, but at the same 
time the temperature dependences get weaker. 

In the limit of narrow short wires at low tempera- 
tures the phonon thermal conductance g ets quantized, 
as wa s experimentally demonstrated by ISchwab et al\ 
1I2OOOI). This limit had been theoretically addresse d by 
Angelescu all l)l998|) . lRego and Kirczenowl l)l998f) and 
Blencow^lfTgg^ ■ somewhat in analogy to the well-known 
Landauer result on electrical c onduction t hrough quasi- 
one-dimensional constrictions ^Landaued [T957). The 
quantum of thermal conductance is go = tt'^ k^T / ('ih) , 
and there are four phonon modes at low temperature 
due to four mechanical degrees of freedom each adding 
go to the thermal conductance of the quantum wire. In 
an experiment, see Fig. Ofour such wires in parallel thus 
carried heat with conductance g — 16.go. There are re- 
markable differences, however, in this result as compared 
to the electrical quantized conductance. In the thermal 
case, only one quantized level of conductance. Ago per 
wire, could be observed, and since the quantity trans- 
ported is energy, the " quantum" of thermal conductance 
carries cx T in its expression besides the constants of na- 
ture. 

The thermal boundary resistance (Kapitza resistance, 
after the Russian physicist P. Kapitza) between two 
bulk materials is cx due to acoustic mismatch 

l)Lounasmaal Il97^ . A remaining open issue is the ques- 
tion of thermal boundary resistance in a structure where 
at least one of the phonon baths facing each other is re- 
stricted, such that thermal phonons perpendicular to the 
interface do not exist in this particular subsystem. Clas- 
sically the penetrating phonons need, however, to be per- 
pendicul ar enough in o rder to avoid total reflection at the 
surface l|Pobelill996fl . In practice though, the phonon 



H. Heat transport in a metallic reservoir 

Energy dissipated per unit time in a biased mesoscopic 
structure is given quite generally as Q = IV ^ where / is 
the current through and V is the voltage across the de- 
vice. This power is often so large that its influence on 
the thermal budget has to be carefully considered when 
designing a circuit on a chip. For instance, the NIS cooler 
of Sec. IV.C.ll has a coefflcient of performance given by 
Eq. H77|) . with a typical value in the range 0.01... 0.1. 
This simply means that the total power dissipated is 10 
to 100 times higher than the net power one evacuates 
from the system of normal electrons. Yet this tiny frac- 
tion of the dissipated power is enough to cool the elec- 
tron system far below the lattice temperature due to the 
weakness of electron-phonon coupling. This observation 
implies that 10 to 100 times higher dissipated power out- 
side the normal island tends to overheat the connecting 
electrode significantly, again because of the weakness of 
the electron-phonon coupling. Therefore it is vitally im- 
portant to make an effort to thermalize the connecting 
reservoirs to the surrounding thermal bath efficiently. In 
the case of normal-metal reservoirs, heat can be conve- 
niently conducted along the electron gas to an electrode 
with a large volume in which electrons can then cool via 
electron-phonon relaxation. In the case of a supercon- 
ducting reservoir, e.g., in a NIS refrigerator, the situation 
is more problematic because of the very weak thermal 
conductivity at temperatures well below the transition 
temperature Tc. In this case the superconducting reser- 
voirs should either be especially thick, or they should 
be attached to normal metal conductors (" quasiparticle 
traps") as near as possible to the source of dissipation 
(see discussion in Sec. IV.C.1|) . The latter approach is, 
however, not always welcome, because, especially in the 
case of a good metallic contact between the two conduc- 
tors, the operation of the device itself can be harmfully 
affected by the superconducting proximity effect. 

Let us consider heat transport in a normal metal reser- 
voir. In the first example we approximate the reservoir 
geometry by a semi-circle, connected to a biased sample 



17 



with a hot spot of radius ro at its origin (see the inset 
of Fig.|^(b)). This hot spot can approximate, for exam- 
ple, a tunnel junction of area Trrp. The results depend 
only logarithmically on tq, and therefore its exact value 
is irrelevant when making estimates. We first assume 
that the electrons carry the heat away with negligible 
coupling to the lattice up to a distance r. According to 
Eq. (|23|l , we can then write the radial flux of heat in the 
quasiequilibrium limit in the form 



Q{r) = -KSdT/dr, 



(50) 



where k is the electronic thermal conductivity, S — ■nrt 
is the conduction area at distance r in a film of thick- 
ness t, and T is the temperature at radius r. According 
to the Wiedemann-Franz law and the temperature inde- 
pendent residual electrical resistivity in metals, one has 
n = LoaT = kqT (see below Eq. With these as- 

sumptions, using Eq. (|^ . one finds a radial distribution 
of temperature 



Tin) = \T{r^) 



QRu 

ttLci 



ln(r2/ri). 



(51) 



where ri and r2 are two distances from the hot spot, 
and we defined the square resistance as R\j ~ p/t. Thus 
making the reservoir thicker helps to thermalize it. The 
model above is strictly appropriate in the case where a 
thin film in form of a semi-circle is connected to a perfect 
thermal sink at its perimeter (at r = r2). A more ade- 
quate model in a typical experimental case is obtained 
by assuming a uniform semi-infinite film connected at 
its side to a hot spot as above, but now assuming that 
the film thermalises via electron-phonon coupling. Using 
Eq. (|50|l and energy conservation one then obtains (see 

secnrni 



dQ 

dr 



i:nrt{TP - TP) = 0, 



(52) 



where To is the lattice temperature, and p is the exponent 
of electron phonon relaxation, typically p — 5. We can 
write this equation into a dimensionless form 



d^ 



1 du 
pdp 



= uP'^ 



1. 



(53) 



Here we have defined u = {T{r)/To)'^ and p = r/rg, 
where rg = ■y/Ko/(2S)/rQ^^ ^ is the length scale of tem- 
perature over which it relaxes towards Tq. Figure El (a) 
shows the solutions of Eq. (|53|) for different values (1.1, 
1.2, 1.5, 2.0, 2.5, and 3.0) of relative temperature rise 
T(ro)/To. We see, indeed, that rg determines the relax- 
ation length. To have a concrete example let us consider 
a copper film with thickness t = 200 nm. For copper, 
using p = 5, we have kq ~ 1 WK~^m~^ and S ~ 2 • 10^ 
WK^''^m^'^, which leads to a healing length of rg ~ 500 
/im at the bath temperature of Tq = 100 mK. The tem- 
perature rise versus input power Q has been plotted in 
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FIG. 9 (Color in online edition): Temperature rise according 
to the presented reservoir heating model. For details see text. 



Fig. El (b) for copper films with different thicknesses. In 
particular for t — 200 nm, we obtain a linear response of 
(r(ro, Q)/To - 1)/Q ~ 7 • 10-VpW. 

A superconducting reservoir, which is a necessity in 
some devices, poses a much more serious overheating 
problem. Heat is transported only by the unpaired 
electrons whose number is decreasing exponentially as 
oc exp(— A/fce?^) towards low temperatures. Therefore 
the electronic thermal transport is reduced by approxi- 
mately the same factor, as compared to the correspond- 
ing normal metal reservoir. Theoretically then k is sup- 
pressed by nine orders of magnitude from the normal 
state value in aluminium at T 100 mK. It is obvious 
that in this case other thermal conduction channels, like 
electron-phonon relaxation, become relevant, but this is 
a serious problem in any case. At higher temperatures, 
say at T/Tq > 0.3, a s ignificantly t hicker superconduct- 
ing reservoir can help llClark et oiJ . l2004D . 



III. THERMOMETRY ON MESOSCOPIC SCALE 

Any quantity that changes with temperature can in 
principle be used as a thermometer. Yet the usefulness 
of a particular thermometric quantity in each applica- 
tion is determined by how well it satisfies a number of 
other criteria. These include, with a weighting factor that 
depends on the particular application: wide operation 
range with simple and monotonic dependence on temper- 
ature, low self-heating, fast response and measurement 
time, ease of operation, immunity to external parame- 
ters, in particular to magnetic field, small size and small 
thermal mass. One further important issue in thermome- 
try in general terms is the classification of thermometers 
into primary thermometers, i.e., those that provide the 
absolute temperature without calibration, and into sec- 
ondary thermometers, which need a calibration at least 
at one known temperature. Primary thermometers are 
rare, they are typically difficult to operate, but neverthe- 
less they are very valuable, e.g., in calibrating the sec- 
ondary thermometers. The latter ones are often easier to 
operate and thereby more common in research laborato- 
ries and in industry. 
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In this review we discuss a few mesoscopic thermome- 
ters that can be used at cryogenic temperatures. Ex- 
cellent and thorough reviews of general purpose cryo- 
genic thermometers, other than mesoscopic ones, can be 
found in many text books and review article s, see, e.g., 
llT.oimasmaaL ll 974t iPobellj ll flflfit iQuinnl ll flS.I) and many 
references therein. 

Modern micro- and nanolithography allows for new 
thermometer concepts and realizations where sensors can 
be very small, thermal relaxation times are typically 
short, but which generally do not allow except very tiny 
amounts of self-heating. The heat flux between electrons 
and phonons gets extremely weak at low temperatures 
whereby electrons decouple thermally from the lattice 
typically at sub-100 mK temperatures, unless special care 
is taken to avoid this. Therefore, especially at these low 
temperatures the lattice temperature and the electron 
temperature measured by such thermometers often de- 
viate from each other. An important example of this is 
the electron temperature in NIS electron coolers to be 
discussed in Section V. 

A typical electron thermometer relies on a fairly easily 
and accurately measurable quantity M that is related to 
the electron energy distribution function /(e) via 

dek{e)g[f{e)]. (54) 

-oo 

Here the kernel ^(e) describes the process which is used 
to measure /(e) and g[f] is a functional of /(e). The 
quantity M typically refers to an average current or volt- 
age, in which case g[/] = / — /o is a linear function of 
/(e) with some constant function /□; or it can refer to 
the noise power, in which case g[f] is quadratic in /. For 
the thermometer to be easy to calibrate, A:(e) should be a 
simple function dependent only on a few parameters that 
need to be calibrated. Moreover, if fc(e) has sufficiently 
sharp features, it can be used to measure the shape of 
/(e) also in the nonequilibrium limit. 



A. Hybrid junctions 

Tunnelling characteristics through a barrier separating 
two conductors with non-equal densities of states (DOSs) 
are usually temperature dependent. The barrier B may 
be a solid insulating layer (I), a Schottky barrier formed 
between a semiconductor and a metal (Sc), a vacuum 
gap (I), or a normal metal weak link (N). We are going 
to discuss thermometers based on tunnelling in a CBC 
structure. C and C stand for a normal metal (N), a su- 
perconductor (S), or a semiconductor (Sm). As it turns 
out, the current- voltage (/ — V) characteristics of the 
simplest combination, i.e., of a NIN tunnel junction, ex- 
hibit no temperature dependence in the limit of a very 
high tunnel barrier. Yet NIN junctions form elements of 
presently actively investigated thermometers (Coulomb 
blockade thermometer and shot noise thermometer) to 



be discussed separately. The NIN junction based ther- 
mometers are suitable for general purpose thermometry, 
because their characteristics are typically not sensitive 
to external magnetic fields. Superconductor based junc- 
tions are, on the contrary, normally extremely sensitive 
to magnetic fields, and therefore they are suitable only 
in experiments where external fields can be avoided or at 
least accurately controlled. 

In SBS' junctions one has to distinguish between two 
tunnelling mechanisms, Cooper pair tunnelling (Joseph- 
son effect) and quasiparticle tunnelling. The former oc- 
curs at low bias voltage and temperature, whereas the 
latter is enhanced by increased temperature and bias 
voltage. In the beginning of this section we discuss quasi- 
particle tunnelling only. 

Let us consider tunnelling between two normal-metal 
conductors through an insulating barrier. I-V charac- 
teristics of such a junction were given by Eq. H29I) . We 
assume quasi-equilibrium with temperatures Ti, i — L, R 
on the two sides of the barrier. Since Ni{E) = 1 to 
high precision at all relevant energies {k-^Ti, eV ^ Ep = 
Fermi energy), Eq. H29|) can be integrated easily to yield 
/ ~ V/Rt- Therefore the I—V characteristics are ohmic, 
and they do not depend on temperature, and a NIN junc- 
tion appears to be unsuitable for thermometry. There is, 
however, a weak correction to th is result, due to the finite 
height of the tunnel barrier (fSimmong| . ri963a|) . which will 
be discussed in subsection III.D. 



1. NIS thermometer 

As a first example of an on-chip thermometer, let us 
discuss a tunnel junction betwe en a normal metal and a 
superconductor (NIS junction) ("Ro well and TsuiL Il976|) . 
I-V characteristics of a NIS junction have the very im- 
portant property that they depend on the temperature 
of the N electrode only, which is easily verified, e.g., 
by writing I{V) of Eq. with Nl{E) = const, and 
Nfi{E) — Ns{E) in a symmetric form 

= / NsiE)[hiE-eV)~hiE + eV)]dE. 

(55) 

This insensitivity to the temperature of the supercon- 
ductor holds naturally only up to temperatures where 
the superconducting energy gap can be assumed to have 
its zero-temperature value. This is true in practice up to 
T/Tc < 0.4. 

Employing Eq. H55|) one finds that a measurement of 
voltage V at a constant current I — Iq yields a direct 
measure of /n(-E), and in quasi-equilibrium, where the 
distribution follows thermal Fermi-Dirac distribution, it 
also yields temperature in principle without fit parame- 
ters. Figure [TUI (a) shows the calculated I-V characteris- 
tics of a NIS junction at a few temperatures T/T^- Figure 
IIUI (b) gives the corresponding thermometer calibration: 
junction voltage has been plotted against temperature at 
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FIG. 10 NIS thermometer characteristics. Calculated I-V 
characteristics of a normal metal - superconductor tunnel 
junction at various relative temperatures T/Tc in (a). The 
corresponding voltage over the junction as a function of tem- 
perature when the junction is biased at a constant current is 
shown in (b). Characteristics at a few values of current are 
shown. 



a few values of (constant) measuring current. 

A NIS thermometer has a number of features which 
make it attractive in applications. The sensing ele- 
ment can be made very small and thereby NIS junctions 
can probe temperature locally and detect temperature 
gradients. Junctions made by electron beam lithogra- 
phy can be much smaller than 1 /itm in linear dimen- 
sion (|Nahum and MartiniaL Il99,'^ . Using a scanning tun- 
nelling microscope with a superconducting tip as a NIS 
junction one can most likely probe the temperature of the 
surface l ocally on nanometer s cales in an instrum ent like 
those of lMoussv ef^ l|200ll) : IVinet et al\ l|200lh . Self- 
heating can be made very small by operating in the sub- 
gap voltage range e\V\ < A (see Fig. llOf) . where current 
is very small. The superconducting probe is thermally 
decoupled from the normal region whose temperature is 
monitored by it. The drawbacks of this technique in- 
clude high sensitivity to external magnetic field, high 
impedance of the sensor especially at low temperatures, 
and sample-to-sample deviations from the ideal theoreti- 
cal behaviour. Due to this last reason, a NIS junction can 
hardly be considered as a primary thermometer: devia- 
tions arise especially at low temperatures, one prominent 
problem is saturation due to subgap leakage due to non- 
zero DOS within the gap and Andreev reflection. 

A fast version of a NI S thermometer was implemented 
bv lSchmidt et al\ ^200^ . They achieved sub -fis readout 
times (bandwidth about 100 MHz) by imbedding the NIS 
junction in an LC resonant circuit. This rf-NIS read-out 
is possibly very helpful in studying thermal relaxation 
rates in metals, and in fast far-infrared bolometry. 

NIS jun ction thermometry h as been applied in x-ray 
detectors llNahum et oil Il993l). far-infrared bolometers 
l|Chouvaey et adll999HMees et a0 . ll99l|) , in probing the 
energy distribution of electrons in a metal (|Pothier et all 
Il997al) , and as a therm ometer in electronic coolers at sub- 
kelvin temperatures ijLeivo et al 1 Il996t iNahum et all . 



Il994j) . It has also been suggested to be used as a far- 
infrared photon counter ijAnghel and Kuzminl |200^ . In 
many of the application fields of a NIS thermometer it is 
not the only choice: for example, a superconducting tran- 
sition edge sensor (TES) can be used very conveniently in 
the bolometry and calorimetry applications (see Sec. IV). 

Recently Schottky contacts between silicon and su- 
perconducting metal have been shown to exhibit sim- 
ilar characteristics a s fully metallic NIS junctions 
(|Buonomo et all 1200,'^ ISavin et all l200ll) . These struc- 
tures are discussed in subsection V.C.3. 



2. SIS thermometer 

A tunnel junction between two superconductors 
supports supercurrent, whose critical value Ic has 
a magnitud e which depends on t e mpera ture ac- 
cording to llAmbegaokar and Baratofil Ll963): Ic — 
tanh( ^^^ ). This can naturally be used to indicate 
temperature 1)ecause of the temperature dependence of 
the energy gap and the explicit hyperbolic dependence 
on T. Yet, these dependencies are exponentially weak at 
low temperatures. Another possibility is to suppress Ic 
by magnetic field, e.g., in a SQUID configuration, and to 
work at non-zero bias voltage and measure the quasipar- 
ticle current, which can be estimated by Eq. H29|l again 
with both DOSs given by Eq. H12() now. The resulting 
current depends approximately exponentially on temper- 
ature, but the favourable fact is that the absolute magni- 
tude of the current is increased because it is proportional 
to the product of the two almost infinite DOSs matching 
at low bias voltages, and, in practical terms, the method 
as a probe of a quasiparticle distribution is more robust 
against magnetic noise. Figure ITTl shows the calculated 
and measured dependences of /(l^) at a few values of tem- 
perature T/Tc- Although the curves corresponding to the 
lowest temperatures almost overlap here, it is straightfor- 
ward to verify that a standard measurement of current 
can resolve temperatures down to below O.lTc using an 
ordinary tunnel junction. 

SIS junctions have not found much use as thermome- 
ters in the traditional sense, but they are extensively 
investigated and used as photon and particle detec- 
tors because of th eir high energy resolution (see, e.g., 
l|Booth and Goldiel il996'l and Section IV). In this con- 
text SIS detectors are called ST J detectors, i.e., super- 
conducting tunnel junction detectors. Another applica- 
tion of SIS structures is their use as mixers (TinkhanJ 
Il996|) . While superficially similar to ST J detectors, their 
theory and operation differ significantly from each other. 
SIS junctions are also suitable for studies of quasipar- 
ticle dynamic s and fluctuations in general (see, e.g., 
llWilson et a/.Ll2n0lh 1. 
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FIG. 11 SIS junction as a thermometer, (a) Calculated and 
(b) measured I — V curves of a SIS tunnel junction at a few 
temperatures. Su percurrent has been sup pressed. Experi- 
mental data from (|Savin and Pekolai l200.'Tt) . 



3. Proximity effect thermometry 

For applications requiring low-impedance (< 1 fi) ther- 
mometers at sub-micron scale, the use of clean NS con- 
tacts may be more preferable than thermometers apply- 
ing tunnel contacts. In an SNS system, one may again 
employ either the supercurrent or quasiparticle current as 
the thermometer. For a given phase (j) between the two 
superconductors, th e former can be expressed through 
l|Belzig et all Il999|) 
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dEjs{E-<j>){l--2f{E)), (56) 



where Rn '^s the normal-state resistance o f the weak link , 
js(E;<j)) is the spectral supercurrent l)Heikkila et all 
12002), and the energies are measured from the chemi- 
cal potential of the superconductors. Hence, the super- 
current has the form of Eq. H54|l . In practice, one does 
not necessarily measure Is{(l>), but the critical current 
Ic = max0 Is ((/)). In diffusive junctions, this is obtained 
typically for (j) near 7r/2, although the maximum point 
depends slightly on temperature. 

The problem in SNS thermometry is in the fact that 
the supercurrent spectrum js(£) depends on the qual- 
ity of the interface, on the specific geometry of the sys- 
tem, and most importantly, on the distance L between 
the two superconductors com pared to the superconduct- 
ing c oherence length = -^//iD/ (2 A) (fHeikkila et all 
|2002|) . Therefore, Ic{T) dependence' is not universal. 
However, the size of the junction can be tuned to meet 
the specific temperature rang e of interest. In the limit 
of sh ort junctions, L < l)Kulik and Omel'vanchukl 
Il978|) . the temperature scale for the critical current is 
given by the superconducting energy gap A and for 
ksT <^ A, the supercurrent depends very weakly on 
the temperature. In a typical case ~ 100 . . . 200 
nm, and thus already a weak link with L of the order 
of 1 /xm, easily realisable by standard lithography tech- 
niques, lies in the "long" limit. There, the critical cur- 
rent is Ic = c(fcBr)3/2 exp{-.y2TTkBT/ET)/{eRNVE^) 



ijZaikin and Zharkovl Il98l[l . This equation is valid for 
ksT > 5Et- Here Et = hP/L-^ and the pr efactor c 
depends on the geometry l|Heikkila et a/.Ll2'002^ . for ex- 
ample for a two-probe configuration c — 64\/27r^/(3 -I- 
2V2). The exponential temperature dependence and 
the crossover between the long- and short-junction limits 
were experimentally investigated bv lDubos et al\ l)200lj) 
and the above theoretical predictions were confirmed. 

Using a four- p robe configuration w ith SNS junctions of 
different lengths. Fjiang et all l|2003f) exploited the strong 
temperature dependence of the supercurrent to measure 
the local temperature. The device worked in the regime 
where part of the junctions were in the supercurrent- 
carrying state and part of them in the dissipative state. 
In the limit where supercurrent is completely suppressed, 
there is still a weaker temperature dependence of the 
conductan ce, due to the proxirn ity-effect correction (c.f., 
Eq. gni) (fCharlat et all Il996l). Th i s can also be used 
for thermometry l)Aumentado et all Il999|) . but due to 
the much smaller effect of temperature on the conduc- 
tance, it is less sensitive. 

Besides being just a thermometer of the electron 
gas, critical current of a SNS Josephson junction has 
been sho wn to p robe the electron energy distribution 
Easelmans et a/j. [1999. 2001b: .Giazotto e t oL .2004bt 
iHuang et all l2002() as indicated by Eq. 



B. Coulomb blockade thermometer, CBT 

Single electron tunnelling (SET) effects were foreseen 
in micro-lithographic s tructu res in the middle of 1980's 
l)Averin and Likharevl Il986j) . Such effects in granular 
stru ctures had been known to exist already much ear- 
lier llGiaever and Ze llcr"l968 ULambe and Jakleviciri969t 
iNeugebauer and Webb, .196^ The first lithographic 
SET device was demonstrated by iFulton and Polar] 
(|l987f) . Since that time SET effects have formed a 
strong sub field in mesoscopic physics. Typically SET 
devices operate in the full Coulomb blockade regime, 
where temperature is so low that its influence on elec- 
trical transport characteristics can be neglected; at low 
bias voltages current is blocked due to the charging en- 
ergy of single electrons. Coulomb blockade thermome- 
ter (CBT) operates in a different regime, where single 
electron effects still play a role but temperature pre- 
dominant ly influences t h e electrical transport charac- 
teristics l|Bergsten et all I2OOII: iFarhangfar et all Il997t 
iMeschke et a^J . 120041: iPekola et a^J .|l994>l. CBT is a pri- 
mary thermometer, whose operation is based on compe- 
tition between thermal energy k^T at temperature T, 
electrostatic energy eV at bias voltage V, and charging 
energy due to extra or missing individual electrons with 
unit of charging energy Ec — e^/(2C*), where C* is the 
effective capacitance of the system. Figure shows an 
SEM micrograph of a part of a typical CBT sensor suit- 
able for the temperature range 0.02. ..IK. This sensor 
consists of four parallel one-dimensional arrays with 40 
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FIG. 12 A typical CBT sensor for the temperature range 20 
mK - IK. The structure has been fabricated by electron beam 
lithography, combined with aluminium and copper vacuum 
evaporation. Both top view and a view at an oblique angle 
are shown; the scale indicated refer s to the top view. Figure 
adapted from llMeschke et a/.ll200^ . 



junctions in each. 

In the partial Coulomb blockade regime the I—V char- 
acteristics of a CBT array with N junctions in series do 
not display a sharp Coulomb blockade gap, but, instead, 
they are smeared over a bias range eV ~ Nk-gT . The 
asymptotes of the I — V sA large positive and negative 
voltages have, however, the same offsets as at low T, de- 
termined by the Coulomb gap. In the partial Coulomb 
blockade regime it is convenient to measure not the I—V 
directly, but the differential conductance, i.e., the slope of 
the I—V curve, G = dl/dV, vs. V. The result is a nearly 
bell shaped dip in conductance around zero bias. Figure 
1131 illustrates a measured conductance curve, scaled by 
asymptotic conductance Gt at large positive and nega- 
tive voltages. The important property of this character- 
istic is that the full width of this dip at half minimum, 
Vi/2i approximately proportional to T, determines the 
temperature without any fit, or any material or geome- 
try dependent parameters, i.e., it is a primary measure 
of temperature. In the lowest order in -^p, one finds for 
the symmet ric linear array of N junctions of capacitanc e 
C in series llFarhangfar et all Il997t iPekola et all Il994|) 

For such an array C* = NC/[2{N - 1)] and g{x) = 
[a;sinh(x) - 4sinh^(x/2)]/[8sinh'^(x/2)]. The full width 
at half minimum of the conductance dip has the value 



Vi/2 5.439iVfcBr/e, 



(58) 



which allows one to determine T without calibration. 

It is often convenient to make use of the secondary 
mode of CBT, in which one measures the depth of the 
dip, AG/Gt, which is proportional to the inverse tem- 
perature, again in the lowest order in 



AG/Gt ^Ec/{6kBT). 



(59) 



Measuring the inflection point of the conductance 
curve provides an al ternative means to det ermine T with- 
out calibration. In ijBeresten et a/.Ll200lj) this technique 





FIG. 13 Performance of a CBT thermometer, (a) Typical 
measurement: G{V) /Gt is the differential conductance scaled 
by its asymptotic value at large positive and negative volt- 
ages, plotted as a function of bias voltage V. V1/2 indicates 
the full width at half minimum of the characteristics, which 
is the main thermometric parameter. The full depth of the 
line, AG/Gt, is another parameter to determine tempera- 
ture. In (b) temperature deduced by CBT has been com- 
pared to that obtained by a "^He melting curve thermometer. 
Saturation of CBT below 20 mK indicates typical thermal 
decoupling between e lectrons and phonons. Figure from Ref. 
llMeschke all\2004) . 



was used for fast thermometry on two dimensional junc- 
tion arrays. A simple and fast alternative measurement 
of CBT can be achieved by detecting the third har- 
monic current i n a pu re AC voltage biased configuration 
llMeschke et qiJ . 120051) . 

There are corrections to results l|57)l . l|55|) and ifS^ 
due to both next order terms in £*<-;/ /crT and due t o 
non-uniformities in the array ijFarhangfar et o/l Il997tl . 
Higher order corrections do sustain the primary nature 
of CBT, which implies wider temperature range of op- 
eration. The primary nature of the conductance curve 
exhibited by Eq. H57|l through the bias dependence 
d^ Nksr ) preserved also irrespective of the capacitances 
of the array, even if they would not be equal. Only a 
non-uniform distribution of tunnel resistances leads to a 
deviation from the basic result of bias dependence of Eq. 
H57|l . Fortunately the influence of the inhomogeneity on 
temperature reading is quite weak, and it leads typically 
to less than 1 % systematic error in T. 

The useful temperature range of a CBT array is limited 
at high temperature by the vanishing signal (AG/Gt 
T^^). In practice the dip must be deeper than ~ 0.1 % 
to be resolvable from the background. The low temper- 
ature end of the useful temperature range is set by the 
appearance of charge sensitivity of the device, i.e., the 
background charges start to influence the characteristics 
of the thermometer. Th is happens when AG/Gt — 0.5 
l|Farhangfar et gilll997^ . With these conditions it is ob- 
vious that the dynamic range of one CBT sensor spans 
about two decades in temperature. 

The absolute temperature range of CBT techniques 
is presently mainly limited by materials issues. At 
high temperatures, the measuring bias range gets wider, 
Vi/2 oc T, and the conductance becomes bias dependent 
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also due to the finite (energy) height of t he tunnel bar- 
rier (|Gloos et Mi l2(m i immon Because of 
this, the present high temperature limit of CBTs is sev- 
eral tens of K. The (absolute) low temperature end of the 
GET technique is determined by self-heating due to bi- 
asing the device. One can measure temperatures reliably 
down to about 20 mK at present with better than 1 % 
absolute accuracy i n the range 0.05. . . 4 K, and about 3 
% down to 20 mK l|Meschke a^J.l2004j) . 

Immunity to magnetic field is a desired but rare prop- 
erty among thermometers. For instance, resistance ther- 
mometers have usually strong magnetoresistance. There- 
fore one could expect CBT to be also " magnetoresis- 
tive" . On the contrary, CBT has proven to be im- 
mune to even the strongest magnetic fields (> 20 T) 
llvan der Linden and BehniaL l2004 iPekola et all l2002l 
Il998|) . This happens because CBT operation is based on 
electrostatic properties: tunnelling rates are determined 
by charging energies. Hence, CBT should be perfectly 
immune to magnetic field as long as energies of electrons 
with up and down spins do not split appreciably, which 
is always the case in experiments. 

Coulomb blockade is also a suitable probe of non- 
equilibriu m energy distributions as demonstrated by 
lAnthore et al. ( 2003,) . 



C. Shot noise thermometer, SNT 

Noise current or voltage of a resistor (resistance R) has 
been known to yield absolut e tempera t ure si n ce 1920's 
llBlanter and BiittikeJ. 120001: iJohnsonl 119281: iNvauistl 
Il928|) . Till very recently only equilibrium noise, i.e., 
noise across an unbiased resistor, was employed to mea- 
sure temperature. In this case Johnson noise voltage 
and current in squared have expectation values (v^) — 
R^{in) — Ak^TRl\f , where A/ is the frequency band of 
the measurement. In other words, the current spectral 
density is given by = Ak^T / R, and voltage spectral 
density by S\j = AksTR. There are several critical issues 
in measuring temperature through Johnson noise. First, 
the bandwidth has to be known exactly. Secondly, the 
small signal has to be amplified, and the gain must be 
known to high accuracy. Thirdly, the noise signal be- 
comes extremely small at cryogenic temperatures, which 
in turn means that the measurement has to be able to de- 
tect a very small current or voltage, and at the same time 
no other noise sources should cause comparable voltages 
or currents. 

Success in noise thermometry depends critically on the 
performance of the amplifiers used to detect the tiny 
noise current (or voltage). Constant progress in improv- 
ing SQUIDs (Superconducting Quantum Interference De- 
vice) and in optimizing their operation has pushed the 
lowest temperature measurable by a Johnson noise ther- 
mometer down to below 1 mK (jLusher et all l200l|) . 

It is not only the equilibrium noise that can be useful 
in thermometry. In the opposite limit, at eV ^ fceT, the 
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FIG. 14 Cross-over from Johnson noise to shot noise when 
increasing the bias voltage of a tunnel junction. The data 
shown by open symbols has been measured at 300 K with 
Rt ~ 0.32 GQ, and data shown by solid symbols at 77 K 
with Rt = 2.7 GQ. Adapted from llBirk et aZ.lll995ri . 



dominating noise, e.g. , in tunnel junctions, is shot noise 
l)Blanter and Biittikeil.l2000l:ISchottk^ll918(l . whose cur- 
rent spectral density is given asymptotically by Sj = 
F2e\I\. Fano factor F equals 1 for a tunnel junction. 
The way Johnson noise transforms into shot noise upon 
increasing the bias obeys the relation (see Eq. H34a|) with 
Tn 0, and, e.g., l)van der Zi el. 1986)) 

5,(T) = 2e/coth(^). (60) 

At \V\ :s> k^T/e, Eq. ifHU]! obtains the Poisson expression 
above, whereas at \V\ <C k-QT/e it assumes the thermal 
Johnson form Si ~ Ak-oT/RT- 

Temperature dependent cross-over characteristics from 
Johnson to shot noise according to Eq. ifBUIl have been 
demonstrated experiment ally, e.g., in scann ing tunnelling 
microscope experiments d Birk etlj,i as shown in 

Fig. Recently, iSpietz et al\ l)2003) employed this 

cross-over in a lithographic tunnel junction between 
metallic films for thermometry (shot noise thermometer, 
SNT) in the temperature range from few tens of mK up 
to room temperature. At both ends of this range there 
are systematic errors in the reading, which can possi- 
bly be corrected by a more careful design of the sensor. 
In a rather wide range around 1 K (about 0.1. . . 10 K), 
the absolute accuracy is better than 1 %. The crossover 
voltage depends only on k^T/e, which means that the 
thermometer is indeed primary; ideally the temperature 
reading does not depend on the materials or geometry 
of the sensor . Figu re 1151 shows the experimental data of 
ISpietz et all l|2003|) at several temperatures: normalised 
(by Johnson noise) current noise has been plotted against 
normalised voltage x = eV/ (2fcBT), whereby all curves at 
different temperatures should lie at ^J^™™ = xcoth(2:). 
As seen in the bottom half of the figure, the residuals are 
smallest around 1 K. 

As discussed above, noise measurements are difficult to 
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FIG. 15 Normalised junction noise versus normalised voltage 
at various temperatures. The residuals from the expected 
X cothfa:) law are sh own in the bottom half of the figure. From 
llSpietz et aLll2003l) . 



FIG. 16 The suspe nded thermal sensor employed in 
ijBourgeois et aZ.! l200Fi) . The Nb^^'Ni.j./ thermometer can be 
seen in the lower part of the rectangular silicon membrane. 
The 450 000 Al superconducting rings are located in the mid- 
dle part of the membr ane; examples of them are shown in (b) 
and (c). Figure from (iBourgeois et 0^.1120051) . 



perform especially at low temperature where signal gets 
very small. At any temperature it is critically important 
to know the frequency window of the measurement and 
the gains in the amplifiers in Johnson noise thermometry. 
Yet the SNT avoids some of these problems. It is based 
on the cross-over between two noise mechanisms, both 
of which represent white noise, whereby the frequency 
window is ideally not a concern, since the same readout 
system is used in all the bias regimes. Moreover, the 
gains of the amplifiers are not that critical either because 
of the same argument. One can also use relatively high 
bandwidth which increases the absolute noise signal to be 
measured, and thereby makes the measurement faster. 

The SNT technique has a few further attractive fea- 
tures. It is likely that its operation can be easily ex- 
tended up to higher temperatures despite the deviations 
observed in the first experiments. The sensor consists 
of just one, relatively large size tunnel junction, which 
means that it is easy to fabricate with high precision. 
Also it is hkely, although not yet demonstrated, that the 
SNT is not sensitive to magnetic field, since its operation 
is based on tunnelling characteristics in a NIN tunnel 
junction as in CBT. 

Finally, noise measurements can in principle be used 
to measure the distribution function in non-equilibrium 
as well, as proposed by fPistolcsi et al... (,2004,) . 

D. Thermometry based on the temperature dependent 
conductance of planar tunnel junctions 

The effect of temperature on the current across a 
tunnel barrier with finite height is a suita ble basis for 
thermometry in a wide temperature range ijGloos et all 
120001) . ISimincm^ (|l963a) showed that the tunnelling con- 
ductance at zero bias across a thin insulating barrier de- 



pends on temperature as 

G(r)-Go[l + (r/To)2], (61) 

where Gq is the temperature independent part of con- 
ductance and the scaling temperature Tq depends on the 
barrier height (jjQ. For a rectangular barrier of width s 
one has Tn = <t>o -^here m is the effective mass 
of the electrons within the insulating barrier. Experi- 
ments over a temperature range from 50 K up to 400 K 
on Al-AlOx-Al tunnel junctions have demonstrated that 
Eg. lICTll is obeyed rem arkably well l|Gloos et all l200Ct 
ISuoknuuti et all l200ll) . Moreover, in these measure- 
ments the scaling temperature was found to be Tq ~ 720 
K in all samples, without a clear dependence on the spe- 
cific (zero temperature) conductance of the barrier, which 
varied over three orders of magnitude from 3 /zS/ fim^ up 
to 3000 /iS//im^. This property makes the method at- 
tractive in wide range thermometry, and Tq can indeed 
be considered as a material specific, but geometry and 
thickness independent parameter up to a certain accu- 
racy. 

E. Anderson-insulator thin film thermometry 

As regards to temperature read-out of micro- 
calorimetric devices, resistive thin film thermometers 
near the metal-insulator transition (MIT) are rela- 
tively popular. Electrical resistivity properties on 
both sides of the MIT are rather well understood 
ijBelitz and KirkpatrickLri994j) . and in general resistance 
of such thin films shows strong temperature dependence, 
suitable for thermometry and in particular for calorime- 
try. On the insulator side resistivity p is determined by 
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hopping, and it has typically p oc eS'^°/'^^ temperature 
dependence, with Tg and n constants. On the metal- 
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have been successfully employed. The suitable conduc- 
tion regime can be tailored b y adjusting x [x') in elec- 
tron beam co-evaporation ( De nlinger et al\ . Il994(l or in 
dc magne tron sp uttering of Nb in a nitrogen atmosphere 
IFomina va et alV Qj)97) . 

Bolometric and calorimetric radiation detectors are 
discussed in detail in Sec. IIVI Here we briefly mention 
the application of a Nbj:'Ni_j./ thermometer in a mea- 
surement of the heat capacity of 450 ODD superconducting 
thin film loops on a silicon membrane ijBourgeois et al\ . 
l2005l) . see Fig.^l The heat capacity of the loops is pro- 
portional to their total mass, which was about 80 ng in 
this case. Vortices entering simultaneously into the 450 
000 loops under appHcation of ma gnetic field could b e 
observed. A similar measurement (jLindell et a^J . l2000(l . 
employing a NIS thermometer could resolve the specific 
heat jump at of 14 thin film titanium disks with total 
mass of 1 ng on a silicon nitride membrane. 



IV. THERMAL DETECTORS AND THEIR 
CHARACTERISTICS 

The absorption of electromagnetic radiation by mat- 
ter almost always ends in the situation where the inci- 
dent optical power has, possibly via cascades of different 
physical processes, transformed into aggravated random 
motion of lattice ions, i.e. into heat. Here lies the fun- 
damental principle behind the thermal detection of ra- 
diation: the transformation of the input electromagnetic 
energy to heat. In many cases, this state of maximum 
entropy has lost all the coherent properties that the in- 
cident radiation might have possessed, but yet there is 
information in this messy final state of the system: the 
rise in the system temperature. In this section, we give a 
short overview of thermal detectors, their theory and op- 
eration, and discuss some examples of thermal detectors 
and their applications. 

There is a very small difference between thermom- 
etry and bolometry, the thermal radiation detection. 
Yet thermometry often implies measurement of tem- 
perature changes over a large fractional temperature 
range, whereas in the typical bolometric application, 
the observed temperature variations are extremely small. 
This allows us in the following to concentrate on the 
limit of small temperature changes around some mean 
well-defined value, i.e., we assume the operation in the 
quasiequilibrium regime (see Sec. II). In essence, bolom- 
etry is really high-precision thermometry, nothing more, 
with the added ingredient of finding efficient ways of cou- 
pling incident radiation to the device. 

Although thermal detectors have developed enormous 



diversity since their introduction in 1880 by Samuel P. 
Langley, and major advances in the way the temperature 
rise is measured and the detectors are constructed have 
been made, the basic principle remains the same. The 
operating principle of a thermal detector can be traced 
back to the generalized thermal model, shown in Fig. ^ 
The four main parts of a thermal detector are a thermally 
isolated element, a thermal link with a thermal conduc- 
tance Gth (here Ge-ph or Gph~sub), a thermal sensing 
element (i.e., a thermometer), and a coupling structure 
(e.g., an impedance matching structure for electromag- 
netic radiation) that serves to maximize the absorbtion 
of the incident radiation, be it in the forms of alpha par- 
ticles or microwaves. 



A. Effect of operating temperature on the performance of 
thermal detectors 

Regardless of the exact architecture of the thermal 
detector, lowering the operating temperature will im- 
prove the performance significantly. This fact intro- 
duces the cooling techniques to the use of thermal de- 
tectors. Because of the diversity of the technological ap- 
proaches it is hard to summarize the effect of the op- 
erating temperature in a fully universal fashion. How- 
ever, some general trends can be evaluated. The fig- 
ure of merit for thermal detectors is the noise equivalent 
power (NEP), which relates to the signal to noise ratio 
by SNR = Qopt/(NEP-\/2Tint), where Qopt is the incident 
optical power and Tint is the post-detection integration 
time. The limiting NEP for an optimized thermal detec- 
tor is given by the thermal fiuctuation noise (TEN) aris- 
ing from random fluctuations of energy across the ther- 
mal link with a thermal conductance Gth = dQ/dT which 
result in variations in the temperature of the device (see 
also discussion in Subs. ITLK2|) . The TEN Hmited noise 
equivalent power (NEP) in the li near order is given by 
the fluctuation-dissipation result l|Matheilll982() 

oc e-p decoupling 

cx T^/^ lattice isolation 



NEPtfn « v/4fcBT2Gth 



(62) 



with the two cases relating the temperature dependence 
to the location of the thermal bottleneck. The expression 
for the lattice isolated case is generally valid at temper- 
atures T < Od/IO with 9j:i the Debye temperature of the 
insulating material. The temperature T is the highest 
of the temperatures present in the system, i.e., the tem- 
perature of the thermally isolated element or that of the 
heat sink. 

Refrigeration, combined with the advances in nano- 
lithography techniques have recently opened a whole new 
realm for the application of electron-phonon decoupling 
to improve the thermal isolation of bolometric detec- 
tors. The operation of these so-called hot electron de- 
vices is usually limited to low (< 1 K) temperatures 
and very small thermally active volumes V as the TEN 
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between the electron gas and the lattice is given by 
^5fcBSe-pV(T| + T6J (lOnlwala, f.t all where S 

varies between 1-4 nW//ini^/K^ in metals (See TablePl. 
The use of the lattice for the thermal isolation lends itself 
to operation over much broader temperature range as the 
geometry of the thermal link can be used to increase the 
thermal resistance. 

Typically, the bath temperatures of cryogenic thermal 
detectors are centered around four temperature ranges: 
4.2 K, the boiling temperature of liquid He at 1 atm, 
around 300 mK a temperature attainable with '^He sorp- 
tion refrigerators, around 100 mK, easily attainable us- 
ing a compact adiabatic demagnetization refrigerator 
(ADR), or below ~ 50 mK, when a dilution refrigerator 
is used. In the recent years, the use of electronic refriger- 
ation has become appealing with the development of the 
SINIS coolers. These coolers could enable the operation 
of thermal detectors at a much lower temperature than 
is 'apparent' to the user. For instance, the ^He sorption 
coolers are quite attracting due to their compactness, low 
cost and simple operation. Used together with a SINIS 
cooler would allow for an affordable cryogenic detector 
system without having to sacrifice in performance. 

Since the thermal fluctuation noise in a thermal detec- 
tor is proportional to -|- T^^j^, maximum performance 
improvement is achieved by the reduction of the bath 
temperature ijAnghel et all liooH) . Modest performance 
increase is possible with direct electronic cooling, such as 
is the case in many of the SINIS bolometer experiments. 

In addition to improved noise performance, direct cou- 
pling of a SINIS cooler to a thermal detector could allow 
for increasing the dynamic range of for example bolome- 
ters based on transition edge sensors (TES): A SINIS 
cooler can be used to draw a constant power from a TES 
that would otherwise saturate due to an optical load. 



B. Bolometers: Continuous excitation 

Thermal detectors that are used to detect variations 
in the incident flux of photons or particles are called 
bolometers (from the Greek word bote - beam). This 
condition is generally met when the mean time between 
incident quanta of energy that arrive at the detector is 
much shorter than the recovery time of the bolometer 
^. The bolometric operating principle is very simple: 
Change AQopt in the incident optical power creates a 
change in the temperature of a thermally isolated element 
by AT = AQopt/Gth- A sensitive thermometer is used 
to measure this temperature change. The recovery time 
To is determined by the heat capacity C of the bolome- 
ter, and the thermal conductance Gth with tq = C/Gth, 
analogously to an electrical RC circuit. 



Bolometric detectors remain popular today, 125 years 
after their first introduction. Probably the most impor- 
tant advantage of thermal detectors is their versatility: 
Bolometers can detect radiation from a -particles to radio 
waves, their dynamic range can be easily adapted for a 
variety of signal or background levels. As an extreme ex- 
ample, bolometers have been used to detect infrared radi- 
ation from nuclear fireballs (Stubbs and Phillip s, 196 0) , 
and the cosmic microwave background. ijLamarre et al\ . 
I2003D . 

In the early days, bolometers typically utilized hand- 
crafted construction (dental floss, cigarette paper and 
balsa wood are ex amples of typica l materials used in the 
construction) (Da vis et all Il964|) . If detecting electro- 
magnetic waves, typically the absorber consisted of a 
metal with a suitable thickness yielding a square resis- 
tance of 377 fl, i.e. matching the impedance of the vac- 
uum. Further improvements on matching were achieved 
by placing the bolometer in a resonant cavity. One ma- 
jor setback with bolometers in their early days was un- 
avoidably slow speed, caused by the large heat capacity 
resulting from the macroscopic size of the components 
used. The dawn of modern microfabrication techniques 
has all but eliminated this shortcoming, with bolometers 
of high sensitivity achieving time constants as short as a 
few hundred nanoseconds. 

Today, the most common type of cryogenic resistive 
bolometers utilize transition edge sensors for the ther- 
mometry. In a TES bolometer, a superconducting film 
with a critical temperature Tc is biased within its super- 
conductor - normal metal transition where small changes 
in the film temperature result to changes in the current 
through the device (or the voltage across the film). In 
most cases, the TES consists of two or more sandwiched 
superconductor-normal metal layers. The relative thick- 
nesses of the S and N layers are used to tune the tran- 
sition temperature to a desirable value by the proximity 
effect. 

Transition-edge sensors are by no means a novel type of 
a thermal detector, as first sugge stions for t heir use came 
out a lready in the late thirties ijAndrewsl llQSSt iGoetz , 
1936), and first experimental results by 1941 jAndrews . 



1941j) . Two principal problems prohibited the wide use 



The opposite (calorimetric) limit will be discussed in section llV.CI 



of this type of thermal detectors for some fifty years: 
Typically the normal-state resistance of the supercon- 
ducting films was too low in order to obtain adequate 
noise matching with field effect transistor (FET) pream- 
plifiers, and the lack of good transimpedance amplifiers 
usually required the films to be current biased with a volt- 
age readout. This introduced a requirement to tune the 
heat bath temperature very accurately within the narrow 
range of temperatures in the superconducting transition. 
This also made the devices exceedingly sensitive to small 
variations in the bath temperature, introducing stringent 
requirements for the heat bath stability. 

These limitations can be overcome by the use of an 
external negative feedback circuit that maintains the 
film within its transition temperature and above the 
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ba th temperatu re. ^Such an approach was adopted 

by IClarke et al. who were able to demonstrate 

NEP=1.7 • 10~^^ W/-\/Hz at an operating temperature 
of 1.27 K, using a transformer-coupled FET as the volt- 
age readout. Introducing a negative feedback has similar 
advantages as in the case of operational amplifiers: lin- 
earity is improved, sensitivity to internal parameters of 
the amplifier (or bolometer) is reduced, and the speed is 
increased. Interestingly, the use of an external negative 
feedback in conjunction with superconducting transition 
edge sensors never found widespread use, possibly due 
to the (slightly) more complicated read-out architecture, 
and the need for a matching transformer. 

The breakthrough of TESs came in 1995, when super- 
conducting quantum interference device (SQUID) amme- 
ters, which are inherently well suited in matching to low 
load impedances , were introduced as the readout devices 
for TESs l|lrwinLll995Hlrwin et al\ .\l995a]. This allowed 
for the use of voltage biasing, which introduces strong 
negative electrothermal feedback (ETF) that causes the 
thermally isolated film to self-regulate within its super- 
conducting transition. The local nature of the ETF 
makes the operation of these detectors very simple as 
no external regulation is necessary. As with an external 
negative feedback, an important advantage of the voltage 
biased TES is the fact that once the bath temperature 
is below ~ Tc/2, the need for bath temperature regula- 
tion is significantly relaxed. Finally, the increased speed 
due to the strong negative ETF increases the bandwidth 
of the detector, allowing for either detecting faster sig- 
nal changes in a bolometer, or higher count rates in a 
calorimeter. A comprehensive review of the theory and 
operation, of voltage biased TESs has been recently pub- 
lished ijirwin and HiltonL 120(5^ . 

The behaviour of thermal detectors is generally well 
understood. In the following discussion we shall summa- 
rize the main results for the theory of thermal detectors. 
The results are quite generally applicable to any resistive 
bolometers, but the treatment is geared towards voltage 
biased TES, firstly as they are currently the most popular 
type of thermal detectors, and secondly because the elec- 
trothermal effects are very prominent in these devices, 
having a major impact on the device speed and linearity. 

We start by writing the equation governing the ther- 
mal circuit (here we limit the discussion to a simple case 
of one thermal resistance and heat capacity). Gener- 
ally, the power flow to the heat sink is given by Qout = 
K{T'^—Tq), where if is a constant which depends on ma- 
terials parameters and the geometry of the link. The time 
dependence of the bolometer temperature can be solved 
from the heat equation for a bolometer with a bias point 
resistance of i? = V/I absorbing a time-varying optical 
signal Qopt(i) = Qoc"^' (c.f., Eq. 
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where ST is used to denote the temperature change due 
to the signal power and Qbias describes the incoming heat 
flow due to the detector bias. Equating the steady state 
components of the equation yields Qbias = K{T'^ — Tq), 
from which one can obtain the result for the aver- 
age operating temperature of the bolometer, given by 
T — Qbias/Gth + To where an average thermal conduc- 
tance Gth is defined by Gth = is:(r" - T^)/(T - Tq) = 

gbia./[(Qbias/if + T?)'/" " To]- 

The electrothermal feedback manifests itself through 
the fact that the change in input signal power modifies 
the bias dissipation, an effect described by the last term 
in Eq. H63(l . Taking a closer look at the temperature 
change we obtain 
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where Gth = dQout/dT « nKT"~^ is the dynamic ther- 
mal conductance. Now, considering the electrothermal 
term 
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where a — d\ogR/d\ogT describes the sensitivity of 
the detector resistance to the temperature changes and 
(3{uj) = [R - Zs{uj)]/[R + Zs{uj)] is the effect of the 
bias circuit (with an embedding impedance of Zg) on 
the ETF. Taking into account the thermal cut-off of the 
bolometer, the frequency-dependent loop gain is defined 
as 



C{uj) = Co 
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where Co = QuasCt/ (GthT) and tq = G /Gth is the intrin- 
sic thermal time constant of the bolometer. The electro- 
thermal loop gain describes the effect of varying incident 
optical power to the bias power dissipated in the detec- 
tor. For positive bolometers with a > the loop gain is 
positive for current bias (since Re[/3(tj)] > 0) and nega- 
tive for voltage bias (as Re[/3(ti;) < 0]). For bolometers 
with a negative temperature coefficient of resistance the 
situation is reversed. For metallic bolometers operated 
at room temperature a ^ 1 and the loop gain is typically 
small [Co <1) so that the role of ETF is negligible. On 
the contrary, superconducting detectors with a ^ 100 
and Gth some three orders of magnitude smaller than 
for room temperature devices can have large loop gain 
{Cq > 50) , so that ETF plays a significant role in the de- 
tector characteristics. A major impact of negative ETF 
is that the bolometer time constant is reduced from tq 
to Toff To/il + /5(0)£o]- This reduction in the time 
constant is one of the major benefits of strong negative 
ETF. 

Voltage biasing conditions are typically reached by 
driving a constant current /q through a parallel com- 
bination of the bolometer and a load resistor Z^- As 
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long as <C R, the bolometer is effectively voltage bi- 
ased. The responsivity of a voltage biased bolometer can 
be derived as follows: The current responsivity is de- 
fined as 7^I = dl/dQo- Using Eqs. and (|^ . 
V = IoZsR/{Zs + R), and / = IoZs/{Zs+R), the result 
for the current responsivity becomes 
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(67) 



In the limiting case at w = with (3 — 1 (perfect voltage 
bias) and £o S> 1, 



7^1(0)---^. 



(68) 



In order to evaluate the NEP for a bolometer, a dis- 
cussion on the noise sources in bolometers is merited. In 
general, the NEP, noise spectral density S, and respon- 
sivity are related by NEP = y'Svj/l'R-vjl, where the sub- 
scripts refer to voltage or current noise and responsivity, 
respectively. The noise in bolometers is due to several 
uncorrelated sources. The most important of them is 
the thermal fluctuation noise, mentioned already above. 
Generally, this contribution is given by 



NEPtfn = V^^fkBTfG 
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Here 7 describes the effect of the temperature gradient 
acros s the the rmal link between the sensor and the heat 
sink l|Matheil[l982,) in the diffusive limit (i.e., no ballistic 
heat transport present) as 
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Here we assumed that the thermal conductivities obey 
K oc T^. For resistive bolometers, another important con- 
tribution to the NEP is due to the Johnson noise, given 
by 



NEPj = 
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In addition to these noise sources, the current noise of 
the amplifier, ^i^^^p, adds a contribution NEPamp = 
y/S^^^^/Hi- The total NEP of a bolometer is then 
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where the slightly ambiguous contribution NEPcxcoss en- 
compasses various contributions from additional external 
and internal noise sources. Typical excess external noise 
contributions arise from heat bath temperature fluctua- 
tions and pickup in the cabling to name a few. In addition 
to the external excessive noise, it has in recent years be- 
come clear that there are also internal noise sources that 



are not fully accounted for. For instance, TESs with 
significant internal thermal resistance can no longer be 
treated using the simple lumped element model, and they 
suffer from internal th ermal fluctuation n oise (ITFN), 
adding a contribution (|Hoevers et all l200(]|) 



NEP 



ITFN — 



ikBTR 



Lo 



G 



th 



'0 ' 



(73) 



where Lq is the Lorenz number. 

For some devices, these noise sources are sufficient to 
explain all the observed noise. However, many groups are 
developing X-ray microcalorimeters (see below) which 
are often operated at a small fraction of the normal state 
resistance exhibit noise that increases rapidly as the bias 
point resistance is decreased, and has a significant influ- 
ence on the performance of the detectors. Several possi- 
ble explanations have been put forth, e.g., noise arising 
from the fluctuations in magnetic domains or phase-slip 
lines l)Knoedleil Il983 IWollman et adll997|) . A system- 
atic study of the excess noise in different TES geometries 
has been published recently ( UUom et a/.l . 12004 '). showing 
that there exists a clear correlation between a and the 
observed excess noise, with the magnitude of the excess 
noise scaling roughly as 0.2^/a. A quantitative agree- 
ment with the measured excess noise spectrum has so 
far b een achieved in one experiment ijLuukanen et all 
I2003D . where the TES consisted of an annular (so called 
Corbino) geometry with a superconducting center con- 
tact, and a concentric current return at the outer perime- 
ter of the annular TES. This geometry results in strictly 
radial current fiow, enabling a simple analytical expres- 
sion for the current density in the TES. The 1/r depen- 
dence in the current density together with a small radial 
temperature gradient causes the TES to separate to two 
annular superconducting and normal state regions. 

For any system undergoing a second order phase tran- 
sition, order parameter fluctuations will take place. The 
excess noise arises from the thermally driven fluctua- 
tion of the phase boundary which manifest as resistance 
fluctuations. The volume associated with the order pa- 
rameter fluctuation can be obtained by noting that the 
Ginzburg-Landau free energy SF associated with the 
fluctuation is SF ~ fceTc- As the TES is biased towards 
smaller R, the relative contribution of the order parame- 
ter fluctuations becomes larger until it fully dominates 
over the other noise contributions. For the Corbino- 
geometry TES, the contribution due to the fluctuation 
superconductivity noise (FSN) to the NEP is 
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where F m 10~^ K/vHz is a constant dependent on the 
TES parameters. 

So far, this noise model has not successfully been ap- 
plied to TESs in the more conventional square geome- 
try, mainly due to the fact that the current distribution 
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varies with bias point and is not easily calculable. A so- 
lution could be obtained by solving the full 2D Ginzburg- 
Landau equations for a square geometry. 



1. Hot electron bolometers 

In principle, modern bolometers can be artificially di- 
vided into two major sub-classes depending on where 
the dominant thermal bottleneck lies. The so-called hot- 
electron bolometers (HEBs) utilize the decoupling of the 
electron gas in a metal or a semiconductor from the 
phonon heat bath. The earliest HEBs were based on 
InSb, where the weak coupling of the electrons to the 
lattice at temperatures around 4 K allows the electrons 
to be heated to a temperature significantly above that 
of the lattice even in a bulk sample. Mobility in InSb 
is limited by ionized impurity scattering which results 
in decreasing resistivity with an increasing elect ric field . 
This is the basis o f the detection mechanism l|BrownL 
119841: iRollinL Il96l[l . Often HEBs have sufficient speed 
for mixing, and in fact HEB mixers are a current topic 
of considerable interest. Most HEBs mixers operate at 4 
K, and in order to maintain our focus on phenomena and 
devices at temperatures below 1 K, we unfortunately will 
not discuss them within this review . 

HEB direct detecto rs l|Ali et ali l2003t iKarasik et aZI . 
1200,'^ iRichardj . Il994l) have been a subject of consider- 
able interest, but in general a full optical demonstration 
remains to be carried out. The attractive features of 
such devices include very short time constant (well be- 
low 1 jjLS), potentially very good NEP performance (be- 
low 10~^^ W/\/Hz when operated at or below 0.3 K), 
and simple construction that does not require surface mi- 
cromachining steps. The architecture is essentially very 
similar to that of the HEB mixers: a small supercon- 
ducting TES film coupled to the feed of a lithographic 
antenna. The application of the SQUID readout scheme 
utilized in typical hot phonon microbolometers and mi- 
crocalorimeters might not be as straightforward as one 
could expect as the introduction of the SQUID input coil 
inductance to the voltage biasing circuit can make the 
system unstable due to the interaction of the poles of the 
electrical circuit and the thermal circuit. The approxi- 
mate criterion for the stability of a voltage biased TES 
is that the effective time constant of the TES should be 
about one order of magnitude lo nger than the elec trical 
time constant of the bias circuit ijirwin et 

A hot-electron bolometer that has been demonstrated 
is based on electron thermometr y with NIS junctions 
l|Nahum and Martini j . fl993lll995j) . Here the incident op- 
tical power elevates the electron temperature in a small 
normal metal island weakly thermally coupled to the lat- 
tice phonons, and the change in the electron tempera- 
ture can be sensed as a change in the tunneling cur- 
rent (see Subs. UTT.A.ll) . Noise equivalent powers below 
10-19 W/VHz have been predicted ijKuzminl 1200^ but 
remain to be experimentally verified. An additional at- 



tractive feature of the SINIS bolometer is that a DC bias 
on the device can be used to refrigerate the electrons 
to a temperature below that of the bath temperature 
(see Subs. lV.C.l|l . Another benefit over TESs is the fact 
that the SINIS bolometer saturates much more gently 
compared to the TESs, which basically have no response 
at all once the device is overheated above Tc. The self- 
cooling property of the SINIS bolometer can also be used 
to compensate for excessive background loading, thus ef- 
fectively giving it a larger dynamic range. The main ob- 
stacle towards constructing large arrays of (SI)NIS based 
HEBs is that their impedance (typically 1 kJl - 100 kfi) 
is hard to match to the ex i sting c ryogenic SQUID mul- 
tiplexers llChervenak efall. 119991: Ide Korte et aH 120031: 
iLanting et al', 2005^ iReintsema et all 120031 : lYoon et alL 
|2001). In principle, one could apply superconducting 
transformers to match the SQUID noise, but transform- 
ers with sufficient impedance transformation range are 
quite large, which makes this approach unpractical. A 
novel readout method that lends itself for array readouts 
is a microwave reflectometric measurement, in which the 
SINIS b olometer is connected in ser ies with a tuning in- 
ductor l|Schmidt et a/.U2003ll2004b(l . The LC resonance 
frequency of the inductor and the stray capacitance of 
the junctions is tuned to fall within the bandwidth of 
the cryogenic microwave amplifer (400-600 MHz), facili- 
tating good impedance match. The dynamic resistance 
of the device is highly sensitive to the electron temper- 
ature, and thus temperature changes cause modulation 
of Q of the resonance circuit. This modulation is sensed 
by sending a small RF signal to the resonant circuit, and 
measuring the refiected power. The electrical NEP in- 
ferred from noise measurements was in these experiments 
1.6 X 10-1^ W/Vb^. 



2. Hot phonon bolometers 

The second major class of bolometers are hot phonon 
bolometers (HPBs). They rely on a geometrical de- 
sign of the heat link Ggoom so that the thermal bot- 
tleneck lies between two phonon populations. This 
approach is the most common, and allows for opera- 
tion at temperatures up to and beyond room temper- 
ature. The earliest and most widely used of the con- 
temp orary HPBs are the s o-called spider-web bolome- 
ters ijMauskopf et all Il997|) . operated at 300 mK and 
below, where a free-standing Si3N4 mesh is used to sup- 
port a thermal sensing element. Narrow Si3N4 legs pro- 
vide the thermal isolation for the mesh. Before the intro- 
duction of TESs and SQUIDs the thermometer of choice 
was a small crystal of neutron-transmutated (NTD) Ge 
due to its relatively high temperature coefficient of resis- 
tance ((ilog/i/(ilogT|T=o.3K ~ —6) and large resistance 
(~ 25 Mfi) that allowed for good noise matching with 
FET preamplifiers. On the other hand this made the 
devices very microphonic. Efficient optical coupling was 
possible since radiation at wavelengths smaller than the 
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mesh period are absorbed to a resistive film deposited on 
the Si3N4 mesh. An example of a spider-web bolome- 
ter using a NTD Ge is shown in Fig. El Later ver- 
sions of the spider-web bolometers have adopted the use 
of TES as thermom eters, coupled to a SQUID readout 
l|Cildemeister et allFM. 




FIG. 17 (Color in online edition): A micrograph of a "spider- 
web" bolometer. A NTD Ge thermistor is located at the 
centre of the web. Image courtesy of NASA/JPL-Caltech. 



Alternatively to the spider-web absorber, a A/4 reso- 
nant cavity can be used to maximize the optical efficiency 
over a limited bandwidth. The SCUBA -2 instrument 
ijDuncan et a^J . l2003HHolland et al\ . \200^ on the James 
Clerk Maxwell telescope is an ambitious overtaking with 
a pixel count of over 12 000. The bolometers consist 
of thermally isolated, (1 mm^)x 60 nm silicon 'bricks' 
with the front surface of the Si degenerately doped with 
phosphorous to 377 per square. The resonant cavity 
is formed by the Si in between the doped layer, and a 
Mo-Cu TES deposited to the back side of the pixel. 

Even though the superconducting TES thermometers 
are becoming increasingly popular, bolometers utilizing 
lithographic semiconducting thermistors still yield im- 
pressive performance and are more forgiving in terms of 
the saturation power. In the so-called pop-up bolome- 
ters, a lithographic, doped Si is used as the thermis- 
tor, while an ingenious torsional bending method of the 
Si3N4 legs is used to bend the legs and wiring layers per- 
pendicular to the absorbers. This architecture enables 
the construction of ar rays with a very high filling factor 
l|Voellmer et a?J . l2003|) . Another promising lithographic 
resistive thermometer technology is based on thin films of 
NbSi (,,Camus et al., ,2000,) ■ The high resistivity of these 
films allows for impedance levels that provide a good 
noise matching to room temperature JFET amplifiers. 

Instead of the surface absorbing approach, another 
method is ba s ed on the use of a lithog raphic antenna 
l)Hwang et al 1 1X9791: iNeikirk et oii l 1983), ^-iiti terminat- 
ing the induced currents to a thermally isolated bolome- 
ter, an approach often used in the HEB mixers discussed 
above. The attractive feature of this method is that 



the thermally sensing volume can be made much smaller 
compared to the case of the surface-absorbing bolome- 
ters, making these devices much faster. The low heat 
capacity often allows for lower NEP, as in many cases 
the NEPtfn is limited by the maximum time constant 
Tmax — C/Gth,min allowcd by the application. In this case 
NEPtfn — \J ^ik-oT'^C / Tn^a.yi- Moreover, the small size 
of the thermally sensitive volume makes it far less sen- 
sitive to out of band stray light, relaxing filtering and 
baffling requirements of the incoming radiation. An- 
tenna coupling also lends itself to the construction of 
integrated, on-chip filters for defining the bands for an 
array of bolometer s ijHunt et g^J . 120031: iMees et all.ll99lt 
iMvers et a/.l . l2005l) . Arrays of antenna-coupled bolome- 
ters can also utilize the inherent polarization selectivity 
of the antennas. An example of an antenna-coupled TES 
bolometer that incorporates on-chip transmission line 
impedance transformers and band-pass filters is shown 
in Fig. ^1 No NEPs have yet been measured on this 
device, but the expected NEP is 10^^^ W/VHz for a 
device with Tc = 450 mK. 




FIG. 18 (Color in online edition): Micrograph of an antenna- 
coupled TES bolometer. A dual-polarized double-slot an- 
tenna (a) is coupled via microstrip transmission lines (b) with 
band-pass and low-pass filters (c) to two Al/Ti bilayer TES 
thermometers located on suspended Si3N4 membranes (d). 
The inset shows in detail the termination of the microstrip 
line (e) to a resistor (f), and the Al/Ti TES (g). Figure cour- 
tesy A. T. Lee, UC Berkeley. 



While not published yet, the best performance of an 
antenna-coupled HPB utilizing a similar design has been 
obtained by the JPL-Caltech bolometer group, with an 
electrical NEP=5 • lO'^^ W/Vli at 230 mK l|KenvonL 
1200,'tD . The architecture is shown in Fig. 

In the simplest case, the bolometer is simply a strip of 
metal placed to the feed of the lithographic antenna. Al- 
though this approach is by far the simplest, it introduces 
limitations for the performance of a HPB: Unlike in the 
case of HEB mixers, maximizing the thermal isolation of 
the bolometer is always desirable. However, the bolome- 
ter resistance should be matched to the impedance of the 
antenna (typically lOO for broadband lithographic 
antennas on Si substrates) , and thus the thermal conduc- 
tance to the bath through the antenna, Gant j is fixed by 
the Wiedemann- Franz law, Gant = ioT/Re(Z'a), where 
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riety of applications, from the detection of weakly in- 



FIG. 19 SEM micrograph of an isolated Si3N4 platforms with 
Mo/Au TESs. As an extreme example of thermal isolation, 
the meandering Si3N4 legs have an aspect ratio of ~ 2800:1, 
yielding a thermal conductance Gth ~ 100 fW /K at a bath 
temperature of 210 mK. Figure courtesy M.E. Kenyon, JPL. 



Za is the antenna impedance. A parallel heat loss path 
is to the substrate below the bolometer film, with con- 
ductance Gsub oc LW where L and W are the length 
and width of the bolometer film, respectively. For sub- 
strate mounted antenna-coupled HPBs it is thus bene- 
ficial to minimize the size of the bolometer. This re- 
quirement can be relaxed if the bolometer strip is re- 
leased from the underlying substrate to form an air- 
bridge (JNeikirk and Rutlcdgc, 1984). For air-bridged de- 
vices at temperatures higher than 4 K, NEPtfn oc T'^/^. 
Recently, an air-bridge bolometer operating at 4.2 K 
was demonstrated, showing potential for background- 
limi ted performance when obs erving 300 K blackbod - 
ies (iLuukanen et "oil. 12005: .Luukanen and PekolaL l2003|) . 
The potential applications for these devices include pas- 
sive detection of concealed weapons under clothing, re- 
mote trace detection, and terrestrial submillimetre-wave 
imaging. 



C. Calorimeters: Pulsed excitation 

In the limit opposite to the bolometric detection, i.e. 
when the mean time between the quanta of energy ar- 
riving at the detector exceeds the device relaxation time, 
thermal detectors are known as calorimeters. While the 
topic of calorimetry also encompasses heat capacity mea- 
surements especially in mesoscopic samples, the following 
discussion concentrates on the detection of radiation only 
in order to keep the scope of our review in reasonable 
limits. For those interested in microcalorimetry in the 
sense of heat ca pacity measureni e nts, w e direct the reader 
to references ([Bourgeois a/.l. l2 005t iDe nlinger a^J . 
119941: iFominava et a^J Il997t iLindell et all |2000() and 
ijMarnieros et al\ . \l999(\ . 

Cryogenic calorimeters are used today in a large va- 
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calorimeters are usually operated at temperatures below 
1 K. The theory and operation of calorimetric thermal 
detectors is very similar to bolometers, and generally the 
theoretical treatment above is valid. However, the opti- 
mization of calorimetric detectors can be quite different. 
The quantum of energy E deposited by either a photon, 
charged particle, WIMP etc. can be determined from the 
temperature rise AT = E/C, where C is the heat capac- 
ity of the calorimeter. This temperature rise then decays 
exponentially with a time constant TcS to its equilibrium 
value. Thermometry in the calorimeters is most often 
done either using semiconducting or TES thermometers, 
as is the case with bolometers. In addition, thermometry 
based on the change of magnetization of a paramagnetic 
sensor is appealing due to its non-dis sipative nature, and 
has yielded some proin ising results ijFleischmann et aZI . 
l2003tlSchonefeld et allEom . 

The figure of merit for a calorimeter is the energy res- 
olution, AE, of full- wid th at half-max i mum (FWHM), 
related to NEP through l|Moselev et a.ll\l984) 
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For a 'classic' calorimeter with a white noise spectrum the 
energy resolution is in terms of the operating temperature 
and the heat capacity given by 



AE = 2%/21n2v/fcBr2C, 



(76) 



indicating that AE scales oc T^/^ or oc T^/^ depending 
whether the heat capacity of the sensor is dominated by 
the lattice or the electronic system, respectively. 

In most of the devices under development today the 
lattice temperature is being measured, as sufficient cross 
section to the incoming energy often requires a rather 
large volume of a device. An exception to the norm are 
TES calorimeters optimized for optical single photon de- 
tecti on for applications in secure quantum key distribu- 
tion ijNam et a/.ll2004|l . shown in Fig. EOI In these de- 
vices, the thermal isolation is via electron-phonon decou- 
pling. A trade-off is made between energy resolution and 
the speed of the detectors. Energy resolution of the de- 
tectors is sufficient to determine the photon-number state 
of the incoming photons, while maintaining a speed that 
is adequate for fast information transfer. 

As with bolometers, the TESs currently outper- 
form the competition in terms of the sensitivity. The 
best reported energy resolution for any energy dis- 
persive detector was recently obtained with a Mo-Cu 
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FIG. 20 (Color in online edition): The energy spectrum of 
a pulsed 1550 nm laser, measured with an optical TES mi- 
crocalorimeter. The peaks in the plot correspond to the 
photon-number state of the incoming pulses. The inset shows 
a micrograph of the devices. Figure courtesy of A.J. Miller, 
NIST. 



calorimeter (see Fig. I21|l . yielding an energy resolu- 
tion of 2.38ib 0.1 1 eV a t a photon energy of 5.89 keV 
l|Irwin and HiltoiA l2005|) . The driving application in 
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FIG. 21 An energy spectrum of the ^^Mn Ka complex, ob- 
tained with a Mo/Cu TES microcalorimeter. The width of 
the measured characteristic x-ray lines is a convolution of the 
intrinsic x-ray linewidth and the gaussian detector contribu- 
tion with a FWHM of 2.38 ±0.11 eV. The TES is fabricated 
on a free-standing Si3N4 membrane. The inset shows a SEM 
micrograph of the TES. The normal metal bars extending 
partially across the TES have been experimenta lly verified to 
impr ove the energy resolution of the detector llUllom et all 
\2004) . Figure courtesy of J. Ullom, NIST. 



the development of X-ray microcalorimeters are two ma- 
jor X-ray astrophysics missions planned by the Euro- 
pean Space Agency (ESA) and the U.S. National Astro- 
nautics and Space Administration (NASA). Both mis- 
sions, AT-ray Evolving f/niverse S'pectroscopy [XEUS) 



l)The European Space Agenc y. '2005") mission and the 
Constellation-X mission ||nASA . .2005^) will employ X- 
ray microcalorimeters as their primary instrument. The 
energy resolution of the state of the art microcalorimeters 
have already reached the requirements stated in the sci- 
ence goals for these missions. Thus, the current primary 
focus of the technical research is on the development of 
large arrays of microcalorimeters, and more importantly, 
SQUID multiplexing read-outs for the detector arrays. A 
prototype 5x5 array of Ti/Au TES microcalorimeters is 
shown in Fig. [221 
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FIG. 22 (Color in online edition): A prototype Ti/Au mi- 
crocalorimeter array. The dark regions within the Si3N4 mem- 
brane are holes in the membrane. Figure courtesy of SRON- 
MESA. 



D. Future directions 

While cryogenic thermal detectors are the most sen- 
sitive radiation detectors around, further significant per- 
formance increases are envisioned in the future. The push 
is from single pixels to large staring focal plane arrays 
of detectors for numerous astrophysics applications. For 
this reason, novel thermal detector concepts must be able 
to be integrated into large arrays. 

The performance of hot electron bolometer mixers has 
been improving rapidly over the past years. Noise tem- 
peratures approaching 10 x hv/(2kB) are being reported, 
and the emphasis is in pushing the operating frequency 
deeper into the THz region. In principle, frequencies in 
the infrared range are not out of the question. Litho- 
graphic antenn as have demonst r ated good performance 
up to 30 THz (|Grossman et all Il99ll) . while it is clear 
that the fabrication becomes increasingly challenging as 
the required feature size is reduced. 

In the single pixel direct detector development, the 
quest for ever better sensitivity is still ongoing. For 
bolometers, the improvement in NEP directly translates 
to the capability of observing over less pre-detection 
bandwidth without sacrificing signal to noise ratio. This 
would enable the construction of arrays capable of yield- 
ing spectroscopic information. It is here where the func- 
tions of HEB mixers, bolometers and calorimeters conjoin 
- in the capability of performing single photon spectrom- 
etry at far-infrared wavelengths. 
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V. ELECTRONIC REFRIGERATION 

Thermoelectric effects and, in particular, tficrmoelec- 
tric cooling have been discovered more than 170 years 
ago l)Peltieilll834j) . During the last 40 years considerable 
progress has been made in developing practical thermo- 
electric refrigerators for industrial and scientific appli- 
cations (Nolas et at, 2001; Rowc and Bhandari, 198^. 
The temperature range of interest has been, however, 
far above cryogenic temperatures. Yet, during the last 
decade, solid state refrigerators for low temperature ap- 
plications and, in particular, operating in the sub-kelvin 
temperature range have been intensively investigated. 
The motivations of this activity stem from the success- 
ful development and implementation of ultrasensitive ra- 
diation senso rs and quantum circ uits which require on- 
chip cooling jPekola et all l2004hD for proper operation 
at cryogenic temperatures. Solid state refrigerators have 
typically lower efficiency as compared to more traditional 
systems (e.g., Joule-Thomson or Stirling gas-based re- 
frigerators). By contrast, they are more reliable, cheaper 
and, what is more relevant, they can be easily scaled 
down to mesoscopic scale. All this gives a unique op- 
portunity to combine on-chip refrigeration with different 
micro- and nanodevices. 

The aim of this part of the review is to describe the ex- 
isting solid state electron refrigerators operating at cryo- 
genic temperatures (in particular below liquid helium 
temperature) , and to give an overview of some novel ideas 
and refrigeration schemes. 



A. General principles 

The physical principle at the basis of thermoelectric 
cooling is that energy transfer is associated with quasi- 
particle electric current, as shown in Sec. II. Under suit- 
able conditions thermal currents can be exploited for heat 
pumping, and if heat is transferred from a cold region to 
a hot region of the system, the device acts as a refrig- 
erator. The term refrigeration is associated throughout 
this Review to a process of lowering the temperature of 
a system with respect to the ambient temperature, while 
cooling, in general, means just heat removal from the 
system. It is noteworthy to mention that the maximum 
cooling power of a refrigerator is achieved at a vanishing 
temperature gradient, while the maximum temperature 
difference is achieved at zero cooling power. The effi- 
ciency of a refrigerator is normally characterized by its 
coefficient of performance (rj), i.e., the ratio between the 
refrigerator cooling power (Qcooier) and the total input 
power (Ptotai)- 
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Irreversible processes (e.g., thermal conductivity and 
Joule heating) degrade the efficiency of refrigerators, and 



FIG. 23 Basic Peltier thermoelement. 



are essential elements that need to be carefully evaluated 
for the optimization of any device. 

The basic principles of the design are in general sim- 
ilar for different types of refrigerators. The increase of 
temperature gradient can be achieved by realizing a mul- 
tistage refrigerator. In this case, the stage operating at 
a higher temperature should be designed for larger cool- 
ing power in order to efficiently extract the heat released 
from the lower-temperature stage ("pyramid design"). 
The enhancement of the refrigerator cooling power can 
be achieved by connecting several refrigerators in par- 
allel. The parallel design is more effective both for an 
efficient heat evacuation from the hotter regions of the 
device and for the application of higher electric currents 
to the refrigerator. It also allows more freedom to prop- 
erly bias all cooling elements. 

The temperature dependence of the electric and ther- 
mal properties of the active parts in the refrigerator may 
limit their exploitation at low temperatures. The reduc- 
tion of thermal conductivity by lowering the temperature 
has both positive (better thermal insulation between cold 
and hot regions) and negative (difficulty in removing heat 
from the system) effects. 

At cryogenic temperatures different types of supercon- 
ductors can be efficiently exploited. They can be used 
both as passive and active elements: in the former case, 
owing to their low thermal conductivity and zero elec- 
tric resistance (e.g., as one of the two arms in Peltier 
refrigerators), in the latter as materials with an energy 
gap in the density of states for energy-dependent electron 
tunneling (e.g., in NIS coolers). 

Currently, the development of refrigerating techniques 
follows two main directions: search of new materials with 
improved characteristics suitable for existing refrigera- 
tion schemes, and development of new refrigeration meth- 
ods and principles. 



B. Peltier refrigerators 

Thermoelectric (Peltier) refrigerati on is widely used 
for cooling di f ferent electronic devices llNolas et a^J . l200il 
iPhelan et all. l2002t iRowe and BhandariL Il983|) . Nowa- 
days Peltier refrigerators providing temperature reduc- 
tion down to 100. ..200 K and cooling power up to 100 
W are available. Peltier cooling (or heating) occurs 



33 



when an electric current is driven through the junction of 
two different materials. The heat released or absorbed, 
Q Peltier, depending on the direction of the electric cur- 
rent at the junction, is proportional to the electric cur- 
rent (/) driven through the circuit, Qpeiuer = ^abI, 
where IV ab — ciabT, and IIab and uab are the Peltier 
and Seebeck coefficients of the contact, respectively (see 
also Eq. H39|) '). In order to obtain enhanced values of 
the Peltier coefficient, conventional Peltier refrigerators 
generally consist of p- and n-type semiconductors with 
opposite sign of 11 coefficients (see Fig. 1^ . The ef- 
ficiency of a Peltier refrigerator is not determined only 
by the coefficient Hab — — lis, but also by thermal 
conductivities (k) of both materials across the contact. 
Furthermore their electric resistances (p) are responsible 
of Joule heating affecting the coefficient of performance. 
The maximum temperature difference (A Tpn.-„) achiev- 
able with a Peltier refrigerator is given by l|Nolas el alV 

I200I AT^ax = ZTU2, where Z = ^ is a figure 
of merit of the refrigerator, and Tcoid is the tempera- 
ture of the cold junction (see Fig. [^ . More often, 
however, the refrigerator efficiency is characterized by 
the dimensionless figure of merit ZT. We recall that 
ZT oc (■^^)^. Most of the materials used in thermo- 
electric applications have ZT ^ 1 at room temperature 
l)DiSalvolll999tlMin and Rowell200(]() . In general, at low 
temperatures only one single thermoelectric material is 
needed, because a superconducto r can be used as one of 
the two arms of t he refrigerator ijGoldsmid et all Il988t 
iNolas etallEim^. 

In spite of the drastic reduction of cooling efficiency 
at low temperature, there is, however, some develop- 
ment of new materials and devices suitable for opera- 
tion at cryogenic temperatures. Recently, Peltier cool- 
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0.2 K was demonstrated below 10 



mg vviLiii 

K (H arutvunvan et al\ . l2003|) (see Fig. Crystals 
of CeBe were used to exploit the strong thermoelectric 
coefficients arising from the Kondo effect. The dimen- 
sionless figure of merit of this material is 0.1. ..0.25 in the 
temperature range between 4 K and 10 K. The authors 
claim that a proper optimization of the refrigerator de- 
sign would allow more than 10% temperature reduction 
below 4K with a single-stage refrigerator. 

At millikelvin temperatures lattice specific heat and 
thermal conductivity decrease drastically, and thermo- 
electric refrigeration might become feasible. Following 
this idea and taking into account the Wiedemann-Franz 
relation, values of ZT as high as 20 at temperatures be- 
low 10 mK and the possibilit y to achieve the rmoelectric 
refrigeration were predicted ijGoldsmid and Grav. .19791 
iNolas et QiJ.l2001^). 

iKapitulniS l)l992 l) proposed to exploit a metal close to 
its metal-insulator transition for the implementation of 
a thermoelectric refrigerator operating below liquid-He 
temperatures. The basic concept behind this proposal 
is that near the metal-insulator transition, the transport 
coefficients acquire anomalous power laws such that their 
relevant combination describing the figure of merit for 
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FIG. 24 Thermoelectric refrigeration at cryogenic tem- 
eratures using cerium hexaboride. Adapted from 
iHarutvunvan et aJlH)c!^. 
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FIG. 25 Schematic diagram of the thermoelectromechanical 
cooler, time sequences of the pulsed current applied to the 
device, and the two modes of canti lever contact: synchronized 
and optimized iMiner et aZ.llT999|) . 



efficient cooling also becomes large. 

Further improvement of the efficiency of conven- 
tional thermoelectric coolers could be, in principle, 
achieved through a thermoelectromechanical refrigerator 
l)Miner et all Il999t iThonhauser et all l2004j) . In such a 
device, a periodic variation of the electric current through 
a Peltier element combined with a synchronized mechan- 
ical thermal switch would allow to improve the overall 
cooling performance (see Fig. I25|l . The enhancement of 
refrigeration is due to the spatial separation of the Peltier 
cooling and Joule heating. 

The analysis of thermoelectric devices is usually based 
on the parameters typical of bulk materials. Significant 
progress in low temperature Peltier refri geration might 
be a c hieved by using ex o tic ni aterials l)Goltsev et all 
\200A iRontani and ShamL 12 000*) and low-dimensional 
struc tures (Balan din and Laz arcnkova. 2003; DiSalvol 
119991: iHicks et 

^,119917 

such as composite thin films, 
modulation-doped heterostructures, quantum wires, nan- 
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otubes, quantum dots, etc. These systems offer, in gen- 
eral, more degrees of freedom to optimize tfiose quantities 
that affect the efficiency of thermoelectric refrigerators. 



C. Superconducting electron refrigerators 

1. (SI)NIS structures 

Although heat transport in superconducting mi- 
crostructures origina lly dates back more than 40 years 
ago in SIS junctions ("Chi and Clarld Il979l: iGravl Il978t 
iMclt on et ai, 1981; Parmcntei'. 1961|) , it is instructive to 
start our description of this topic from NIS tunnel junc- 
tion structures. FigureEHa) shows the calculated Q for a 
NIS tunnel junction (see Sec. II. G. 2) versus bias voltage 
at different temperatures (T = T^^n ~ Te.s)- When Q 
is positive, it implies heat removal from the N electrode, 
i.e., hot excitations are transferred to the superconduc- 
tor. For each temperature there is an optimal voltage 
that maximizes Q and, by decreasing the temperature, 
the heat current results to be peaked around V ~ A/e. 
Figure I26r bl displays the heat current versus tempera- 
ture calculated at each optimal bias voltage. The quan- 
tity Q{T) is maximized at T w 0.25 A/ks = Topt (as 
indicated by the arrow in the figure) where it reaches 
Q — 6 X 10~^A^/e^i?T, decreasing both at lower and 
higher temperatures. Also shown in the figure is Q{T) ob- 
tained assuming a temperature-independent energy gap. 
Such a comparison shows that this latter assumption is 
fully justified for T < 0.2 A//cb. In the low temperature 
limit (Te^jv < 7e,s ^ A/fce) it is possible to give an ap- 
proximate expression (Anghcl and Pckola, 2001) for the 
optimal bias voltage {Vopt), Vopt ~ (A — 0.66fcsre.Ar)/e, 



as well as for the maximum cooling power at Vc 
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Equation H78|l is useful for getting quantitative esti- 
mates on the performance of realistic coolers. In the 
same temperature limit and for V — Vopt the cur- 
rent through the NIS junction can be approximated as 

/ « 0.48^ \J^^^^- The NIS junction coefficient of 

performance -q is given by ri{V) = Q{V) /[I{V)V]. For 
V ~ A/e and in the low temperature limit, -q thus ob- 
tains the approximate value 



riopt 



0.7- 



T, 



e,N 



T, 



(79) 



where we assumed A = 1.764 /csTc, and Tc is the critical 
temperature of the superconductor. Equation (j79|l shows 
that the efficiency of a NIS junction is around or below 
20% at the typical operation temperatures. The full be- 
havior of rj versus temperature calculated at each optimal 
bias voltage is displayed in Fig. I26r c') . The simple results 
presented above point out how the optimized operation of 



a superconducting tunnel junction as a building block of 
microrefrigerators stems from a delicate balance among 
several factors such as the contact resistance, the opera- 
tion temperature, the superconducting gap A as well as 
the bias voltage across the junction. 

The first observation of h eat extraction from a normal 
metal dates back to 1994 (|Nahum et oil Il994l . where 
cooling of conduction electrons in Cu below the lattice 
temperature was demonstrated using an AI/AI2O3/CU 
tunnel junction. A significant improvement was made 
two years later, still in the AI/AI2O3/CU system, by 
iLeivo et al\ ()l996j) . which recognized that using two NIS 
junctions in series and arranged in a symmetric configu- 
ration (i.e., in a SINIS fashion) leads to a much stronger 
cooling effect. This fact can be understood keeping in 
mind that Q is a symmetric function of V so that, at 
fixed voltage across the structure, quasi-electrons above 
A are extracted from the N region through one junc- 
tion, while at the same time quasi-holes are filled in 
N below —A from the other junction (see Fig. EHI^d)). 
In this configuration, a reduction of the electron tem- 
perature from 300 mK to about 100 mK was obtained. 
Later on, several other experimental evidences of elec- 
tron cooling in SINIS metallic structures were reported 
(jArutvunov et a/.. .2000; .Clark et aLt . .2004: .Fisher g/l 
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FIG. 26 Characteristics of (SI)NIS refrigerators, (a) Calcu- 
lated cooling power Q of a NIS junction vs bias voltage V 
for several temperatures T = Te,N = Tc.s. Also shown is the 
behavior of a NIN junction, (b) Q calculated at the optimal 
bias voltage as a function of temperature, assuming both a 
temperature-independent energy gap (dash-dotted blue line) 
and the real BCS dependence (black line). Topt indicates the 
temperature value that maximizes Q. (c) Coefficient of per- 
formance Tj calculated at the optimal bias voltage versus tem- 
perature, (d) Scheme of the energy band diagram of a voltage 
biased SINIS junction. The electric current flows into the nor- 
mal region through one tunnel junction and out through the 
other, while the heat current Q flows out of the N electrode 
through both tunnel junctions. 
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normal electrode l|Rowell and TsuiL 11976') . Figure HT^c 
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FIG. 27 SINIS refrigeration, (a) Scanning electron micro- 
graph of a typical SINIS microrefrigerator. The structure 
was fabricated by standard electron beam lithography com- 
bined with Al and Cu UHV e-beam evaporation. The probe 
junction resistance usually satisfies the condition Rp S> Rt 
in order to prevent self-cooling, (b) Sketch of a typical mea- 
surement setup for electron cooling, (c) Electron tempera- 
ture in the N region versus Vrefr measured at different bath 
temperatures. Dots represent experimental data, while solid 
lines are theoretical fits. The inset shows the cooling power 
(pW) against temperature obtained from the fits (dots), (d) 
Contour plot of the theoretical electron temperature T^^n vs 
bias voltage V and Te,s for a SINIS cooler, assuming thermal 
load due to the electron-phonon interaction (se e text), (a) 
is ad apted from IPekola et al. ■2004al : (c) from iLeivo et all 
Il996r). 



1995", '1999^ 'Leivo et al', '1997'; 'Leoni et al}, 199?, '2003 
Luukanen et a/...,.2000;.Pekola et aL..2000a.,2004a.. 20001 



Tarasov et a/.l . l2003HVvstavkin et all\l999^ . In these ex- 
periments NIS junctions are used to alter the electron 
temperature in the N region as well as to measure it. In 
order to measure the temperature, the N region is nor- 
mally connected to additional NIS contacts (i.e., "probe" 
junctions) acting as thermometers (previously calibrated 
by varying the bath temperature of the cryostat), and 
operating along the lines described in Sec. III.A.l. More- 
over, the differential conductance of the probe junctions 
gives also detailed information about the ac tual quasipar- 
ticle d istributi on fun ction in the N region ijPekola et all 
l2004ai: Pothie r et ad [l997b). 

Figure I27r a) shows the SEM micrograph of a typi- 
cal AI/AI2O3/CU SINIS refrigerator including the super- 
conducting probe junctions. The schematic of a com- 
monly used experimental setup for electron refrigeration 
and temperature measurement is shown in Fig. I27f b). 
The voltage bias K-e/r across the SINIS structure al- 
lows to change the electron temperature in the N re- 
gion; at the same time, a measure of the voltage drop 
across the two probe junctions {Vth) at a constant bias 
current (/q) yields the electron temperature T^^n in the 



illustrates the experimental data of iLeivo et al. 
of the measured electron temperature T = Te^N versus 
Vrefr, taken at different bath temperatures (i.e., those 
at Vrefr = 0)- As Can be readily seen, the electron 
temperature rapidly decreases by increasing the voltage 
bias across the SINIS structure, reaching the lowest value 
around Vrefr ~ 2A ~ 360 /icV (note that now two junc- 
tions in series are involved in the process). Then, by 
further increasing the bias, the electron temperature rises 
due to the heat flux induced into the metal. Furthermore, 
the minimum electron temperature strongly depends on 
the starting bath temperature. The inset shows the ex- 
tracted maximum cooling power (dots) that obtains val- 
ues as high as 1.5 pW at T = 300 mK, corresponding 
to about 2 pW//im^ (in the present experiment, sub- 
micron NIS junctions were exploited, with barrier resis- 
tances Rt — 1 kfJ). 

The general operation of a SINIS microrefrigerator like 
that shown in Fig. I27r a) can be understood by recalling 
that the final electron temperature in the N region stems 
from the balance among several factors that tend to drive 
power into the electron system (i.e., power losses), as dis- 
cussed in detail in Sec. II. Most of metallic SINIS refrig- 
erators fabricated so far operate in a regime where strong 
inelastic electron-electron interaction tends to drive the 
system into quasiequilibrium, where the electron system 
can be described by a Fermi function at a temperature 
Tgjv that differs, in general, from that of the lattice 
(Tp/i). At the cryogenic temperatures of interest (i.e., 
typically below 1 K), the dominant contribution comes 
from electron-phonon scattering that exchanges energy 
between electrons and the lattice phonons. The refriger- 
ator cooling power {Qrefr) can be defined as the max- 
imum power load the (SI)NIS device can sustain while 
keeping the N region at temperature Te^N, and can be 
generally expressed as 

Qrefr = '^Q{V/2; Te^N, Te,s) — Qe-ph{Te^N, Tph), (80) 

where the factor 2 takes into account the presence of two 
identical NIS junctions, and it is often assumed that in 
the superconductors Te^s = Tph- The minimum tem- 
perature Te^N reached by the N metal thus fulfills the 
condition Qrefr = 0. An example of Te^N versus V be- 
havior along the lines of Eq. 180() is represented by the 
solid curves plotted in Fig. I27r c). Furthermore, any ad- 
ditional term that adds thermal load into the electron 
system can be properly taken into account by includ- 
ing its specific contribution to the right-hand side of Eq. 

(Fisher et al, 1999; UUom and Fisher, 2000). Fig- 
ure 123 (d) shows an example of the calculated electron 
temperature Te^N vs voltage and Te,s as obtained from 
the solution of Eq. (|5n|) . Qrefr — 0. Here we consid- 
ered a typical AI/AI2O3/CU SINIS cooler with volume 
V = Q.lb^lTsr•, i?T = 1 kfl, E = 2 nW/^im^K^ (Cu), and 
Te,s — Tph- As it can be readily seen, the Te.N value 
strongly depends on the bias voltage as well as on the 
lattice temperature. Depending on the thermal load due 
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to the electron-phonon interaction and the operating V, 
the SINIS device can thus behave as a cooler or as a 
refrigerator. 

Let us now highhght some of the basic requirements 
of SINIS structures for electron coohng operation (their 
use for lattice coohng is addressed in Sec. IV.C.Gp . In 
particular, the working temperature range of these struc- 
tures depends mainly on the type of the superconductor 
(through its energy gap A and hence its critical temper- 
ature Tc), on the strength of the electron-phonon inter- 
action S, on the junction resistance Rt, and on the N 
region volume V. First of all, a reduction of the active 
volume to be cooled is the most straightforward method 
to increase the efficiency of the refrigeration process, this 
being more relevant at high lattice temperatures (accord- 
ing to the assumed Tp^-dependence of electron-phonon 
interaction). In addition, cooling power maximization in 
the high-temperature regime (let us say 1...4.2 K) would 
require both to lower the electron-phonon coupling con- 
stant and to increase A through a proper choice of the 
superconductor. The first issue can be solved up to some 
extent with a variety of materials with different S (see 
Table nj. As far as the second issue is concerned, there is 
a vast choice of superconducting metals with Tc covering 
the temperature range up to about 20 K (Kittcl, 1996). 

On the other hand, the question of junction resistance 
requires a careful discussion. In general, from Eq. (|78|l . 
Q enhancement is expected from a reduction of Rt- This 
latter issue can be accomplished both by making thinner 
barriers and by increasing the junction area. The first 
issue is not so straightforward, as already discussed in 
Subs. III. F. 21 due to the intrinsic difficulty in fabricating 
high-quality low-i?c barriers, although optimized barri- 
ers are currently under investigation (see Sec. IVI.F.ip . 
Latter option was experimentally addressed by Fisher 
et al. (1999) in Al/Al203/Ag refrigerators, where large 
cooling powers of a few tens of pW were obtained with 
junctions of 20 x 20 surface area. The reduction in 
electron temperature was, however, much inferior to that 
achievable with sub-micron sized junction. The prob- 
lem intrinsic to junctions with large overlap (especially 
at the lowest temperatures) stems from the larger den- 
sity of quasiparticles present in the superconductor, due 
to the fact that quasiparticles require a larger time to 
exit the junction region and escape from the supercon- 
ductor. Therefore, this excess of quasiparticles alters in 
general the refrigerator performance by returning energy 
to the normal electrode, mainly due to back-tunneling 
from the superconductor as well as due to recombina- 
tion processes , where phonons can enter and h e at up 
the N region Oochum et all Il998t iKaolan et aH Il976t 
lUllom and Fi shed . l200(Tfr In addition, inelastic scatter- 
ing with phonons and dynamic impurities can also lead to 
an excess of quasiparticles. These contributions can eas- 
ily overcompensate the junction cooling power, so that 
it is crucial to remove this excess of quasiparticles from 
the superconductor. Toward this end a number of tech- 
niques exist among which we mention the exploitation 
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FIG. 28 Optimization of SINIS refrigerators, (a) SEM micro- 
graph of an AI/AI2O3/CU SINIS microrefrigerator exploiting 
large-area junctions (~ 10 ^m^) with quasiparticle traps, (b) 
SEM image of a part of comb-like SINIS structure with 10-1-10 
junctions for cooling and a SINIS thermometer, (c) Electron 
temperature Te.jv in the N region of an AI/AI2O3CU SINIS 
refrigerator versus Vrefr measured at different bath temper- 
atures. The lowest curve (red dash-dotted line) shows the 
anomalous heating effect observable at the lowest tempera- 
tures and attributed to the presence of quasiparticle states 
within the superconducting gap. (d) Theoretical ultimate 
minimum electron temperature of a SINIS cooler Tmin/Tc at 
V ~ 2 A/e versus F/A in quasiequi librium. (a) is adapte d 
from U Pekola et a P. '20003*): (b) fro mjLuukanen et aZ.ll2000D : 
(c) and (d) from LPckola et aL. ..2004Ejr 



of defect-free and thick S electrodes (that allow quasi- 
particles to escape ballistically from the junction area 
as well as to decrease their density near the barrier) 
and t he exploitation of quasiparticle "traps" jParmenteil 
I1961D . i.e., normal metal films connected to the super- 
conductor in the junction region th rough a tunnel bar- 
rier or in direct m etallic contact l|lrwin et al I Il995bt 
iPekola et all . l2000a|) (see Fig. EHIa)). Such traps act 
as sinks of quasiparticles absorbing almost all the ex- 
cess quasiparticle energy present in the superconductor, 
and have proved to efficiently help heat removal from 
the superconductor leading to a signific ant improvement 
of SINIS device cooling performa nce ijLuukanen et all 
hood iPekola et all l2000al l2004a|) . All these are com- 
monly exploited tricks for the thermalization of the su- 
perconducting electrodes, and in these conditions, it is 
a reasonable approximation to set Tf, 



thermore, the experiments l)Pekola et al 



f l2nonii d( 



Fur- 
demon- 
strated that the trap performance is in general superior 
when it is in direct metallic contact at short distance from 
the cooling junction (typically below 1 /im, although the 
optimal distance mainly depends on the superconductor 
coherence length), nevertheless in small-area junctions 
even a contact through an insulating barrier seems suffi- 
cient for the purpose. The effectiveness of trapping in 
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SI NIS structures wa s theor etically add ressed in detail 
bv IVoutilainen et al\ ll2005() and .Golubev and Vasenkd 
il2nn2D . Another possibility to achieye high cooling by 
maximizing the ratio of the junction area to the size of 
the region to be cooled is to use several small-area junc- 
tions (with size in the sub-micron range) connected in 
parallel to the N electrode (to limit the drawbacks typical 
of large jun ctions) (Arutyuno v et al. , 2000; Le oni e t al., 
ligggtrLuukanen et all . l200(i [Manninen et aR Il999(l . as 
shown in Fig. 1Mb). 

The issue of tunnel juncti on asymmetr y in SIN IS re- 
frigerators was addressed bv lPekola eTdl l|2000H) . This 
effect is fortunately weak: these authors theoretically 
showed that the maximum cooling power is reduced by 
7% in the case the junction resistances differ by a fac- 
tor of two, as compared to a symmetric structure with 
the same total junction area. Furthermore, the reduc- 
tion is only about 25% even when the resistance ratio 
is four. This effect stems from the "self aligning" char- 
acter of the double junction structure: the voltage drop 
across each junction is simultaneously close to A/e, thus 
corresponding to the maximum cooling power, when the 
voltage across the whole SINIS system is close to 2 A/e. 
This fact is due to the high non-linearity of the current- 
voltage characteristics of the two junctions which carry 
the same current. The experiments have confirmed such 
a weak depe ndence of the cooling power on the structure 
asymmetry iPek ola et al\ . \206Oti(\ . 

In the low-temperature regime the situation is rather 
different. While power load from electron-phonon inter- 
action becomes less and less dominant by decreasing the 
temperature, and typically below 0.1 K the cooler behav- 
ior can be described as if the lattice would not exist at 
all, other factors can limit the achievable cooling power 
as well as the lowest attainable electron temperature. 
First of all, at the lowest temperatures the maximum 
achievable cooling power of a NIS junction is intrinsically 
limited, given by Q oc ie'^RT)~^A^/^{kBTe,N f/^. Then 
non-idealities in the tunnel barriers, where the possible 
presence of Andreev-like transport channels (e.g., pin- 
holes) may strongly degrade the cooling power (in the 
high temperature regime this contribut ion is in general 
overc ome by thermal activation) (Ba rdas and AverinL 
Furthermore, nonequilibrium effects in the N re- 
gion as well in the S electrodes may be a limitation. 
In the N region, a suppression of the cooling power is 
expected by increasing t he quasiparticle relaxation time 
l)Frank and Krechlll997(l . i.e., by driving the electron gas 
far from equilibrium. In the superconductors, a non- 
thermal distribution can alter the cooling response of the 
refrigeration process as well as the presence of hot quasi- 
particle excitations (like with large-area junctions) may 
be responsible of additional heat load into the N region. 
Yet, quasiparticle states within the superconducting gap 
represent a crucial problem at the lowest temperatures, 
yielding anomalous heating in the N regi on and limit- 
ing th e achievable minimum temperature ( Pekola et all 
Such quasiparticle states, due mainly to inelastic 



scattering in the superconductor l)Dynes et all Il984|) or 
to inverse proximity effect from the nearby N region, are 
generally taken into account by adding an imaginary part 
(r) to the superconductor DOS, as in Eq. (|12|) . Figure 
I28f c') shows a representative set of cooling curves taken 
in an AI/AI2O3/CU SINIS refrigerator at different lat- 
tice temperatures, where at the lowest temperature (red 
dash-dotted line) (this typically happens around or below 
To « 100 mK) the electron gas first undergoes heating, 
then it is strongly cooled at the expected bias around 
Vrefr — 2 A/e. A similar anomaly w as reported in some 
experiments on SINI S refrig e rators llFisher et all Il999t 
IPekola et all l20n4ai l2000bl: ISavin et all 120011) . The 
anomalous heating was attributed to the presence of such 
states within the S gap that give rise to a sort of dissi- 
pative chan nel which dominates h e at current in a certain 
bias range l|Pekola et a/.L l2004a(l . IPekola et al\ l|2004a|) 
theoretically demonstrated that the minimum achiev- 
able electron temperature (Tmin) in SINIS refrigerators is 
set by the amount of quasiparticle states present within 
the superconducting gap, and is given by Tmin/Tc — 
2.5(r/A)^/^ at Vrefr — 2A/e in quasiequilibrium (see 
Fig. EHI^d))- We note that a measure of the F/A value in 
real NIS junctions is approximately given by the ratio of 
the zero-bias to the normal state junction conductance 
at low temperature. The existence of quasiparticle states 
within the superconducting gap thus sets a fundamental 
limit to the minimum achievable temperature, and great 
care has to be devoted to get rid of their presence to opti- 
mize the refrigeration process at the lowest temperatures. 



2. SilS2(ISi) structures 

The enhancement of superconductivity b y quasipar- 
ticles extraction was proposed in 1961 by iParmenteil 
l)l96lf) in the context of S1IS2 tunnel junctions, where 
Si and S2 represent s uperconducto r s with different en- 
ergy gaps. Later on, iMelton et al\ l)l98l[l theoretically 
discussed the possibility to exploit such a system to re- 
alize a refrigerator (addressing both the basic features 
and performance). F rom an experimental point of view, 
IChi and ClarkS (|l97fl) observed in Al films an enhance- 
ment of the energy gap of t he order of 40% due to quasi- 
particle extraction. Then, iBlamire et al ] (I1991D . usmg 
Nb/A10a;/Al/A10a;/Nb structures, were able to obtain 
an enhancement of the critical temperature of the Al 
layer by more than 100%. The physical mechanism giving 
rise to this effect w as dis cussed by Hcslinga and Klaowiiy 
(1993. ) and bv fZaitsevI <jl992|) . More recently, also 
iNevirkovetsI l)l997|) observed in similar structures an en- 
hancement of the Al gap by quasiparticle extraction. 

Figure I29f a) shows the calculated heat current ver- 
sus bias voltage for a S1IS2 tunnel junction at T = 
Te,si = Te.s2 = 0.2 Ai(0)/A:b for various A2/A1 ratios 
(see Eq. (|30|) '). The quantity Q{V) is an even function of 
V, thus allowing connection of two junctions in a sym- 
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FIG. 29 SiIS2(ISi) refrigerator, (a) Calculated cooling power 
Q of a S1IS2 junction vs bias voltage V a.t T = Te,si = 
Te,S2 ~ 0.2 Ai(0)/fcB for several A2/A1 ratios. Red-dotted 
line represents Q when S2 is in the normal state, (b) Mea- 
sured cooling of a Ti island to and in the superconducting 
state by quasiparticle tunneling. _Ro is related to the electron 
temperature Te,s2 in Ti (see text). Dashed lines: Ti is in the 
normal state. Thin solid lines: Ti is cooled from the normal 
to the superconducting state. Thick solid line: Ti is already 
superconducting at Kefr = 0. (c) Measured minimum elec- 
tron temperature (Tmin) of Ti versus bath temperature To- 
(b) and (c) are adapted from iMannincn et a/.. .1999i) . 



FIG. 30 Semiconductor-superconductor structures for elec- 
tronic cooling, (a) Optical micrograph of a typical Al/Si/Al 
microrefrigerator and schematics of the measurement setup. 
The current /q and voltage Vth are used to determine the elec- 
tron temperature in the n"'"^-Si layer, and the bias voltage V 
is used to change its value, (b) Cross section of the structure 
along the line AB in (a), (c) Measured electron temperature 
Te^AT in Si versus bias voltage across the Al/Si /Al cooler at 
differ ent substrate temperatures. Adapted from iSavin et all 
I2OO3D. 



metric configuration as in the NIS case, and is positive 
for \V\ < (Ai(r) -f A2(r))/e where hot quasiparticle ex- 
citations are removed from S2. Moreover, the heat cur- 
rent is maximized at = zt ( Ai (T) — ^21^) ) ; where 
it is logarithmically diverging ifHarrisl [l974t IXinkhaiil 
Il996j) (note that this is somewhat broadened by the 
smeari ng in the density of st a tes in a realisti c situ- 
ation llFrank and Kre'dl Il997t iManninen et all Il999t 
iPekola et all 1200481) 1. From Fig. it follows that 

heat extraction from S2 only occurs if ^2{T) < Ai(T) 
holds. Then, atV ^ ±(Ai(r) + A2(T))/e a sharp transi- 
tion brings Q{V) to negative values (more details about 
the heat transpor t in S i IS 2 junctions can be found in 
l|Frank and Krechl Il997|| ). The dotted line represents 
Q{V) when the electrode S2 is in the normal state (i.e., 
NIS case). Notably, when both electrodes are in the su- 
perconducting state, Q(V) can exceed significantly that 
in the normal state. This peculiar characteristic of heat 
transport in S1IS2 junctions makes them promising in 
realizing a "cascade" refrigerator that might operate at 
bath temperatures higher than those for a SINIS cooler, 
where several combined superconducting stages are used 
to efficiently cool a normal or a superconducting region. 

The experimental observation of electron cooli n g in a 
superconductor was reported bv IManninen et aZI l|l999() 
using Al/insulator/Ti tunnel junctions. In this experi- 
ment, aluminum (i.e, the larger-gap superconductor) was 
used to cool the electron system of a Ti strip from 1.02 Tc2 
to below 0.7 Tc2, where Tc2 = 0.51 K was the Ti criti- 



cal temperature. Figure l^^ b) shows a representative set 
of measurements taken at different bath temperatures 
Tq. In particular, it shows the zero-bias resistance of 
the thermometer junctions i?o against the bias voltage 
across the refrigerators (T^cfr)- Dashed lines indicate the 
behavior at bath temperatures where Ti is still in the nor- 
mal state: the electron temperature in Ti decreases (i.e., 
Rq increases) by increasing Vrcfr, reaching a minimum 
slightly below Vrch — 2Ai/e ~ 420 fieV, as expected for 
a SINIS refrigerator (Ai is the energy gap of Al). Note 
that dashed lines are the upside-down equivalent, for ex- 
ample, of the measurements shown in Fig. I27r c). Thin 
solid lines indicate the temperature behavior when Ti 
is cooled from the normal to the superconducting state. 
This happens for Tq larger than Tc2 but below 625 mK, 
at T^efr values for which Rq has a deep minimum. Here, 
for instance, starting from Tq = 520 mK, the minimum 
electron temperature reached in Ti was Tmin ~ 320 ± 40 
mK, thus demonstrating the effectiveness of the cool- 
ing mechanism also in all-superconducting refrigerators. 
Thick solid line corresponds to the case where the Ti 
strip is already superconducting, i.e., at a bath temper- 
ature below Tc2- The results for the minimum electron 
temperature reached in Ti against Tq is represented by 
the open circles in Fig. I29f c). Comparison to the the- 
ory gave good agreement with the experiment (lines in 
Fig. I29r c)) with an electron-phonon interaction constant 
Sxi w 1.3 nW/zim^K^. 
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3. SSmS structures 

Superconductor-semiconductor structures have re- 
cently attracted interest in the fie l d of electro n cooling 
llBuonomo et "aIl. l2003tlSavin et a^J . l2003ll200l|) . The ba- 
sic idea is to exploit, for the active part of the cooler, 
heavily-doped semiconductor layers instead of normal 
metals, with the natural Schottky barrier forming at the 
contact to the metal electrode s llBreim ani Il999t iLiithL 
l2nnit iRhoderick and Wil1ia.Tn4 ll98St ISzeL ll98lL ll98fJl . 
The peculiar characteristic intrinsic to this scheme stems 
from the possibility to alter up to a large extent the semi- 
conductor electronic properties (like, for example, the 
charge density, etc.) and to change the transmissivity 
of the Schottky barrier through proper doping or choice 
of the semiconducting layer (see also Sec. VI. F. 2). 

Heavily-doped silicon is a natural choice both from a 
practical and technological point of view, mainly owing 
to the widespread and well-developed Si technology. Up 
to now, the only evidence of superconducting electron 
cooling with semiconductors comes from microrefrigera- 
tors realized with the silicqn-on-insulator (SOI) te chnol- 
ogy l|Buonomo et a;J . l2003lSavin et a]] . l2003ll200l|) . The 
first demonstration of electron refrigeration in Si was re- 
ported bv .Savin et al. (2001), where the authors obtained 
a maximum temperature reduction under hot quasiparti- 
cle extraction of the order of 30% at the bath temperature 
To = 175 mK. Figure lSUT a') shows the optical micrograph 
of a typical Al/Si/Al cooler. The two bigger Al electrodes 
are used for cooling, while the two smaller for detecting 
the electronic temperature in Si. A sketch of the struc- 
ture cross-section is displayed in Fig. I30f b). where the 
n~''~'~-Si region appears just on top of the silicon diox- 
ide insulating layer. The structures were characterized 
by a doping level Nn = 4 x 10"'^^ cm~^. A set of cool- 
ing curves ext racted from one of such devices is shown 
in Fig. EOJc) l|Savin et alV l2003ft . It displays the elec- 
tron temperature Tejv in the n"'""'"-SOI layer versus bias 
voltage across the refrigerator at different starting bath 
temperatures. The observed electronic temperature re- 
duction is of the order of 60% at Tq « 150 mK. Although 
the device volume and the cooler resistances were rather 
high, the authors attributed this significant cooling to a 
very small electron-phonon coupling constant in Si: as a 
matter of fac t, the latter was de termined to be Sgi ~ 0.1 
nW/^m^K'^ feavin et ad l200l[l . i.e., about one order of 
magnitude smaller than in Cu. The maximum achieved 
cooling power was Qmax ~ 1-3 pW at Tq = 175 niK, 
mainly limited by the high value of the specific contact 
resistances {Rc ~ 67.5 kfi/xm^). These authors also ad- 
dressed the e ffect of carrier conc entration in Si on cooling 
performance (|Savin et oi.Ll2003i) . 

Buonomo et al. l)Buonomo et all 1200,'^ also obtained 
a cooling effect in Al/Si/Al refrigerators, with a device 
configuration similar to that of Fig. I30r a) . Their struc- 
tures were characterized by a doping level Nd = 8 x 10^^ 
cm~'^ and contact specific resistances Rc ~ 100 kil/im^. 
Notably, they confirmed the same value reported by 



Savin et al. (2001) for the electron-phonon coupling con- 
stant in Si. In addition, these authors also explored the 
Nb/Si/Nb combination for electron cooling. Their re- 
sults, however, showed no cooling effect, owing proba- 
bly to a too transmissive interface at the contact with 
Si. This effect was ascribed both to the slightly lower 
Schottky b arrier height of Nb/Si contact (~ 0.5. ..0.6 eV 
ijChang et a l. 1971; Heslinga and Klapwii k, 1989)) with 
respe ct to the Al/Si (~ 0.6. ..0.7 eV |chang et all 
Il97lt ISinghl . 11994')) an d to diso rder-enhanced subgap 
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, 1199,1) . Fur- 
ther results on the electron-phonon coupling constant 
as well as on the electroni c thermal conductivity in n- 
typc Si were reported in IjHeslinga and KlanwiikL Il992t 
Kivinen et al., 2003). 

Because of the low electron-phonon coupling con- 
stant in heavily doped Si, one may expect to cool is- 
lands with larger volumes than in the normal metal 
case, and perhaps improve the high-temperature limit 
of cooling by some factor. The issue of high Rc val- 
ues that limit the maximum achievable cooling power 
still represents a drawback of these refrigerators that 
has to be overcome. High doping of the semiconduc- 
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3 Il985f) . see also Sec. VI.F.2) as well as engi 



neering of the Schottky barrier height through incor- 
poration o f interface l^^y ^ rs at the r netal-semico n ductor 
junction (ICant ile et al', 1994a b; Costa et al, 199U 

De Franceschi et al 1998a^b . 2000; Grant and \y aldrQj, 

1988l:lKovanagi et a/.LIl99a lMarincUi eral[.l200q) ( 



can be 

exploited in order to tailor the interface transmissivity. 



4. SF systems 



As discussed at the beginning of Sec. IV.G.ll decreasing 
the contact resistance (Rt) is not a viable route to en- 
hance heat current in NIS junctions. A possible scheme 
to surmount the problem of Andreev reflection at the 
metal-superconductor interface is to exploit, instead of 
an insulating barrier, a thin ferromagnetic (F) la yer in 
good electric contact with S jGiazotto ad 12002^. The 
physical origin of the SF cooler operation stems from the 
spin-band splitting characteristic of a ferromagnet. The 
electron involved in the Andreev reflection and its phase- 
matched hole belong to opposite spin bands; as a conse- 
quence, depending on the degree of spin polarization V 
of the F layer, strong suppressio n of th e Andreev current 
may occur at the SF interface l|de Jon g and Bccnakkej 
Il995j) . In the limit of large V, the subgap current is 
thus strongly suppressed while efficient hot-carrier trans- 
fer leads to a sizeable heat current across the system. In 
the following we give the main results of such a proposal. 

The impact of partial spin polarization {V < 100%) in 
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FIG. 31 Superconductor-ferromagnet refrigerator, (a) Cal- 
culated heat current Q of a SF junction vs bias voltage V 
at T = Te,F = Te,s = 0.4 A{0)/kB for several spin polar- 
izations v. The inset shows the same quantity calculated at 
the optimal bias voltage against temperature for some values 
of v. (b) Calculated maximum cooling power surface den- 
sity Qa versus V for various temperatures. The inset shows 
the coefficient of performance rj calculated at the optimal bias 
voltage versus temperature for some V values. Adapted from 
llGiazotto et a;.ll200l I2OO5II . 



with a N region (i.e., a SFN refrigerator) it turned out 
that cooling effects comparable to the SF case can be 
achieved with a F-layer thickness of a few nm (i.e., of the 
order of the magnetic lenght). We note that the large 
Qa typical of the SF combination makes it promising as 
a possible higher-temperature first stage in cascade cool- 
ing (for instance, over or around 1 K), where it could 
dominate over the large thermal coupling to the lattice 
characteristic for such temperatures. 

The results given above point out the necessity 
of strongly spin-polarized ferromagnets for a proper 
operation of the SF refrigerator. Among the pre- 
dicted half-metallic candidates it is possible to indi- 
cate Cr 02 (Brcncr et at, 2000; Kamocr et at, 198^ 
ISchwartzl Il98fit) . for which valu es of V in the 
range 85. ..100% have been reported llCoev et a/.l flQQSt 
iDedkov et ai\ . 120021: Iji al\ . l200lt iParker et a Jl2002t) 
(Coi- ^Fe.x)S2 fMazin', "2000), NiMnSb f de Groot et all 
Il983l). Sr^FeMoOfi ( Kob avashl et a/., . 1998j) and NiMnV9 
llWeht and PickettL 11999^ So far no experimental real- 
izations of SF structures for cooling applications have 
been reported. 



the F region is displayed in Fig. l3ir a') where the heat cur- 
rent Q versus bias voltage across the junction is plotted 
for some values of P at T = Te,F = Te,s = 0.4 A(0)/fcB 
{Te^F is the electron temperature in F). For each V 
value there exists an optimum bias voltage (Vopt) which 
maximizes Q. In the limit of a half-metal fcrromag- 
net ( i.e., V = 100%) (Cocv and Sanvito, 2004; Mazin, 
Il999j) this value is around Vopt — A/e. Moreover, for 

V = 94% there is still a positive Q across the junc- 
tion. The inset of Fig. |^ (a) shows Q calculated at 
each optimized bias voltage against temperature for var- 
ious values of V. The heat current is maximized around 
T = Topt ~ 0.4 A(0)/fcB, rapidly decreasing both at 
higher and lower temperatures. 

The junction specific cooling power is shown in Fig. 
I31f b'l. where Qa (evaluated at each optimal bias volt- 
age) is plotted versus V at different bath tempera- 
tures. Notably, for a half-metallic ferromagnet at T = 
0.4 A{0)/kB, cooling power surface densities as high as 
600 nW/^m^ can be achieved, i.e., about a factor of 30 
larger than those achievable in NIS junctions at the op- 
timized interface transmissivity (i.e., T ~ 3 x 10^^ at 
T ~ 0.3A(0)/fcB, see Fig. gjb)). This marked difference 
stems from SF specific normal-state contact resistances 
as low as some 10^^ fl^m^ tha t are currently achieve d 
in highly-transmissive junctions l|Upadhvav et aZ.lll998|) . 
The inset of Fig. I31f b) shows the junction coefficient of 
performance (77) calculated at the optimal bias voltage 
versus temperature for various V values. Notably, for 

V = 100%, 7] reaches ~ 23% around T = 0.3 A(0)/A:b 
and exceeds 10% for V = 96%. In addition, in light of 
a possible exploitation of this structure in combination 



5. HTc NIS systems 

Heat transport in high-critical temperature (HT ^) NIS 
junctions was theoretically addressed bv iDevvatov et al\ 
(2000). In these systems the cooling power depends on 
interface transmissivity as well as on orientation of the 
superconductor crystal axes and temperature. In par- 
ticular, these authors showed that the maximum pos- 
itive heat current in these structures can be achieved 
in junctions with zero superconducting crystallographic 
angle, at temperature T = 0.45 Tc and bias voltage 
V « 0.8 A(0)/e. The behavior of Q{V) turns out to 
be qualitatively similar to that of NIS junctions based 
on low-critical temperature superconductors (see Fig. 
I26f a)). and with comparable values (in relative units). 
From this it follows that the cooling power of electronic 
refrigerators based on HTc materials is approximately 
two orders of magnitude larger than in NIS junctions 
based on low-critical temperature superconductors (at 
much lower temperatures). 

A somewhat different cooling effect in H T^ supercon- 
ductors was predicted bv'Svidzi nskvl l)200^ . who showed 
that an adiabatic increase of the supercurrent in a ring 
(or cylinder) made from a HTc superconductor may lead 
to a cooling effect. The maximum cooling occurs if the 
supercurrent is equal to its critical value. For a clean 
HTc superconductor, the minimum achievable tempera- 
ture (Tmin) was found to be around Tmin = /Tc, with 
To the initial temperature of the ring, thus meaning that 
substantial coolin g can be ac hieved using large Tc values. 

Experimentally, iFeel |l993l) realized a Peltier refrigera- 
tor junction exploiting a HTc superconductor and oper- 
ating around liquid nitrogen temperatures. In particular, 
his device consisted of a BiSb alloy and YBa2Cu307_5 
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superconducting rods connected by a small copper plate 
which acted as the cold junction of the device. The lat- 
ter showed a maximum cooling of 5.35 K below the bath 
temperature Tq = 79 K. The figure of merit of the junc- 
tion, Z, was estimated to be as large as 2.0 x 10^'^ K^^. 



6. Application of (SI)NIS structures to lattice refrigeration 

One well established appUcatio n of SINIS refrig- 
erators concerns lattice cooling llClark et all . l2005t 
iLuukanen et all |2p00'; Manninc n et aZ.lll997l) . As a mat- 
ter of fact, while NIS tunneling directly cools the elec- 
tron gas of the normal electrode, the phonon system 
can be refrigerated through the electron-phonon coupling 
(see Eq. (HDl). The latter, however, is typically very 
small at the lowest temperatures, thus limiting the heat 
transfer from the surroundings to the electrons. This 
situation normally happens whenever the metal to be 
cooled is in direct co ntact with a thick substrate, for 
instance, oxidized Si l|Leivo et all 1199(1 iNahum et a.l\ . 
^94) : only the electrons of the N region cool down while 
the metal lattice presumably remains at the substrate 
temperature. The metal lattice can be refrigerated con- 
siderably if the thermal resistance between the phonons 
and the substrate is not negligible compared to that be- 
tween the electrons and the phonons. One effective choice 
to meet this requirement, that was i ndeed suggested at 
the beginning of cooling experiments ijFisher et al. ',19951 
|Nahum et al, 199^^ as well as for detecto r applications 
iDciker et al, 2004* Doricsc et al, 2004; Irwin et^aJl 
Il99(k .Nakm and Martinis, . .1995; .Pekola et all l2004bt 
lUllom et a,ll \2004\ . is to exploit a thermally isolated 
thin dielectric membrane on which the N region of the 
cooler is extended. In this way, tunneling through the 
NIS junction will cool down the electrons of the metal, 
then the phonons of the metal (via electron-phonon 



brane phonons ( 


Clark et a/.Ll2005HLuukanen et allhOOd 


iManninen et al. 


. Il997l) according to 



QsiNIs{y',Te,N,Te,s) + Qph-sub{Tph,Tsub) — 0, (81) 



where Tg 



N 



Tsub, Tph is the lattice tem- 



perature in the dielectric membrane, Tsub is the lattice 
temperature in the substrate, and Qph-sub is the rate 
of exchanged energy between the membrane phonons 
and substrate phonons. Eventually, if additional de- 
vices are standing on the same dielectric platform (for 
instance, detectors, etc.), the latter will cool down first 
the phonons of the device and then its electrons through 
the electron-phonon interaction. 

Dielectric membranes made of silicon nitride (Si3N4) 
have proved to be attractive for this purpose in light 
of their superior thermal isolation properties at low 
temperatures. Low-temperature heat transport char- 
acterization as well as thermal relaxation in low-stress 
Si3N4 m embranes and films were quite recen t ly ad - 
dressed ()Holmes et all Il998t iLeivo and Pekolal Il998|) . 
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FIG. 32 Lattice refrigeration with SINIS. (a) SEM image 
of a Si3N4 membrane (in the center) with self-suspended 
bridges. Two normal-metal cold fingers extending onto the 
membrane are used to cool down the dielectric platform. The 
AI/AI2O3/CU SINIS coolers are on the bulk (far down and 
top) and the thermometer stands in the middle of the mem- 
brane, (b) Maximum temperature decrease (Tmin/To) versus 
bath temperature To measured at different positions in the 
cooler shown in (a), (c) Maximum lattice temperature de- 
crease on the membrane versus bath temperature measured 
in two other samples similar to that shown in (a). In this 
case, the refrigerators exploited many small-area junctions 
arranged in parallel in a comb-like configuration, (d) SEM 
micrograph of a NIS refrigerator device with a neutron trans- 
mutation doped (NTD) Ge resist ance thermometer atta ched 
on to p of it. (c) is adapt ed from iLuukanen et all\200(t) . (d) 
from JClark et all 12005^ . 



The first demon stration reported of lattice cooling 
ijManninen et all ^97) exploited such membranes in 
combination with AI/AI2O3/CU SINIS refrigerators. In 
this experiment the authors were able to achieve a 2% 
temperature decrease in the membrane at bath temper- 
atures To « 200 mK. 

Figure I32r a) shows a SEM image of a typical new- 
generation lattice cooler fabricated on a Si3N4 membrane 
with self-suspended bridges. The membrane consists of 
a low-stress Si3N4 film deposited by low pressure chem- 
ical vapor deposition (LPCVD) on Si, and subsequently 
etched (with both wet and dry etching) in order to create 
the suspended bridge structure. In such thin membranes 
phonon propagation is essentially two-dimensional. The 
condensation of the phonon gas into lower dimen- 
sions i n ultrathin membranes was also theoretical l y dis- 
cussed fAnghel and ManninenLll999HAnghel et aZ.l.ll998t 
iKuhn et al, 2004). The self-suspended bridges improve 
thermal isolatio n of the dielect ric platform from the heat 
bath (,Leivo and Pekolal Il998|) . The image also shows 
the AI/AI2O3/CU refrigerators of the SINIS type that 
are placed on the bulk substrate (i.e., outside the mem- 
brane) to ensure good thermal contact with the bath. 
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The N cold fingers extend onto the sihcon nitride mem- 
brane, whose temperature is determined through an ad- 
ditional SINIS thermometer placed in the middle of the 
structure. 

The lattice refrigeration effect achieved in this SINIS 
refrigerator is shown in Fig. I32f b). Here the maximum 
temperature decrease of the membrane (Tmin/To) against 
bath temperature Tq is displayed (red circles), and shows 
that temperature reduction as high as about 12% was 
reached in the 400. ..500 mK range. The electron refrig- 
eration effect in the Cu region was also measured at two 
different positions in the device, i.e., nearby the cooling 
junctions and on the membrane. Notably, the Tmin/To 
behavior is almost the same for the different sets of data; 
this basically means that: a) good thermalization was 
achieved in the cold fingers; b) the electron-phonon cou- 
pling was sufficiently large while Kapitza resistance be- 
tween Cu and Si3N4 was sufficiently small to ensure the 
lattice temperature on the membrane to be nearly equal 
to the Cu electron temperature on the membrane itself. 
The best results of lattice refrigeration by SINIS coolers 
reported to date are shown in Fig.|22Jc) ( Luukan cn et all 
I2OOO) for two other devices (labeled C and D in the fig- 
ure) similar to that of Fig.|23a). These devices exploited 
three Cu cold fingers and several small-area junctions ar- 
ranged in parallel in a comb-like configuration for the 
SINIS cooling stage. The junction specific resistances 
were Rc — 1.39 kil/im^ and Rc = 220 ^Ifim^ for device 
C and D, respectively. Lattice temperature reduction as 
high as 50% at 200 mK was achieved in the sample with 
lower Rc, thus confirming the effectiveness of small-area 
junctions in yielding larger temperature reductions. The 
achieved cooling power in these devices was estimated 
on the pW level. The reduction of the refrigeration ef- 
fect at the lowest temperatures can be explained in terms 
of larger decoupling of electrons and phonons, but also 
the effects discussed in Sec. IV. C"T1 should play a role. 

Figure 021 (d) demonstrates the realization of a 
complete re frigerator device including a thermometer 
l)Clark et a,ll 12005^ . where four pairs of NIS junctions 
are used to cool down a 450 x 450 /Ltm^ suspended Si3N4 
dielectric membrane. Each NIS junction area is 25 x 
15 /im^, and the N electrode consists of Al doped with Mn 
to suppress superconductivity while Al is used for the su- 
perconducting reservoirs. Also shown is a neutron trans- 
mutation doped (NTD) Ge resistance thermometer (with 
volume 250 x 250 x 250 ^m"^) glued on the membrane. In 
such a refrigerator the authors measured with the NTD 
Ge sensor a minimum final temperature slightly below 
240 mK starting from a bath temperature Tq — 320 mK, 
under optimal bias voltage across the cooling junctions. 
This result is promising in light of a realistic implementa- 
tion of superconducting refrigerators, and shows the pos- 
sibility of cooling efficiently the whole content of diel ec- 
tric membranes through NIS junctions l)Pekolal|2005(l . 

Possible strategies to attain enhanced lattice refriger- 
ation performance in the low temperature regime (i.e., 
below 500 mK) could be a careful optimization in terms 



of the number and specific resistance of the cooling junc- 
tions as well as to exploit superconductors with the gap 
larger than in Al. Making the dielectric platform thin- 
ner and reducing thermal conduction along it should 
also increase the temperature drop across the mem- 
brane. The exploitation of the described method around 
or above 1 K, however, is still now not so straightfor- 
ward, mainly due to the strong temperature dependence 
of the electron-phonon interaction that thermally shunts 
more effectively the N electrode portion standing on the 
bulk substrate to the lattice (note that also the thermal 
conduction along the niembrane is strongly temperatu re 
dependent ijKuhn et al l l2004l:lLeivo and Pekola[ll998|) 'l. 
Toward this end, a reduction of the N volume placed out 
of the membrane should help; in addition, SiIS2(ISi) re- 
frigerators (see Sec. IV.C.2|I as well as SF junctions (see 
Sec. IV.C.4|I might significantly improve the niembrane 
cooling in the higher temperature regime. 



7. Josephson transistors 

A further interesting field of application of SI- 
NIS structures concerns supe r condu cting weak links 
l|Golubov et al!.l . l2004^^Likharevlll979^ . In particular, in 
diffusive SNS junctions, i.e., where the junction length 
largely exceeds the elastic mean free path, coherent se- 
quential Andreev scattering between the superconduct- 
ing electrodes may give rise t o a continuum spectrum of 
resonant levels (jBelzig et all Il999l: iHelkkila et all l2002t 
IyIp. 1998) responsible for carrying the Josephson cur- 
rent across the structure. The supercurrent turns out 
to be given by this spectrum weighted by the occupa- 
tion number of correlated electron-hole pairs, the lat- 
ter being determined by the quasiparticle energy dis- 
tribution in the N region of the weak link. In con- 
trollable Josephson junctions, the supercurrent is mod- 
ulated by dr iving the quas i partic l e distribution out of 



HVolkov 


and Pavlovskii. 
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^\kov and Takavanagi', '1997'; 'van Wee s et all 11991 
Wilhelm et al., 1998; Yip, 1998 ) via dissipative current 
injection in the weak link from additional normal-metal 
terminals. This operation principle leads to a con- 
trolled supercurrent supp ression and was successfully ex- 
ploited both in a l l-metal i j Ba selmans et al. 2002a, 199£ 
[2001albl. l20 02bl; i Huang eTaL. .2002; .Morpurgo et al 
[199^ Ishai khaida roT'era/T l2000fl (where a transi- 
tion to a TT-state was also reported) as well as 
in hybrid semiconduct o r-superconduct or weak links 



llKutchinskv et all Il999t iNeurohr et all Il999t iRichteil 
l200rt 'Sch aoers et all l2008ai 119981 l2003b^ . The situa- 
tion drastically changes if we allow current injection from 
additional s uperconducting terminals arranged in a SI- 
NIS fa shion jBaselmansl l2002t iGiazotto etal] . l2004albl 
I2OO3IJ) . In this way, thanks to the SINIS junctions, crit- 
ical supercurrent can be strongly increased as well as 
steeply suppressed with respect to equilibrium, leading 



43 




FIG. 33 Josephson transistor, (a) Simplified scheme of a 
SINIS-controlled Josephson transistor. The Josephson cur- 
rent in the SjNSj weak link (along the white dashed line) 
is controlled by applying a bias Vsinis across the SINIS line 
connected to the weak link, allowing to increase or decrease 
its magnitude with respect to equilibrium, (b) Theoretical 
normalized critical current Ij versus Vsinis calculated in the 
quasiequilibrium limit for several lattice temperatures To for 
a Nb/Cu/Nb long junction. The inset shows the supercurrent 
vs electron temperature characteristic calculated ai (p — -k 12. 
(c) SEM image of an Al/Cu/Al Josephson junction includ- 
ing the AI/AI2O3/CU symmetric SINIS electron cooler. The 
SNS long weak link is placed in the middle of the structure. 
Also shown is a scheme of the measurement circuit, (d) Mea- 
sured critical current 7j versus Vsinis at three different bath 
temperatures To for the device shown in (c). The inset dis- 
plays the measured differential curre nt gain Gi versus /sini s 
at To = 72 mK. (b) is adapted from (Giazo tto tt ttZ.ll2003l3) . 
while (c) and (d) from dSavin et a/., .2004) . 

to a tunable structure with large current and power gain. 

A simplified scheme of this class of transistors is dis- 
played in Fig. |33Ja). A diffusive SjNSj long Josephson 
junction of length Lj (i.e., with Lj ^ where ^0 is 
the S,] coherence length) shares the N region of a SI- 
NIS control hne of length Lsinis- The superconductors 
Sj and S can be in general different (i.e., with gaps Aj 
and A5, respectively); in addition, the Sj electrodes are 
kept at zero potential, while the SINIS line is biased with 
l^SiNis- The supercurrent (/j) behavior in response to a 
bias voltage Vsinis stems from the degree of nonequilib- 
rium the SIN IS line is able to generate in the weak link 
acco rding to l)Heikkila et alV l2002t IWilhelm etdV Il998t 

IjiVsims) I dEjsiE, 0)(1 - 2fiE, ^sinis)), 

e^N Jo 

(82) 

where Rn is the weak link normal-state resistance, 
is the phase difference across the superconductors, and 
is{E,4') is the spectral supercurre nt, obtainable from 
the solution of the Usadel equations ijUsadej 11970(1 . The 



most straightforward situation occurs at sufficiently low 
lattice temperatures, i.e., typically below 1 K, when 
ie-e < -^siNis < f^e-ph (sec Scc. II). In such a case, strong 
electron-electron relaxation forces the electron system 
in the N region to retain a local thermal quasiequilib- 
rium (see also Sec. IV.dip . so that 1 — 2/(£', Vsinis) — 
tanh[£^/(2fcBre,Af(VsiNis))]- The transistor effect in the 
structure thus simply depends on the electron tempera- 
ture AT (Vsinis) established in the weak link upon bi- 
asing the SINIS hne according to Eq. The full cal- 
culation of the behavior of a prototype Nb/Cu/Nb long 
Josephson junction with integrated Al/ Al2 0.'^/Cu SINIS 
electron cooler is displayed in Fig. ESfb) ( Giazotto et all 
l2003b). Here the normalized /,/ is plotted against Tasinis 
for several bath temperatures Tq. For all chosen val- 
ues of To, the supercurrent value increases monotonically 
up to about Vsinis ~ 1.8 As/e, where the SINIS hne 
provides the largest cooling power and allows to attain 
the lowest electron temperature. Then, by further in- 
creasing the voltage, an efhcient suppression of Ij occurs 
due to electron heating: the device behaves as a tun- 
able superconducting junction. The above given results 
can be easily understood recalling that for a long SNS 
junction at low lattice temperature (i.e., fc^To <C Aj) 
and for ksT^.N ^ Eth = f^D/Lj {D is the N-region 
diffusion coefficient), Ij depends exponentiall y on the ef- 
fective electron temperatu re T^^n ifwilhelm et a/.l . [l997t 
IZaikin and Zharkovlll98l|) and is almost independent of 
Tq. Hence, long junctions are more appropriate for the 
device to obtain large Ij changes with small T^^n varia- 
tions (as indicated by the red hatched region in the inset 
of Fig. EHb)). In the short junction limit {Lj <C ^o), 
conversely, the Ij temperature dependence is set by the 
energy gap Aj, thus implying a much reduced effect from 
the cooling line. Furthermore, it is expected that the 
power dissipation values from two to four order of mag- 
nitude smaller than with all-normal control lines can be 
obtained in these structures, thus making them promis- 
ing as mesoscopic transistors for low dissipation cryogenic 
applications. 

So far, the only successful demonstra tion of this 
trans istor-like operation was reported by l|Savin et all . 
l2004j) in Al/Cu/Al SNS junctions integrated with 
AI/AI2O3/CU SINIS refrigerators. The SEM image of 
one of these samples along with a scheme of the mea- 
surement setup is shown in Fig. I^^f c). The structure pa- 
rameters were the following: cooler junction resistances 
Rt — 240 il, Josephson weak link normal state resis- 
tance Rn = 11.5 n, and minimum SNS interelectrode 
separation Lj ~ 0.4 fxia that compared with the su- 
perconducting coherence length (^0 ~ 62 nm) provides 
the frame of the long junction regime. The critical cur- 
rent Ij versus Vsinis ^^t three different Tq is shown in 
Fig. ISST d). For each bath temperature, Ij increases 
around Vsinis ~ l-8As/e, as expected from the reduc- 
tion of Tg jv by electron cooling, while it is steeply sup- 
pressed at larger bias voltages. The resemblance of the 
experiment with the curves of Fig. |2SIb) is evident. In 
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the present case, Ij enhancement under hot quasipar- 
ticle extraction by more than a factor of two was ob- 
served at To = 283 mK. The transistor current gain 
G/ = dIj/dlsiNis, shown in the inset of Fig. ISST d). 
obtained values in the range —11... 20 depending on the 
control bias. As far as power dissipation is concerned, 
these authors reported low dissipated power on the 10"^'^ 
W level in the extraction regime while of some tens of pW 
in the regime of Ij suppression. 

The transistor behavior for arbitrary inelastic scat- 
tering strength in the SINIS line (or, equivalently, for 
arbitrary values of isiNis) as well as for different SNS 
junction l engths , was the oretically addressed in detail in 
(iGiazotto et all l2004allJ l. The role of geometry, ma- 
terials combination, phase dependence, and the input 
noise power were also discussed. Notably, a marked 
supercurrent transition to a 7r-state under nonequilib- 
rium (about two times larger than that achievable with 
an all- normal control channel (Wil helm et a/. . .1998: yTdI 
11998)) was predicted to occur for negligible or moderate 
electron-electron interaction. We recall that in the tt- 
state (Bulaevskii et ai, 19"^ the supercurrent flows in 
opposite direction with respect to the phase difference 
(j) between the two superconductors, i.e., such a sign re- 
versal is equivalent to the addition of a phase factor tt 
to the Josephson current-phase relation. Furthermore, 
current gain in the range 10^... 10^ and power gain up to 
lO'^ where predicted to occur depending on the control 
voltage. 

Finally, transistor operation in a Josephson tunnel 
junction integrated with a S1IS2IS1 refrigerator was also 
theoretically addressed in the quasiequilibrium regime 
IGiazotto and Pekola, ,200.5,). In this case, the device 
benefits from the sharp characteristics due to the pres- 
ence of superconductors with unequal energy gaps (see 
also IV.(r2)l . that leads to improved overall characteris- 
tics as compared to the SINIS-controUed SNS junction in 
the same transport regime. 



D. Perspective types of refrigerators 

Any electric current that is accompanied by the extrac- 
tion of hot electrons (or holes) can be used, in principle, 
for refrigeration purposes. For instance, this may happen 
in thermionic transport over a potential barrier as well as 
in energy-dependent tunneling through a barrier. Since 
an exhaustive analysis of all possible predicted refrigera- 
tion methods is far beyond the limits set to this Review, 
in the following we give a brief description of a few exam- 
ples that we believe are relevant in the present context. 
In such devices, several different effects may contribute 
to the refrigeration process (e.g., thermionic transport, 
quantum tunneling as well as thermoelectric effects), so 
that making a proper " classification" of the refrigeration 
principle is, strictly speaking, rather difficult. As a con- 
sequence, we shall try mainly to follow the definitions as 
they were introduced in the original literature. 



1. Thermionic refrigerators 

A vacuum thermionic cooling device consists of two 
electrodes separated by a vacuum gap. Cooling oc- 
curs when highly energetic electrons overcome the vac- 
uum barrier through thermionic emission, thus reduc- 
ing the electron temperature of one of the two elec- 
trodes. In such a situation, the refrigerator operation 
is mainly affected and limited by radiative heat transfer 
between the electrodes. The thermionic emission cur- 
rent density (jn) is given by the Richardson equation, 
Jb. = AqT"^ exp(— ^^), where $ and T are the work func- 
tion and the temperature of the emitting electrode, re- 
spectively, Aq = ATTemkg / is the Richardson constant, 
and m is the electron mass. From the expression above, 
a strong redu ction o f jji is expected upon lowering the 
temperature. iMahan ( 1994) developed a simple model 
for thermionic refrigeration, and demonstrated that its 
efficiency can be as high as 80% of the Carnot value. Vac- 
uum thermionic refrigerators are generally characterized 
by a higher efficiency as compared to thermoelectric cool- 
ers, and are considered to be an attractive so lution for 
future refrigeration devices l)Nolas and Goldsm id. 1999). 
However, the high values typical of currently available 
materials make thermionic cooling efficient, at present, 
only above 500 K. 

Several ideas on how to increase the cathode emis- 
sion current and to improve the o peration of th ese re- 
frigerators have been propos ed (Hishinuma et al!, 200^ 
l2002l:lKorotkov and LikhareVt.l999;.PurceU et aL . .1994tl . 
Korotkov and Likharev suggested to cover the emitter 
with a thin layer of a wide-gap semiconductor, and to 
exp loit the resonant emission cur rent to cool the emit- 
ter ijKorotkov and Likharevl Il999() . The analysis of such 
a thermionic cooler predicted an efficient refrigeration 
down to 10 K. 

The issue of high $ values can be overcome by a re- 
duction of the distance between the electrodes (in the 
submicron range) and by the applica tion of a strong 
electri c field. Following this scheme, iHishinuma et all 
l)200lj) theoretically analyzed a thermionic cooler where 
the two electrodes where separated by a distance in the 
nanometer range. In such a situation, the potential 
barrier is essentially lowered (see Fig. I34r a)). allowing 
both thermionic emission and energy-dependent tunnel- 
ing. As a consequence, rather small (if compared to vac- 
uum thermionic devices) external voltages (^ 1...3 V) are 
required in order to produce significant electric currents. 
For suitable values of the applied voltage and distances 
between the electrodes, electrons above the Fermi level 
dominate the electric transport (both thermionic emis- 
sion over the barrier and tunneling through the barrier), 
thus leading to cooling of the emitter. As it can be in- 
ferred from Fig. I34r b'). the contribution of the energy- 
dependent tunneling to total cooling is essential, and this 
refrigerator could be classified as a vacuum tunneling 
device. The cooling power surface density in this com- 
bined thermionic-tunneling refrigerator was predicted to 
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FIG. 34 Refrigeration by combined tunneling and thermionic 
emission, (a) Schematic diagram of the potential barrier pro- 
file V{x) for $ = 1 eV and with electrode separation of 60 
A. (b) Heat flow distrib ution at T = 300 K. Adapted from 
(iHishinuma et aZ.L 1200 ll^. 



obtain values as high as 100 W/cm^ at room tempera- 
ture. However, in the only experimental demonstration 
of this device, a moderate emission current (below 10 
nA) was reported at room temperature , with an observed 
temperature reduction of about 1 mK ijHishinuma et cd\ . 
|2p03). So far, no experimental demonstration of vacuum 
thermionic refrigeration at cryogenic temperatures has 
been reported. 



2. Application of low-dimensional systems to electronic 
refrigeration 

The exploitation of low-dimensional systems gives ad- 
ditional degrees of freedom in order to engineer ma- 
terials that may lead to enhanced operati on o f ther- 
moelectric and thermionic de v ices ( Hicks et all 11995 ; 
Hishinuma et all 120021: ISalesl l2002t iSofo and MahaE , 
19941) . Some of the limitations intrinsic to vac- 



uum thermionic refrigerator s can be overcome with 
solid state thermionic coolers l|Mahan and Woodslll998t 
IShakouri and Bowersl Il997j) . As a matter of fact, 
modern growth techniques easily allow one to con- 
trol both the barrier height and its width within a 
wide range of values. One disadvantage of solid state 
thermionic coolers stems from the thermal conductiv- 
ity of the barrier which is essentially absent in vacuum 
devices. Nevertheless, a large temperature reduction 
can b e achieved by usin g a m u ltilayered heterostruc- 
ture llMahan and Woo d^ llQQSt IShakouri and Bowersl 
Il997t IZhou ad |l999j). According to theor y 
l|Mahan and Woods!. 11998: "Sha kouri and Bowe^ Il997(l . 
heterostructure-based thermionic refrigerators perform 
somewhat better as compared to thermoelectric coolers. 
The predicted temperature reduction of a single-stage de- 
vice at room temperature can be as high as 40 K, and 
this value can be significantly increased in a multilayer 
configuration. So far, howev er, the experimental im ple- 
mentations of single barrier llShakouri et all Il999^ and 
superlattice l)Fan et all l200lt IZhang et all \200^ refrig- 
erators have reported temperature reductions of a few 
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FIG. 35 Quantum-dot refrigerator, (a) Scheme of a quantum- 
dot refrigerator, (b) Energy-level diagram of the structure. 
The reservoir R is cooled as its quasiparticle distribution 
function is sharpened by resonant tunneling through quan- 
tum dots Dl and Dr to the electrodes Vl and Vr. From 
lEdwards et al... .l99^^ . 



degrees at room temperature. 



and thermionic cooling (p 
IShakouri and Bowerl Il997j) 



Further improvement of thermionic refrigeration can, 
in principle, be achieved by combining laser c ooling 

jMal'shukov a,nd Chanl 

In such an opto-thermionic 
device, hot electrons and holes extracted through 
thermionic emission lose their energy by emitting pho- 
tons rather than by heating the lattice. The the- 
oretical investigation of a GaAs/AlGaAs-based opto- 
thermionic refrigerator predicted specific cooling power 



densities of the order of s everal 
(iMal'shukov and ChaoLl20(il . 



W/cm-^ at 300 K 



The presence of singularities in the energy spectrum 
of low-dimensional systems can be used in solid state re- 
frigeration. For instance, the discrete energy spectrum 
in quantum dots can b e exploited for refrigerat ion at 
cryogenic temperatures l|Edwards et adll995l !T9931. In 
such a quantum-dot refrigerator (QDR) (see Fig.|33a)) 
a reservoir (R) is coupled to two electrodes via quan- 
tum dots (Dl and Dr ) whose energy levels can be tuned 
through capacitively-coupled electrodes. The QD energy 
levels can be adjusted so that resonant tunneling to the 
electrode Vl is used to deplete the states in R above 
and, similarly, holes below in R tunnel to Vr, (see 
Fig. E3b)). As a consequence, the net result will be to 
sharpen the quasiparticle distribution function in R, thus 
leading to electron refrigeration. In spite of a rather mod- 
erate achievable cooling power, the QDR was predicted 
to be effective for cooling electrons of a micrometer-sized 
two-dimensional electron gas reservoir at mK tempera- 
tures, and even of a macroscopic reser voir at lower tem- 
peratures llEdwards et oil Il99fl IT991 . 
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VI. DEVICE FABRICATION 

A. Structure typologies and material considerations 

This section is devoted to the description of the main 
techniques and experimental procedures used for the fab- 
rication of typical superconducting electronic refrigera- 
tors and detectors. Owing to the great advances reached 
in the l ast decades in micro- and na nofabrication tech- 
nology ('Bhushanl 12004 iTimd . Il999fl . the amount of in- 
formation related to fabrication methods is too large to 
be covered here and beyond the scope of the present re- 
view. Therefore, we only briefly highlight all those issues 
that we believe are strictly relevant for this research field. 
In particular we first of all focus on the two typologies 
of existing superconducting structures, namely all-metal 
and hybrid devices, and on the main differences between 
them, both in terms of materials and fabrication tech- 
niques. The former concern structures where the active 
parts of the device, i.e., both the superconducting ele- 
ments and the normal regions are made of metals; in 
hybrid structures, the non superconducting active part 
of the device is made of doped semiconducting layers. 

As far as all-metal-based structures are concerned, 
they are realized with low-critical-temperature thin-film 
superconductors (normally Al, Nb, Ti and Mo), while the 
normal regions usually consist of Cu, Ag or Au. Their 
fabrication protocol includes patterning of a suitable 
radiation-sensitive masking layer through electron-beam 
or optical lithography in combinati on with a sh adow- 
mask (angle) evaporation technique ijPolanL FlQTTyi . The 
final device is thus realized in a single step in the de- 
position chamber, where additional tunnel barriers be- 
tween different regions of the structure are in-situ created 
by suitable oxidation of the metallic layers. Although 
all this leads to an efficient way for fabricating metallic 
structures, however the electronic properties typical of 
metals are only weakly dependent on their growth con- 
ditions and on the specific employed technique. As a 
consequence, it is hard to tailor the metallic properties 
in order to finally match some specific requirements. 

On the other hand, the situation is rather different 
with semiconductors that offer some advantages in com- 
parison to metals. The large magnitude difference of 
Fermi wave-vector between metals and semiconductors 
allows in general to observe quantum effects in struc- 
tures much bigger than with metals. In addition, the 
availability of several techniques for growing high-purity 
crystals (e.g., molecular beam epitaxy) yields the capac- 
ity to tailor the semiconductor electronic properties and, 
at the same time, enables the fabrication of structures 
characterized by coherent transport especially in reduced 
dimensionality (^Cao asso. 1990,') . The exploitation of low- 
dimensional electron systems as active elements of elec- 
tronic coolers was predicted to be th e next possible break- 
through in this research fie ld (Edwa^^^I], ^951^93 
Hicks and Drcssclhaus, 1993a b: IHicks et ad . ll993[IT99a 
Koea et aL . .2000t . .l99&.1999.) . Moreover, the realization 



of structures where charge carriers experience arbitrarily 
chosen effective potentials is a f urther advantage of mod - 
ern engineered heterostructures llYu and Cardonall20ni^ . 
Last but not least, the possibility of changing the carrier 
density through electrostatic gating allows to strongly al- 
ter some semiconductor parameters (like the mobility as 
well as the coherence length), thus enabling to have ac- 
cess to dif ferent electronic transport r egimes in the same 
structure ijFerrv and Goodnickl Il997j) . 

Like all-metal devices, hybrid structures generally ex- 
ploit the same superconducting materials. However, only 
heavily-doped Si layers in combination with supercon- 
ductors were exploited up to n ow for the real i zation 
of supercon d ucting refrig erators l|Buonomo et all l2003t 
ISavin et al[ l2003l l200l|) . The fabrication procedure 
of hybrid structures involves typically two lithographic 
steps, where the first allows to pattern the semiconduc- 
tor active layer through various techniques (such as wet 
or dry etching). The second one the patterning of the 
superconducting electrodes. Differently from all-metal 
devices, at the metal-semicon ductor interf ace forms the 
well-known Schottky barrier ijLutli l200l[l . The latter 
can be tailored through suitable doping of the semicon- 
ductor, allowing to control the interface transmissivity 
(i.e., the junction specific resistance) over several orders 
of magnitude. 

Both all-metal and hybrid structures are normally 
fabricated on semi-insulating semiconductor substrates 
which provide both the mechanical support and rigidity 
as well as the thermalization of the device at the bath 
temperature. It is noteworthy to mention that undoped 
thermally-oxidized Si substrates are widely used toward 
this end for metallic structures, while hybrid coolers typ- 
ically exploit the silicon-on-insulator (SOI) technique in 
order to provide the Si active layer. 



B. Semiconductor growth techniques 

Depending on the physical principle exploited, semi- 
conductor growth technologies can be either physical or 
chemical, such as gas and liquid-phase chemical pro- 
cesses. In the following we give a brief survey of the 
most common semiconductor growth methods. 

The lCzochralskil |)1918|) crystal pulling method is prob- 
ably the most common technique for growing semicon- 
ductor bulk single crystals and the largest amount of 
Si used in the semiconductor industry is obtained using 
this technique. This technique allows to grow large semi- 
conductor single crystals. However, it is often sufficient 
and less demanding from an economical point of view, 
to grow a thin ( ~ 1 /im) perfect crystal layer on top of 
a bulk crystal of lower quality. This is accomplished by 
growing epitaxially the top layer, i.e., the atoms forming 
the latter build a crystal with the same crystallographic 
structure and orientation as the starting substrate (in 
contrast, non-epitaxial layers can be amorphous or poly- 
crystalline). 
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Gas-phase epitaxy is a widespread technique and rep- 
resents nowadays the most important process for the in- 
dustrial fabrication of Si an d GaAs devices. In chemical 
vapour deposition (CVD) jAdamsl Il983t iGrovd Il967t 
ISzc. 1985 ). the constituents of the vapor phase chemically 
react at the substrate surface and the product of the reac- 
tion is a film deposited on the substrate. Nowadays, CVD 
is widely employed for growing several elemental semi- 
conductors and alloys and, depending on the specific ap- 
plication, it has been implemented and adapted in a num- 
ber of different configurations (Jensen, 1989). Among 
these we mention 1) metalorganic CVD (MOCVD), ex- 
ploits a thermally heated reactor and organometallic gas 
sources; 2) plasma-enhanced CVD (PECVD), exploits a 
plasma discharge to provide the energy for the reactions 
to occur; 3) atmospheric pressure CVD (APCVD), does 
not require a vacuum environment. 

While typical CVD processes are carried out in low 
vacuum (i.e., in the pressure range from 10~^ to s everal 
Torr), molecular b eam e pitaxy (MBE) (lArthuil Il96g ; 
Chang and Ludekel. Il97fit Ichdl 1197(1 ICho and Art~ . 
19751) IS an ultrahigh vacuum (UHV) growth technique, 
performed at pressures usually lower than 10"^" Torr. In 
such vacuum conditions the amount of residual gas con- 
taminants present in the growth chamber is minimized, 
thus allowing to deposit epitaxial layers of high purity 
and quality. The main features of this method are a pre- 
cise control of both chemical and stoichiometric composi- 
tion, and a perfect tuning of doping profiles on the scale 
of a single atomic layer (,Cho . 1979). The film growth 
concept is fundamentally an evolution of UHV evapora- 
tion, where thermal or electron-beam sources are used 
to create a flux of molecular species. Significant MBE 
work has been achieved with Si and Ge as well as with 
IV-IV and II- VI semiconductor compounds and several 
metals, but probably the largest amount of research has 
been devoted to III-V semiconductor alloys. 



C. Thin-film metals deposition methods 

1. Thermal evaporation 

Vacuum evaporation is one of the oldest an d sim- 
plest thin film depositio n techniques (HoUandi Il958t 
iMaissel and Gland. IT970|) . It is a physical vapor depo- 
sition (PVD) method widely used to deposit a variety 
of materials, from elemental metals to alloys and insu- 
lators. Evaporation is based on the boiling off or subli- 
mation of a heated source material onto a substrate sur- 
face. The result of the evaporant condensation is the 
final film. The rate of atoms or molecules (A/'e) lost un- 
der evaporation from a source material per unit area per 
unit time can be expressed by the Hertz-Knudsen relation 
a; = a[P^(r) -P][27rMfcBr]-i/^ where a is the evapo- 
ration coefficient {a ~ 1 for a clean evaporant surface), P 
is the ambient hydrostatic pressure acting on the evap- 
orant in the condensed phase, Py(T) is the equilibrium 



vapor pressure of the evaporant, and M its molecular 
weight. This expression shows that the evaporation rate 
strongly depends on the evaporant vapor pressure. Most 
common metals typically deposited by thermal evapo- 
ration (such as Al, Au, Ga, and In) usually have vapor 
pressures in the range between 10~^ to 1 Torr in the tem- 
perature window of 600°. ..2000° C; conversely refractory 
metals (such as Nb, Mo, Ta, W and Ft) or ceramics (such 
as BN, and AI2O3) reach such vapor pressures at much 
higher temperatures, thus making more difficult the ex- 
ploitation of this technique for the deposition. 

Usually evaporation is performed in high or ultra- 
high vacuum (in the range 10~^...10~^'' Torr), where the 
mean free path i for the evaporant species is much larger 
than the substrate-source distance. This translates in 
an almost line of sight evaporation which prevents cov- 
ering of edges perpendicular to the source, t he latter 
also referred to as the lack of step coverage ijMadoul 
Il997j) (not e that this pro perty is at the basis of lift-off 
processes llMoreauL Il988|) as well as of angle (shadow) 
evaporation technique using suspended masks ijDolanL 
^77; Dolan and Dunsmuir, 1988)). Furthermore, vac- 
uum evaporation is a low-energy process (implying a neg- 
ligible damage to the substrate surface), where the typ- 
ical energy of the evaporant material impinging on the 
substrate is of the order of 0.1 eV. Nevertheless, radia- 
tive heating can be high. 

Resistive heating and electron-beam deposition arc the 
two most common methods of evaporation. The for- 
mer relies on direct thermal heating to evaporate the 
source material. This method is fairly simple, robust and 
economic but suffers from a limited maximum achiev- 
able temperature (of the order of 1800° C), which pre- 
vents the evaporation of refractory metals and several 
oxides. On the other side, electron-beam evaporation 
represents a crucial improvement over resistive heating. 
This method exploits a high-energy electron beam that is 
focused through a magnetic field on a localized region of 
the source material. A wide range of materials (includ- 
ing refractory metals and a wide choice of oxides) can 
be deposited owing to the generation of high tempera- 
tures (in excess of 3000° C) over a restricted area. Its 
main drawback relies on the generation of X-rays from 
the high- voltage electron beam which ma y damage sensi- 
tive subst rates (such as semiconductors) (jMoreaul Il988t 
ISzel Il985l) . Achievable deposition rates are up to several 
hundreds A/sec (e.g., 0.5 /zm/niin for Al) (,Madou,. 1997) . 



2. Sputter deposition 

The sputtering process has been known and used 
f or over 150 years iChapman and Mangano.. 1988t 
IChapmanL Il980t iRossnagej ll998t~IWasa and Havakawal 
Il992|) . It is a PVD method widely used nowadays for 
many applications, both in the electronic and mechanic 
industry fields as well as in the pure research environ- 
ment. This process is based on the removal of material 
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from a solid target through its bombardment caused by 
incident positive ions emitted from a (rare) gas glow dis- 
charge. The transferred momentum of the ions leads to 
the expulsion of atoms from the target material, thereby 
enabling the deposition (condensation) of a film on the 
substrate surface. Sputter deposition is generally per- 
formed at energies in the range of 0.4 to 3 keV. Further- 
more, the average energy of emitted ions from the tar- 
get source is in the range 10. ..100 eV. At these energies 
bombarding ions can penetrate up to two atomic layers 
in the substrate thus leading to a gr eat improvement of 
the adhesion of the sputte r ed fil m ijMaissel and GlaiiS 
Il97(i IWasa and Havakawal Il992|) . Normally, relatively 
high pressures (from 10"'* to 10~^ Torr) are maintained 
in the growth chamber during deposition. At these pres- 
sures the mean free path is short (of the order of 1 mm 
at 10^^ Torr) so that the material atoms reach the sub- 
strate surface with random incident angles. As a conse- 
quence, a very good step coverage can be achieved. Being 
essentially mechanical in nature, sputtering successfully 
allows the deposition of refractory metals (superconduc- 
tors) like Nb, NbN, Ta, Mo and W at temperatures well 
below their melting points. 



D. Thin film insulators 

The roles of thin film insulators in solid state elec- 
tronics are various. In particular, deposited films are 
often used as interlevel dielectrics for metals, to realize 
lithographic masking for diffusion and implantation pro- 
cesses, as well as for passivation and prot e ctive layer s 
llGhandhi hmA iNicollian and BrewsL Il98.ll: |Szi Il98,'j) . 
In addition they can be exploited as thin amorphous 
membranes on which micro- and nanostructured de- 



vices are realized ( 


Clark et oil l2005HFisher et all\19m 


Irwin et all Il996l 


Lantine: et all l2005l iLuukanen et al 


20001 iNahum and Martinisl Il995|) in light of their spe- 



Il998t iLeoni et all l2003t iManninen a;J . Il997(l . In the 

following we discuss those insulators which are consid- 
ered particularly relevant for microelectronic fabrication 
processing, i.e., silicon dioxide and silicon nitride. 

Silicon dioxide (Si02) is one of the most exploited insu- 
lators in micro- and nanoelectronics based on Si, mainly 
due to the high quahty of the Si02/Si interface. Si02 
films can be grown on Si substrates by thermal oxida- 
tion using oxygen or steam. Thermal oxidation of Si 
is generally carried out in reactors at temperatures be- 
tween 900° C and 1200° C. The resulting Si02 film is 
amorphous and characterized by good uniformity, lack 
of porosity and very good adhesion to the substrate. 
Some typical parameters of thermally grown silicon diox- 
ide at 1000° C are a refractive index of 1.46, a break- 
down stren gth larger than 10^ V/cm , and a density of 



2.2 g/cm^ jNgu ven. 1988; |Szi Il985j) . An alternative 
way to deposit s i licon dioxide laye rs is through CVD tech- 
niques l|NguvenL ITgsl iReil Il99qi . In particular PECVD 



is considered an effective techn ique l|Adam l ll98fiHHessl 
Il984t iKaganowicz et ai\ . Il984|) because of the low depo- 
sition temperature (Si02 is generally deposited in the 
temperature range 200°... 500° C). Silicon oxide can be 
deposited from sil ane (Sill4) with O2, CO2, N2O or 
CO llAdams et all I198H iHollahanl. \W74\ . In addition, 
plasma oxide film properties are strongly dependent on 
the growth conditions such as reactor configuration, RF 
power, frequency, substrate temperature, pressure and 
gas fluxes. Typical parameters for plasma-deposited sil- 
icon dioxide films at 450° C are a refractive index of 
1.44. ..1.50, a breakdown strength of 2. ..8 x 10^ V/cm, 
and a density of 2.1 g/cm^. 

Silicon nitride is another insulating film that forms 
good interfaces with Si. It is nowadays su ccessfully used 
as interlevel d ielectric (ISwan et J . 1196^ . in multilayer 
resist systems ijSuzuki et aZ.lll982j) . as well as a protective 
coating as it provides an efficient barrier against moisture 
and alkah i ons (e.g Na) fSzc, 1985). PECVD is com- 
monly used ijAdamsl . 1986; Hess, 1984) because of the low 
deposition temperature (250°... 400° C) implying low me- 
chanical stress. Silicon nitride is typically formed by re- 
acting silane and ammonia (NII3) or nitrogen in the glow 
discharge. As for Si02, the properties of the final film 
stron gly depend o n the deposition condi tions llAdams 
1986HChow et all Il982t iDun et all Il98lt iNeuven et al 



1984j) . Typical parameters for plasma-deposited silicon 
nitride films at 300° C are a refractive index of 2.0. ..2.1, 
a density of 2. 5. ..2. 8 g/cm^, and a breakdown strength of 
6 X 10^ V/cm. 



E. Lithography and etching techniques 

Fabrication of thin film metallic circuits as well as semi- 
conductor micro- and nanodevices requires the genera- 
tion of suit able p atterns t hrough lithographic p rocesses 
fCamubell!, 1200 It iJaegeil 120021: iPlummer etd] . I2000D . 
Lithography, indeed, is the method used to transfer such 
patterns onto a substrate (e.g.. Si, GaAs, glass, etc.), thus 
defining those regions for subsequent etching removal or 
material addition. 

In photolithography ijLevinson and Arnoldl . Il997f) a 
radiation-sensitive polymeric material (called resist) is 
spun on a substrate as a thin film. The image expo- 
sure is then transferred to the resist through a pho- 
tomask, consisting normally of a glass plate having the 
desired pattern of clear and opaque areas in the form of 
a thin (~ 1000 A) Cr layer. Two types of resists can be 
used in such a process, i.e., positive and negative resist s 
l|Colclaseil fl980HMoreauL 119881: Irhompson et all\l994} . 
In the former, the solubility of the exposed areas in a 
solvent called developer is enhanced, while in the latter 
the solubility is decreased. After exposure, the resist is 
developed and reproduces the desired pattern images for 
the subsequent processing. The radiation source for pho- 
tolithography depends on the desired final resolution, al- 
though the latter is mainly limited by effects due to light 
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diffraction ijLevinson and ArnoldL Il997j) . In particular, 
high pressure Hg lamps (with a wavelength A = 365 nm 
i-line or A = 436 nm g-line) allow a line-width larger 
than 0.250 /im, while in the range 130. ..250 nm deep UV 
(DUV) sources like excimcr lasers are necessary, such as 
KrF (A = 248 nm) and ArF (A = 193 nm). 

Electron-beam lithography ijBrewerl Il980t 
iMcCord and Rookgl 11997^) represents an attractive 
technique for the fa b ricatio n of micro- and nanos- 
tructures ijKern et all Il984j) . This method exploits 
a focused electron beam (with energy in the range 
10. ..100 keV and diameter of 0.2. ..100 nm) to ex- 
pose a polymer-based electron-sensitive resist (such 
as polymethyl-methacrylate (PMMA)). Like in pho- 
tolithography, the resist can be either positive or 
negative. Even if the wavelength of the impinging 
radiation beam can be smaller than 0.1 nm, the 
maximum achievable resolution is set by the electron 
scattering in the resist and backscat tering from substrate 
(known as t h e "pr o ximity effect" llHoward et al. . 11983 ; 
Jackel et all 1984; 'Jamoto and Shimizu nT983t Kvsei , 
19831: iKvser and Viswanathan. _1975)), so that the 



and recipes are available ijCotler and EltaL Il990|) . 



resolution is generally larger than 10 nm (note, however, 
that resolutions as high as 2 nm have been achieved on 
some materials (Mochcl et al., 1983)). 

In addition to lithographic procedures, etching of thin 
films or bulk substrates represents an important step for 
the fabrication of the final structure llMadodlT997Hs^ 
119851) . Toward this end, insulating or conducting thin 
films are exploited as masking layers for subsequent ma- 
terial removal. Two crucial parameters of any etching 
process are directionality and selectivity. The former 
refers to the etch profile under the masking layer. In 
particular, for an isotropic etch, the etching rate is ap- 
proximately the same in all directions, leading to a spher- 
ical profile under the mask. With anisotropic etch, the 
etching rate depends on the specific direction (e.g., a par- 
ticular crystallographic plane) thus leading to straight 
profiles and sidewalls. Selectivity instead represents how 
well the etchant can differentiate between the masking 
layer and the layer that has to be removed. Moreover, 
etching te chniques can be divid ed into wet llGhandhi , 
' l983HKen dall and Shoultz. [l997l) and dry llMadoil 11997 : 
Wasa and Havakawai il992i|) categories. In wet etching, 
the substrate is placed in a liquid solution, usually a 
strong base or acid. The advantage of wet etching stems 
from its higher selectivity in comparison to dry methods. 
Wet etching is in general isotropic for most substrates, 
and various solutions that yield anisotropic etching are 
available for some materials. In dry etching, the sub- 
strate is exposed to a plasma in a reactor where ions can 
etch the substrate surface. The great advantage of dry 
etching with respect to wet etching resides in its higher 
anisotropy (that allows vertical etch walls), and smaller 
undercut (that enables smaller lines to be patterned with 
much higher resolution). Several materials (e.g., insula- 
tors, semiconductors as well as refractory metals) can be 
successfully dry etched, for which a variety of chemistries 



F. Tunnel barriers 

1. Oxide barriers 



Superconducting tunnel junctions represent key ele- 
ments in a number of electronic applications (Solvma^ 
I1972I) that span fr o m sin gle electron transistors 
Grabert and Devoretl Il992tl to Josephson devices 



Barone and Paterndl 1198^) ?^ iust to mention two rel- 
evant examples. In addition, they are crucial 
building blo c ks of superconduct i ng mi crorefrigerators 
l)Leivo et all Il996t iNahum et all. Il994l) as well as of 
ultrasensitive microbolometers l)Castellano et Ml Il997t 
.Nahum and Martinis. 1993^) . The simplest picture of 
a tunnel junction can be given assuming a rectangu- 
lar barrier of height (f)i and width w. The electron 
transport across the barrier can be easily described 
within the Wentzel-Kramers-Brillouin (WKB ) approxi- 
mation. The main results of this analysis are ijSimmonsl 
E963b) i) an exponential dependence of the zero-bias 
junction conductance (Go) on the barrier width. Go cx: 
exp[—2w^/2m*(j)i /h] , where m* is the effective mass of 
electrons in the barrier; ii) a quadratic voltage {V) de- 
pendence of the conductance G{V); iii) a weak insulating- 
like quadratic dependence of G on the temperature T. 
In general, image forces acting on the electrons tunnel- 
ing through the barrier will red uce bot h its height and 
its effective thickness ( Simmons. Il963al) . Criteria i)-iii) 
require that the dominant process through the barrier 
is direct tunneling and can be used to extract the junc- 
tion parameters in a re alistic situation ijBrinkman et all . 
Il97nt ISimmonsl . [l963a() . Among the above given criteria, 
i) seems a necessary, but not sufficient condition to es- 
tablis h that tunneling is the dominant transport me cha- 
nism C Rabson aIl . l200lH Zhang and RabsonLl2004|) due 
to the possible presence of pinholes (i.e., small regions 
where the insulator thickness vanishes) in the barrier, 
and it is generally accepted that only criterion iii) could 
be safely used to rule out the presence of su ch pinholes 
in the barrier l).Tonsson-Akerman et all l200nt) and assess 
the junction quality (for instance, a reduction of junc- 
tion conductance of about 15% on cooling from 295 K to 
4.2 K is believed to be a good practi cal indication of high- 
quality AlOj^ barriers ^Gloos et r7i.l 12003. ,2000.} ; see also 
Sees. Ill and IV). An important figure of merit of tunnel 
contacts is the junction specific resistance Rc = RjA, 
where Rj = Gq^ is the contact resistance and A its area. 
A route to decrease Rc is to reduce w or choose materials 
with lower effective barrier height. 

Among the available barrier materials, aluminum oxide 
(AIO2:) is probably the most widespread insulator used 
to fabricate metallic tunnel junctions because it can be 
easily and reliably grown starting from an Al film. Its 
main parameters are a typical barrier height 4>i Ki 2 eV, 
although values in the range 0.1. ..8. 6 eV have been re- 
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ported llBarner and Ru ggiero'. "198^;' Gundlach and H0I4 
ll97ltlKadlec and Kadlea.1975t.Lau and Coleman. ■1981jl 
and a dielectric constant smaller than that of bulk 
AI2O3 (4.5.. .8.9 at 295 K (iBolz and TuveL Il983')1. Sev- 
eral methods are currently exploited to fabricate high- 
quality (i.e., highly uniform and pinhole- free) aluminum 



(TMatsuda et all 11999: 


Parkin et al.. 


1999; 


Sun et al. 


l200(i iTsuge and Mitsuzukal. Il997l iZhang et al. 


2001) 



oxidation m air fJVliyazaki and iezuka|. ( 1995(1 and plasma 
oxidation ( Gallag her et all 119971: iMoodera et all Il995t 
ISun et all . ll999l|) of a thin Al layer (typically below 2 nm), 
just to mention the most common techniques. The latter 
two methods lead in general to higher Rc values (of the 
order of several kf] /im^ or larger) with respect to natu- 
ral oxidation. The natural oxidation allows to achieve the 
desired Rc by simply changing the oxidation pressure and 
time (i?c also depends on the original thickness of the Al 
layer), and almost any R^ value can be produced by fol- 
lowing such a procedure. Aluminum oxide junctions with 
Rc values as low as some tens of /xm^ or lower are cur- 



rently realized w i th na t ural oxidation |Childrcss_et_cd 
2001; Dcac et al, 2004; Fuiikat a e^'aZI. f2 001: Li u et al 
2003: Parkin et al. .1999; .Sun et all I200Q; .Zhang et al 



However, low-i?c junctions are more prone to 
defects or pinholes in the barrier that dramatically de- 
grade their performance. Promising methods and mate- 
rials have been proposed for the fabrication of oxidized 
barriers with Rc as low as some f2 ^w? among which 
we mention juncti ons made of Z rAlO,; jLiu et al, 2003; 
IWang et a/.Ll2002^ and HfAlO:, llWang et oij . l2003() . 



2. Schottky barriers 

The metal-semiconductor (NSm) ju nction is an issu e 
that dates back more than 60 years llSchottkvl I1939D . 
but still nowadays represent s a releva nt topic both in the 
physics of semiconducto r s (iBrenngm. | 1999- .LiithL l200lt 



iRhoderick and WiUiamsl. 119881: ISzeL 



1981 



^85|) and in 
1987t ISinghL 



device applications I^MiU mai^ii^Grabe 
[l994i'l . The two most important types of NSm junc- 
tions are the Schottky barrier (SB), showing a rectify- 
ing diode-like I-V characteristics, and the ohmic contact, 
whose I-V is almost linear. Most of NSm junctions, with 
a few exceptions such as co ntacts with InAs , InSb and 
In2;Gai_a;As (for x > 70%) llKaiivama et a/.Lll973|) . are 
affected by the presence of the SB that drastically affect 
their electric behavior. For most semiconductors (in the 
following we concentrate on n-type semiconductors but 
similar conclusions can be given for p-type ones), the SB 
height {(psBn) has a rather weak dependence upon the 
metal used for the contact; for instance, for metal/n-Si 
contacts (f)sBn — 0.7. ..0.85 eV, while for metal/n-GaAs 
contacts (t)sBn = 0.7. ..0.9 eV (jRhoderick and Williamd 
IT988: Singh : 1994) . This fact is explained in terms of the 
Fermi- level pinning a t the NSm interface l)Bardeenlll947t 
iHeineL Il965t iMonchL IT99?il) . 



The current across a Schottky junction depends on 
a number of different mechanisms. In the limit where 
thermionic emission dominates the electric transport, the 
rectifying action of a biased NSm junction is described 
as I{V) — Is[exp{eV/kBT) — 1], where the detailed ex- 
pression for the saturation current Is depends on the 
assumptions made on carrier transport l|Brennanl [l999t 
ISze. 1981, 1985). In such a case the junction specific 
resistance is Rc oc T~^e'x.p{e4>sBn/kBT), thus meaning 
that it can be lowered mainly by decreasing (psBn (typ- 
ical Rc values at doping levels Nd < 10^^ cm~'^ for 
metal/n-Si contacts are of the order of 10^^. ..10"'^^ fi/xm^ 
almost independent of No). By contrast, tunneling 
across the SB can be the dominating transport mecha- 
nism if the semiconductor is heavily doped. In such a case 
I{V) cx exp[~a{(t)sBn - V)/^/TT^], with a = 2h-^Vm*e 
where m* is the semiconductor effective mass and e the 
dielectric permittivity, i.e., the junction is not rectify- 
ing and the current is proportional to V for small volt- 
ages. The contact is thus said to be ohmic and yields 
Rc oc exp(Q;^sBn/V-^-D)- This shows that Rc can be 
reduced up to a large extent by lowering the SB height 
and doping as heavily as possible (again, for metal/n- 
Si contacts and Nd > 10^^ cm~'^, Rc can be in the 
range 10^. ..10® fi/xm^). All this shows the advantage of 
using NSm contacts owing to the possibility of tuning 
the contact specific resistance over several orders of mag- 
nitude (from metallic-like to tunnel-like characteristics) 
through a careful choice of metal-semiconductor combi- 
nations and proper doping levels. This trick is commonly 
exploited to control the NSm interface resistance in cur- 
rent semiconductor technology, although heavy doping of 
the semicond uctor just in proximity to the metal is of- 
ten preferred iG iazotto et fl/.L l2001albl: iKastalskv et all 
I1991I: IShannonl I197&: .Taborvski et all\l99^ . 



VII. FUTURE PROSPECTS 

Low temperature solid-state cooling is still at its in- 
fancy, although operation of a number of individual 
principles and techniques have been demonstrated to 
work successfully. Yet combinations of cascaded micro- 
refrigerators over wider temperature ranges employing 
several stages, or combinations of different refrigeration 
principles, e.g., fluidic coolers together with electronic 
coolers, do not exist. In principle, compact low-power 
refrigerators could be fabricated using micro-machined 
helium-based fluidic refrigerators, for instance based on 
Joule- Thomson process (Little, 1984), and these devices 
could then be directly precooling NIS-refrigerators with 
niobium {Tc = 9 K) as a superconductor. A lot of engi- 
neering effort is, however, needed to make this approach 
work in practise as a targeted micro-refrigerator. 

The solid-state micro-circuits have already proven to 
yield new operation principles and previously unknown 
concepts have been discovered in cryogenic devices, as 
demonstrated throughout this review. We believe that 
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what was demonstrated here is just a presentation of a 
beginning of a new era in low temperature physics and 
instrumentation. As possible new classes of devices we 
could mention those utilizing thermodynamic Carnot cy- 
cles with electrons. Brownian heat engines with elec- 
trons are predicted to a chieve efficiencies close to ideal 
l|Humphrev et al\ . \2Q0'^ . It may be possible in the fu- 
ture to make use of other types of gated cycles where 
energy selective extraction of electrons is producing the 
refrigeration effect. As a conceivable example, a com- 
bination of Coulomb effects and superconducting energy 
gap could form the basis of operation of a refrigerator 
where cooling power would be proportional to the oper- 
ating frequency of the gate cycle. Such a device would 
thus be principally different from the static electronic re- 
frigerators presented in this review, where a DC bias is 
in charge of the redistribution of hot electrons. 

At low temperatures, additional relaxation channels 
besides the electron-phonon scattering, such as coupling 
between electrons and photons, become important. More 
knowledge is needed on these mechanisms. 

In Subs. IV.C.7I we describe how the non-equilibrium 
shape of the distribution function sometimes leads to im- 
proved characteristics of the device. It would be inter- 
esting to see if such effects could be employed to improve 
also the properties of the radiation detectors or other 
practical devices. 

The presently obvious application fields of electronic 
micro-refrigerators include astronomical detectors both 
in space as well as those based on the earth, materials 
characterization instrumentation, e.g., those devices em- 
ploying ultra high resolution x-ray micro-analysis, and 
security instrumentation, e.g., concealed weapon search 
on the airports. It is, however, evident that once real- 
ized in a user-friendly and economic way, refrigeration 
becomes very important in high-tech based industry in 
a much broader perspective. Low temperature electron- 
ics and superconducting devices are often characterized 
by their undeniably unique possibilities, but they are of- 
ten superior to the room temperature ones also in speed 
and power consumption. Therefore, mesoscopic on-spot 
refrigerators to attain the low temperatures forming the 
basis of these instruments are urgently needed. 
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